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Abstract 

A  field  deployable,  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  device 
fiber  coupled  to  a  pair  of  12.5”  Ritchey-Chretien  telescopes  was  used  to  study  atmo¬ 
spheric  propagation  for  open  path  lengths  of  100  m  to  1,000  m  to  estimate  atmospheric 
transmission  at  key  High  Energy  Laser  (HEL)  wavelengths.  The  potassium  (K)  ver¬ 
sion  of  the  Diode  Pumped  Alkali  Laser  (DPAL)  operates  in  between  two  of  the  sharp 
oxygen  rotational  features  in  the  pP  and  the  PQ  branches.  Initial  experiments  were 
performed  in  the  vicinity  of  molecular  oxygen  X3Y,~  to  b3T,p  electronic  transition 
lines  near  the  potassium  emission  line  at  770  nm.  More  than  50  rotational  lines  in 
the  molecular  oxygen  A-band  X3T,~  to  blT,p  transition  near  760  nm  were  observed. 
Temperatures  were  determined  from  the  Boltzmann  rotational  distribution  to  within 
1.3%  (less  than  2  K).  Oxygen  concentration  was  obtained  from  the  integrated  spec¬ 
tral  area  of  the  absorption  features  to  within  1.6%  (less  than  0.04  x  1018  molecules 
/  cm3).  Pressure  was  determined  independently  from  the  pressure-broadened  Voigt 
lineshapes  to  within  10%. 

The  water  vapor  absorption  lines  close  to  the  Cs  DPAL  emission  near  895  nm 
were  also  investigated.  Temperature,  pressure,  and  water  vapor  concentration  were 
determined  for  150  m  and  1  km  open-paths  with  statistical  errors  of  0.2%.  Comparison 
with  meteorological  instruments  yields  agreement  for  the  1  km  path  to  within  0.6% 
for  temperature,  3.7%  for  pressure  and  2.4%  for  concentration. 

A  Fourier  Transform  Infrared  (FTIR)  spectrometer  was  also  used  to  observe  the 
absorption  spectra  at  1  cm-1  resolution.  The  TDLAS  approach  achieves  a  minimum 
observable  absorbance  of  0.2%,  whereas  the  FTIR  instrument  is  almost  20  times  less 
sensitive.  Applications  include  atmospheric  characterization  for  HEL  propagation 


IV 


and  validation  of  monocular  passive  raging. 

The  impact  of  atmospheric  absorption  of  DPAL  laser  weapon  performance  was  in¬ 
vestigated  using  the  Line-by-Line  Radiative  Transfer  Model  (LBLRTM)  with  the  High 
Energy  Laser  End-to-End  Simulation  (HELEEOS).  The  TDLAS  device  was  used  to 
anchor  simulations  to  outdoor  atmospheric  open-path  collections.  The  implications  of 
different  laser  gain  cell  configurations  in  DPAL  systems  are  discussed,  including  spec¬ 
tral  lineshape  and  atmospheric  transmittance.  The  results  are  compared  to  existing 
HEL  systems. 
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I.  Introduction 


High  Energy  Laser  (HEL)  systems  for  weapons  applications  have  been  investigated 
since  the  1960’s.  For  example,  the  135  kilowatt  Tri-Service  Laser  (TSL)  and  the  Ex¬ 
perimental  Laser  Device  (XLD)  were  developed  in  the  early  1970’s.  The  megawatt 
class  Navy-funded  Mid-Infrared  Advanced  Chemical  Laser  (MIR ACL)  was  opera¬ 
tional  by  1980.  Most  recently  the  Navy  Laser  Weapon  System  (LaWS)  is  designed  to 
supplement  and  extend  the  Phalanx  Close-In  Weapons  System  (CIWS)  operational 
range.  [3] 

The  Airborne  Laser  (ABL)  mounted  a  Chemical  Oxygen-Iodine  Laser  (COIL)  on 
a  747-400F,  the  Boeing  YAL-1  (PrototYpe  Attack  Laser)  to  destroy  ballistic  missiles 
while  in  boost  phase.  The  overall  difficulty  of  the  mission  rivaled  that  of  the  task  of 
putting  man  on  the  moon,  with  mission  requirements  that  produced  many  technical 
challenges.  The  system  successfully  destroyed  a  liquid  fuel  boosting  missile  in  early 
2010.  [1]  Unfortunately  the  program  was  recently  canceled.  Attitudes  of  top-level 
officials  towards  the  device  were  clearly  summarized  during  a  House  Appropriations 
hearing  on  May  20th,  2010  where  Robert  Gates,  then  U.S.  Secretary  of  Defense,  said 
the  following  in  answer  to  a  question  from  Rep.  Todd  Tiahrt  (R-Kan.): 


“I  don't  know  anybody  at  the  Department  of  Defense,  Mr.  Tiahrt,  who 
thinks  that  this  program  should,  or  would,  ever  be  operationally  deployed...  The 
reality  is  that  you  would  need  a  laser  something  like  20  to  30  times  more 
powerful  than  the  chemical  laser  in  the  plane  right  now  to  be  able  to  get 
any  distance  from  the  launch  site  to  fire...  The  ABL  would  have  to  orbit  in¬ 
side  the  borders  of  Iran  in  order  to  be  able  to  try  and  use  its  laser  to  shoot 
down  that  missile  in  the  boost  phase... and  if  you  were  to  operationalize 
this  you  would  be  looking  at  10  to  20  Ifls,  at  a  billion- and- a-half  dollars 
apiece,  and  $100  million  a  year  to  operate.”  [22] 


Secretary  Gates  indicates  several  issues  for  why  the  program  was  cut.  First,  the 
ABL  was  originally  designed  to  contain  14  gain  modules  and  at  the  time  only  con- 
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tained  six.  Despite  successful  tests,  the  lower  power  output  reduced  the  ABL  effective 
range,  making  it  less  desirable.  With  a  lower  operational  range,  more  units  would  be 
needed  to  fulfill  mission  requirements  exacerbating  the  second  issue,  budget.  If  more 
devices  are  needed  to  complete  the  mission,  initial  costs  to  construct  a  larger  fleet 
drive  up  the  price,  and  also  add  to  the  yearly  operational  costs  incurred  thereafter. 
Even  after  great  technological  achievements,  and  overall  program  success,  the  ABL’s 
funding  was  greatly  reduced  and  renamed  the  Airborne  Laser  Test  Bed  (ALTB).  The 
aircraft  is  stored  at  Davis-Monthan  Air  Force  Base  in  Tucson  and  maintained  by  the 
309th  Aerospace  Maintenance  and  Regeneration  Group. 

Atmospheric  effects  like  molecular  absorption  and  scattering,  aerosol  absorption 
and  scattering,  diffraction,  atmospheric  turbulence,  and  thermal  blooming  can  all 
reduce  the  performance  of  a  laser  weapon  system.  For  example,  a  hydrogen  fluoride 
(HF)  chemical  laser  operates  near  2.7  /mi  -  near  deep  atmospheric  absorption  lines. 
The  deuterium  fluoride  (DF)  chemical  laser  has  an  emission  spectrum  around  3.9 
/mi  where  an  atmospheric  transmission  window  exists.  In  space,  over  vacuum  paths, 
an  HF  laser  would  suffice  and  would  have  superior  diffraction  limited  performance 
compared  to  a  DF  laser  system.  However,  as  soon  the  mission  included  atmospheric 
slant  paths,  the  HF  laser  becomes  highly  undesirable.  All  laser  system  parameters, 
and  their  respective  affects  on  transmission,  must  be  considered  on  a  system-by¬ 
system  and  mission-by-mission  basis,  the  focus  of  this  work  being  the  atmospheric 
transmission,  particularly  the  molecular  absorption  component. 

The  Diode  Pumped  Alkali  Laser  (DPAL)  attracted  great  interest  from  the  HEL 
community  since  2006  when  Krupke  patented  the  idea  [44]  of  a  three  level  laser  system 
by  pumping  the  D2  2Si/2  -  2P3/2  transition  followed  by  collisional  energy  transfer  to 
the  spin-orbit  split  2Pi/2  state,  and  lasing  on  the  D\  2Pi/2  -  2Si/2  transition.  By 
pumping  the  alkali  laser  gain  cell  with  diode  stacks,  electrical  sources  can  be  used 
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as  the  laser  source  of  power,  removing  the  need  for  chemical  safety  and  logistics 
chains.  With  a  gas-phase  gain  medium,  the  overall  beam  quality  of  a  DPAL  has  the 
opportunity  to  be  excellent,  and  also  gives  the  alternative  of  cycling  the  gas  phase 
medium  to  manage  thermal  effects.  [14]  All  three  alkali  variants,  potassium  (K), 
rubidium  (Rb),  and  cesium  (Cs),  operate  in  the  near-infrared  regime.  Good  beam 
quality,  attractive  origin  for  power,  appealing  wavelength  choices,  options  for  thermal 
control,  and  excellent  quantum  efficiency  make  DPALs  of  great  interest  to  the  high 
energy  laser  community. 

The  atmospheric  propagation  characteristics  of  these  emerging  DPAL  systems  are 
largely  unexplored.  Modeling  tools  such  as  the  HITRAN  database,  LBLRTM,  and 
HELEEOS  are  available  to  predict  line-by-line  atmospheric  transmission  assuming  the 
state  of  the  atmosphere  and  aerosol  content  are  specified.  However,  instruments  capa¬ 
ble  of  long  path  atmospheric  transmission  measurements  for  HEL  applications  are  not 
available.  In  the  presented  work  a  held  deployable  Tunable  Diode  Laser  Absorption 
Spectroscopy  (TDLAS)  apparatus  has  been  designed,  built,  and  tested.  The  system 
has  been  deployed  to  characterize  atmospheric  absorbance  for  potassium,  rubidium, 
and  cesium  variants  of  the  DPAL  system.  The  input  of  atmospheric  oxygen  and 
water  vapor  on  beam  propagation  has  been  investigated  to  benchmark  atmospheric 
transmission  and  laser  weapon  performance  models. 
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II.  Background 


The  study  of  atmospheric  transmission  is  not  a  new  one,  and  while  laser  spec¬ 
troscopy  could  not  exist  before  the  invention  of  the  laser  by  American  physicist 
Theodore  Mairnan  in  1960,  [57]  the  study  of  spectroscopy  began  because  of  atmo¬ 
spheric  absorbing  species.  The  German  optician  Joseph  Fraunhofer  observed  atmo¬ 
spheric  absorption  lines  from  a  dispersed  solar  spectrum  in  1814,  now  referred  to  as 
“Fraunhofer”  lines,  and  he  used  the  relative  line  spacings  to  compare  different  va¬ 
rieties  of  glass.  He  also  extended  Thomas  Young’s  double  slit  experiemnt  to  many 
slits,  leading  to  the  diffraction  grating,  an  optical  element  still  in  use  today  giving 
Fraunhofer  the  ability  to  measure  wavelengths  of  spectral  lines.  Thirty-three  years 
after  Fraunhofer’s  discovery  of  spectral  fingerprints,  the  German  physicist  Gustav 
Kirchhoff,  whose  topically  broad  work  is  still  relevant  to  physicists  and  electrical 
engineers  alike,  made  the  connection  between  absorbing  media  and  their  unique  spec¬ 
tra,  noting  differences  in  the  dispersed  spectra,  and  formed  a  relationship  between 
the  spectral  fingerprint  and  the  chemical  being  studied.  In  1848  the  French  physicist 
Jean  Bernard  Lon  Foucault  observed  the  first  laboratory-based  absorption  spectrum, 
a  sodium  flame  absorbing  yellow-orange  light  from  a  bright  arc-lamp  behind  it.  [91] 

Tunable  Diode  Laser  Absorption  Spectroscopy 

The  study  of  atmospheric  transmission  can  be  achieved  very  effectively,  not  sur¬ 
prisingly,  by  using  lasers  as  the  interrogating  source.  Tunable  Diode  Laser  Absorp¬ 
tion  Spectroscopy  (TDLAS)  relies  on  a  narrowband  tunable  laser  source  that  passes 
through  the  sample  and  is  collected  by  a  broadband  detector.  This  is  the  opposite  of 
conventional  (and  the  original)  method  for  absorption  spectroscopy,  where  an  inco¬ 
herent  broadband  source  passes  through  the  sample  and  is  collected  with  a  spectrally 
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selective  device,  such  as  a  spectrometer  or  interferometer.  For  spectrometers  or  inter¬ 
ferometers,  the  spectral  resolution  is  limited  by  the  resolving  power  of  the  dispersive 
device  used,  the  grating  type  of  the  spectrometer  or  the  maximum  path  difference 
(MPD)  of  an  interferometer.  A  tunable  laser  can  investigate  even  beyond  the  Doppler 
limit  with  a  relatively  small  physical  setup,  and  compared  to  the  other  conventional 
systems  listed  above,  has  a  very  high  spectral  power  density  because  tunable  laser 
sources  have  all  their  emission  power  distributed  over  a  small  frequency  band.  Tun¬ 
able  laser  sources  are  available  from  the  ultraviolet  (UV)  to  the  infrared  (IR)  bands 
using  a  variety  of  pump  sources,  and  there  are  many  tunable  lasers  available  that  im¬ 
plement  diodes  as  the  pump  source,  adding  further  convenience  and  removing  many 
of  the  problems  associated  with  more  cumbersome  systems  such  as  dye  lasers.  Thus, 
TDLAS  devices  have  no  need  for  large  monochromators  or  long  path  interferometers 
at  the  detector,  have  a  high  spectral  power  density,  and  have  a  excellent  detection 
sensitivity  with  spectral  resolution  as  narrow  as  the  laser  emission  lineshape.  [20] 
Investigating  atmospheric  paths  with  TDLAS  devices  is  ideally  performed  by  inter¬ 
rogating  the  open-path  to  be  investigated. 

TDLAS  Theory 

TDLAS  relics  on  Beer’s  Law,  formally  the  Beer-Lambert-Bouger  law,  [58]  which 
describes  how  light  is  absorbed  in  terms  of  the  initial  light  J0  entering  a  medium  to 
the  light  that  exits,  A, 

h(i>)  =  Io(i>)  ■  e-^-R-N  (1) 

where  Jo,  in  intensity  or  power,  passes  through  absorber  path  length  R  at  some 
wavenumber,  z>,  with  a  frequency  dependent  absorber  species  cross-section,  cr(z>) ,  and 
the  state  species  number  concentration,  N .  Species  cross-section  is  related  to  the 
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degeneracies,  g\  and  r/2  and  the  absorption  lineshape,  g(v), 

<?i2  (?)  =  —^-Mig{<>)  (2) 

8vr 

The  absorption  cross-section  and  the  number  concentration  can  be  combined  to  define 
the  absorption  coefficient,  a(z>),  [13]  that  yields  another  form  of  Beer’s  law, 

h{i>)  =  IoW  ■  e-°W-R  (3) 

where  the  absorption  coefficient,  a.{y)  is  the  product  of  the  number  concentration  and 
the  absorption  cross-section.  With  the  appropriate  a  priori  knowledge,  each  of  the 
different  variables  can  be  estimated  for  data  collection.  For  example,  if  the  absorption 
cross-section  and  concentration  are  known  the  path  length  can  be  calculated  [38];  if 
the  absorption  cross-section  and  path  length  are  known,  the  concentration  can  be 
estimated;  and  even  estimates  for  cross-section  can  be  made  this  way. 

For  atmospheric  paths,  absorbing  species  are  broadened  by  approximately  760 
Torr  (1  atm)  of  pressure  making  spectral  features  much  wider  than  commercially 
available  tunable  lasers.  Using  a  pressure  broadening  rule  of  thumb  of  3-10  MHz  per 
Torr,  [13]  a  sea  level  experiment  at  760  Torr  gives  a  pressure  broadened  lineshape  of 
roughly  2.3  GHz,  thousands  of  times  larger  than  most  commercially  available  tunable 
diode  lasers  with  line  widths  500  kHz  or  less. 

The  lineshape  function,  g{y  —  uq),  is  a  convolution  of  the  Lorentzian-shaped  ho¬ 
mogeneous  pressure  broadening  and  Gaussian-shaped  inhomogeneous  Doppler  broad¬ 
ening  effects.  The  Lorentzian-shaped  pressure  broadening  [13]  is  represented  by, 

— _ —  \  __  Aul/(2tt) 

U  U0)  (Az>L/2)2  +  (z>-h0)2 
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where  Av^  is  the  Full  Width  Half  Maximum  (FWHM)  of  the  Lorentzian  and  val¬ 
ues  can  be  calculated  from  parameters  in  the  HITRAN  database  on  a  line-by-line 
basis.  [73]  Using  the  HITRAN  parameters  7 ajr  and  7 seif,  the  air-broadened  and 
self-broadened  half-widths,  respectively,  both  at  reference  temperature  296  K,  are 
corrected  for  the  simulated  temperature  using  rw,  the  temperature-dependence  cor¬ 
rection  factor,  all  used  to  find  the  pressure-broadened  width, 


Av^  2  •  Patm 


j/ air  0.8  +  7 self  ■  0.2] 


(5) 


The  Gaussian-shaped  Doppler  component  can  be  represented  as, 
f  ~  ~ \  2  /My 

9D\y  -  U))  =  -r-^V - exp 

Aud  V  7T 

where  the  Doppler  width  can  be  calculated  from  species  mass  and  temperature  de¬ 
scribed  in  [13]  as, 


—4  ln(2) 


v  -v0 


Av 


D 


(6) 


Avd  =  2  •  v0 


2kT  ln(2) 


■ me * 


and  can  be  represented  to  give  a  wavenumber  result  as, 


Avd  =  7.2a;10  7z>0y  ^ 


(7) 


(8) 


where  T  is  in  K,  M  is  in  atomic  mass  units,  and  z>0  is  in  cm-1. 

The  Doppler  broadening  results  from  the  Doppler  shift  caused  by  the  velocity 
distribution  of  the  absorbing  (or  emitting)  atoms.  This  distribution  gives  the  overall 
width  of  the  Doppler  broadening  component.  The  convolution  of  a  Lorentzian  and 
Gaussian  produce  the  Voigt  distribution,  the  overall  lineshapes  seen  in  molecular  ab¬ 
sorption  features  at  typical  atmospheric  pressures.  Depending  on  the  case,  a  Voigt 
distribution  can  be  purely  homogeneous  or  inhomogeneous.  The  higher  the  temper- 
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ature  or  the  higher  the  pressure  the  species  experiences,  the  wider  the  Doppler  or 
Lorentzian  components,  giving  a  wider  overall  lineshape. 

Use  of  Tunable  Diode  Lasers 

Of  the  published  articles  describing  applied  TDLAS  in  the  held,  nearly  all  these 
make  the  use  of  White  or  Herriott  cells  with  vacuums  and  pumps  moving  samples 
of  the  atmosphere  into  the  cell  for  investigation.  [43,  84,  83]  Some  specifically  use 
White  cells  [10,  41,  90]  while  others  specifically  use  Herriott  cells  [31,  32,  40,  49,  50, 
85,  88,  89].  These  absorption  cells  are  small  chambers  filled  with  the  gas  sample 
under  investigation  with  opposing  mirrors  where  the  tunable  laser  is  propagated  back 
and  forth  to  give  a  long  path  length.  Sometimes  the  absorption  cell  is  simply  a  set 
distance  with  a  retro  rehector  opposite  the  laser  source.  [104]  After  a  measurement 
is  complete,  the  cell  is  then  evacuated  and  filled  with  the  next  sample.  One  of  the 
primary  advantages  of  absorption  cells  is  that  the  pressure  inside  the  cell  can  be  low¬ 
ered  to  reduce  pressure  broadening,  thus  making  line  shapes  narrower  and  removing 
ambiguity  for  spectral  features  that  may  otherwise  not  be  spectrally  isolated.  Ab¬ 
sorption  cells  also  offer  advantages  such  as  a  known  path  length,  known  mirror  losses, 
compact  size,  and  temperature  and  pressure  control.  Again,  these  advantages  do  not 
represent  atmospheric  transmission  of  a  HEL,  and  an  experiment  characterizing  HEL 
atmospheric  propagation  using  an  absorption  cell  would  not  be  complete. 

The  few  published  articles  that  demonstrate  open-path  TDLAS  are  concerned 
with  very  short  path  lengths,  usually  over  a  few  meters.  None  of  these  are  concerned 
with  HEL  propagation.  For  example,  in  [11]  an  open-path  TDLAS  system  is  used  to 
monitor  CO2  around  1.58  /mi  (6329  cm-1)  and  could  scan  nearly  a  3  cm- 1  band widt h 
including  at  minimum  two  absorption  lines  with  the  laser  diode  at  a  single  tempera¬ 
ture.  This  paper  used  wavelength  modulation  spectroscopy  (WMS)  at  70  kHz  along 


with  the  absorption  feature’s  second  harmonic.  The  emission  wavelength  was  con¬ 
trolled  using  a  modulation  of  laser  diode  current.  This  is  by  far  the  most  common 
form  of  TDLAS  performed  in  recent  papers  and  represents  a  significant  shift  from 
changing  cavity  length  to  control  the  output  wavelength  of  a  laser  source.  In  ref¬ 
erence  [11]  the  laser  power  was  split  into  three  beams  where  one  served  as  a  power 
reference  detector,  a  second  was  propagated  through  an  absorption  cell  with  a  known 
concentration  of  the  species  being  investigated  for  spectral  calibration,  and  a  third 
was  propagated  through  the  atmosphere  and  returned  to  a  third  detector  for  the  ab¬ 
sorption  signal.  This  setup  is  typical  of  current  open-path  TDLAS  systems.  Another 
open-path  TDLAS  system  described  in  [11]  gives,  at  the  time,  the  longest  open-path 
of  any  TDLAS  system  with  a  path  length  of  88  m.  Although  the  paper  stated  that 
SNR  decreased  with  a  longer  path  length,  the  authors  did  not  state  why.  With  longer 
path  lengths,  SNR  should  increase  according  to  Beer’s  Law  because  there  is  more 
absorbing  species  in  the  cell,  but  here  the  reduced  SNR  could  be  due  to  effects  from 
turbulence,  lower  power  collected  by  the  receiver  because  of  beam  spread,  etc.  The 
general  system  setup  can  be  seen  in  [55]  and  is  very  similar  to  the  open-path  TDLAS 
system  described  above. 

Another  interesting  paper  concerning  open-path  TDLAS  is  [12]  where  the  authors 
demonstrate  a  system  that  is  very  similar  to  the  systems  described  above,  but  the 
application  is  for  the  detection  of  CO2  over  volcanically  active  areas  through  paths 
of  nearly  100  m,  although  the  authors  claim  that  with  the  retroreflector  used,  a  500 
m  path  is  possible  under  certain  atmospheric  conditions.  A  500  m  path  was  not 
demonstrated.  This  paper  uses  an  off-the-shelf  system  to  monitor  CO2  absorption 
lines  near  1.58  /mi  (6,329  crn^1). 

Another  paper,  [35]  describes  a  hand-held,  open-path  TDLAS  device  that  uses 
the  signal  reflected  from  background  objects  to  find  concentration.  This  device  was 


9 


limited  to  an  operating  distance  of  roughly  20  meters  depending  on  the  reflectivity  of 
the  background.  This  device  can  be  configured  to  detect  hydrogen  fluoride,  hydrogen 
sulfide,  or  methane.  Other  related  papers  with  the  same  authorship  include  [33] 
where  standoff  systems  for  trace  gas  detection  are  described,  and  a  similar  paper, 
[34]  demonstrates  a  combination  of  the  two  papers  above.  The  last  paper  gives  more 
details  about  a  finished  handheld  device  for  remote  gas  detection  including  ideas  for 
oxygen  sensors  at  762  nm.  Remote  vehicle  emission  monitoring  has  been  attributed  to 
the  individuals  of  the  next  interesting  paper.  [16]  This  group  started  remote  emission 
monitoring  in  the  late  1980s  with  the  first  remote-sensing  instrument  monitoring 
CO  and  CO2  using  broadband  infrared  spectroscopy  and  then  moving  on  to  NO  by 
including  UV  sources.  In  [16]  several  narrowband  open-path  experiments  concerning 
vehicle  emissions  are  summarized,  including  experiments  concerned  with  light  and 
heavy  duty  trucks,  locomotives,  snowmobiles,  and  airplanes.  These  papers  are  usually 
concerned  with  CO  and  CO2  found  at  their  well  known  fundamental  absorptions  at 
4.651  /jm  (2,150  cm-1)  and  4.255  /mi  (2,350  cm-1)  respectively.  The  articles  reviewed 
with  open-paths  consist  of  5-10  meter  paths,  usually  orthogonal  to  the  direction  of 
travel  so  that  the  constituents  from  the  vehicles  are  detected  as  the  vehicle  passes 
the  testing  site.  This  broadband  spectroscopy  suffers  with  a  path-length  limitation 
of  8  m  and  suffers  from  water  interference  around  the  hydrocarbon  channel,  but  [16] 
does  briefly  mention  an  interesting  problem  for  open-path  devices.  The  paper  refers 
to  the  failure  of  early  open-path  devices  due  to  a  “shimmering  effect”,  most  likely 
turbulence.  This  is  a  problem  often  ignored  by  most  remote  sensing  and  open-path 
TDLAS  papers.  In  summary,  of  the  papers  referenced  and  searched,  none  consider 
HEL  propagation,  very  few  papers  describe  problems  associated  with  turbulence  - 
though  some  have  symptoms  that  point  toward  turbulence  -  and  none  identify  or 
describe  solutions  for  turbulence  or  describe  how  to  quantify  turbulence  parameters. 
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None  of  the  open-path  TDLAS  devices  demonstrate  propagation  over  more  than  a 
few  hundred  meters. 

Simulation  Tools 

For  the  studies  outlined  in  this  document,  two  tools  were  used  extensively,  the 
Line-By-Line  Radiative  Transfer  Model  (LBLRTM)  and  the  High  Energy  Laser  End- 
to-End  Operational  Simulation  (HELEEOS).  AFIT-specihc  codes  were  used  to  esti¬ 
mate  alkali  absorption,  and  therefore  gain  profiles. 

LBLRTM 

The  Line-By-Line  Radiative  Transfer  Model  (LBLRTM)  was  originally  derived 
from  FASCODE.  LBLRTM  is  developed  and  maintained  by  Atmospheric  and  Envi¬ 
ronmental  Research,  Inc.  (AER).  [19]  Radiative  transfer  models  are  used  to  accu¬ 
rately  predict  atmospheric  effects  on  the  propagation  of  radiation  from  the  ultraviolet 
to  microwave  region  of  the  spectrum.  LBLRTM  examines  extinction  mechanisms  like 
molecular  scattering  and  absorption  and  aerosol  scattering  and  absorption  in  order  to 
characterize  the  ambient  atmosphere  and  determine  path  transmission.  LBLRTM  is 
so  heavily  based  in  FASCODE  that  its  inputs  still  use  the  legacy  “TAPEx”  formatting 
very  similar  to  the  input  formatting  originally  used  in  FASCODE,  a  legacy  format 
based  on  the  magnetic  tape  inputs  originally  used  for  input.  AER  validates  LBLRTM 
using  atmospheric  radiance  spectra  from  the  ultraviolet  to  the  sub-millimeter.  Inputs 
from  the  HITRAN  database  provide  line  parameters  used  by  LBLRTM  to  carry  out 
simulations  for  path  optical  depth  and  transmission.  LBLRTM  is  used  for  calcula¬ 
tion  of  infrared  atmospheric  radiance  and  its  derivatives,  used  to  compare  data  from 
sounding  profiles,  as  the  starting  point  for  many  retrieval  algorithms,  to  calculate 
absorption  coefficients  for  use  in  radiative  transfer  models  in  climate  applications,  to 
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train  systems  for  Numerical  Weather  Prediction  (NWP),  and  used  as  a  comparison 
to  satellite  data,  [19]  In  this  work  it  is  used  to  fit  to  spectral  data  collected  by  the 
TDLAS  apparatus. 

HELEEOS 

The  High  Energy  Laser  End-to-End  Operational  Simulation  (HELEEOS)  was  orig¬ 
inally  funded  by  High  Energy  Laser  Joint  Technology  Office  (HEL-JTO)  and  devel¬ 
oped  by  Air  Force  Institute  of  Technology  (AFIT).  [29]  Written  in  Matlab®  and 
distributed  as  an  executable,  the  simulation  models  nearly  every  aspect  of  a  HEL 
engagement  scenario.  Emerging  technologies  are  generally  supported,  and  a  variety 
of  different  engagements  can  be  simulated.  Air-to-surface,  surface-to-air,  air-to-air, 
air-to-space,  rclay-to-relay,  and  side  observations  are  supported  in  engagement  geome¬ 
tries.  Inputs  are  highly  configurable,  including  a  variety  of  laser  types,  geographic 
locations,  profiles,  aerosol  properties,  and  probabilistic  atmospheric  parameters. 
At  the  time  of  this  study,  HELEEOS  shares  many  calculation  roots  with  the  Scaling 
for  High  Energy  Laser  and  Relay  Engagement  (SHaRE)  toolbox,  currently  sponsored 
by  Air  Force  Research  Laboratory  (AFRL)  and  developed  by  MZA  Associates.  Using 
SHaRE  and  internal  calculations,  HELEEOS  considers  atmospheric  effects,  diffrac¬ 
tion  effects,  platform  disturbances,  thermal  blooming,  beam  control,  and  turbulence 
effects  on  beam  propagation  using  a  beam-spread  strchl  ratio  approach  resulting  in 
rapid  calculation  times.  The  Laser  Environmental  Effects  Definition  and  Reference 
(LEEDR)  [71]  is  used  in  HELEEOS  to  create  vertical  profiles  for  meteorological  data 
and  the  associated  environmental  effects  related  to  molecular  and  aerosol  absorption 
like  scattering,  optical  turbulence,  and  cloud  free  line  of  sight.  This  is  is  computed 
from  databases  and  extrapolated  for  use  in  any  global  location.  This  engagement 
level  simulation  has  inputs  related  to  engagement  geometry  and  atmospheric  pa- 
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rameters,  laser  source  and  platform  inputs,  target  lethality  and  platform  inputs  to 
allow  HELEEOS  to  calculate  output  metrics  related  to  the  slant  path  like  temper¬ 
ature,  pressure,  species  and  concentration,  wind  speed,  and  turbulence  to  compute 
final  engagement  metrics  like  molecular  absorption  and  scattering,  aerosol  absorption 
and  scattering,  thermal  blooming,  beam  strehl  ratios,  dwell  time,  and  probability  of 
kill.  Ultimately  HELEEOS  development  is  centered  around  the  goal  of  accurately 
predicting  total  energy  on  a  target  for  a  broad  range  of  engagement  scenarios.  [29] 
HELEEOS  and  its  internal  components  are  currently  maintained  and  upgraded  by 
the  AFIT  Center  for  Directed  Energy  (CDE). 

The  DPAL  systems  have  received  considerable  attention  and  systems  having  slope 
efficiencies  exceeding  80%  have  been  demonstrated  [96],  while  others  have  shown  ex- 
ciplex  assisted  pumping  schemes  [63] ,  and  while  most  DPAL  demonstrators  use  stable 
resonators,  some  have  used  transverse  excitation  with  unstable  resonators  [102],  The 
larger  spin-orbit  splitting  in  cesium  generally  requires  a  molecular  collision  partner  to 
accelerate  the  fine  structure  mixing,  [100]  but  rare  gases  are  also  usually  employed  at 
pressures  of  1-10  atmospheres  to  broaden  the  absorption  profile  and  better  match  the 
diode  laser  spectral  distribution  for  all  DPAL  variants.  [60,  61]  DPAL  systems  offer 
excellent  thermal  control,  good  beam  quality,  and  can  achieve  quantum  efficiencies 
of  95-99%.  A  high  power,  >1  kW,  system  with  48%  optical-to-optical  efficiency  has 
also  been  reported.  [14] 

Scope  of  Research 

This  work  sought  to  develop  an  open-path  TDLAS  system  that  operates  at  key 
HEL  wavelengths,  be  easily  configured  to  operate  at  different  wavelengths  in  the  field, 
and  can  collect  data  over  long  distances.  The  equipment  was  hardened  to  work  in 
outdoor  environments,  including  enclosing  the  laser  system  and  support  equipment, 
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installing  hardware  to  reduce  shock  and  vibration  during  deployment,  and  imple¬ 
menting  methods  for  preventing  damage  during  periods  of  inclement  weather  and 
transport.  The  complete  device  performance  was  characterized  using  molecular  oxy¬ 
gen,  a  simple  diatomic  molecule  that  is  both  abundant  and  relevant  to  DPAL  studies, 
in  an  indoor  and  an  outdoor  laser  test  range.  Additional  validation  was  achieved 
using  alternative  spectroscopic  devices  simultaneously,  like  a  Fourier  Transform  In¬ 
frared  (FTIR)  spectrometer.  Software  has  been  written  to  control  all  hardware  used 
during  experiments,  including  the  laser  source,  the  modulator,  the  lock-in  detectors, 
the  wavemeter,  the  calibration  laser,  the  data  acquisition  system,  and  any  additional 
hardware  that  may  be  required. 

Once  a  reliable  device  was  constructed,  collections  of  atmospheric  data  over  DPAL 
emission  wavelengths  were  attempted  over  open-paths.  With  successful  data  collec¬ 
tions,  the  atmospheric  constituents  most  important  to  DPAL  atmospheric  transmis¬ 
sion  was  investigated,  and  temperature,  pressure,  and  concentration  of  the  associated 
atmospheric  species  was  determined  using  spectroscopic  methods  and  post-processing 
routines,  all  written  from  scratch.  Results  from  open-path  TDLAS  collections  were 
compared  to  the  atmospheric  transmission  models  for  validation. 

Simulations  for  DPAL  gain  cells  were  written,  and  used  with  an  atmospheric 
transmission  model  to  estimate  DPAL  atmospheric  transmission  as  a  function  of 
DPAL  gain  cell  pressure.  Spectral  absorption  profiles  were  integrated  into  an  end- 
to-end  high  energy  laser  simulation  to  draw  conclusions  about  how  DPAL  gain  cell 
configurations  affect  a  variety  of  simulated  operational  scenarios. 

Modeling  of  spectral  gain  of  alkali  cells  in  DPAL  devices  has  previously  been 
investigated  for  cesium.  [5]  Using  the  cesium  code  as  a  model,  similar  codes  were  de¬ 
veloped  for  for  potassium  and  rubidium.  Each  alkali  transition  is  made  up  of  several 
hyperfme  lines,  and  each  line  is  broadened  by  thermal  and  collisional  effects.  The 
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gain  cell  profiles  from  these  do  not  incorporate  asymmetry,  which  has  been  found 
experimentally.  [52]  These  asymmetries  could  result  in  a  DPAL  laser  emission  profile 
that  is  wider  than  expected  based  on  the  simulations  shown  with  extremely  high  gain 
cell  pressures.  When  considering  spectral  DPAL  emission,  additional  complexity  is 
produced  by  the  resonator  used.  Given  that  the  resonator  will  likely  narrow  the  spec¬ 
tral  emission,  the  spectral  absorption  is  used  to  estimate,  and  perhaps  overestimate, 
the  width  of  DPAL  spectral  gain. 
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III.  A  Tunable  Diode  Laser  Absorption  System  for  Long 
Path  Atmospheric  Transmission  and  High  Energy  Laser 

Applications 


An  open-path  Tunable  Diode  Laser  Absorption  Spectroscopy  (TDLAS)  system 
composed  of  narrow  band  (300  kHz)  diodes  fiber  coupled  to  a  12”  Ritchey-Chretien 
transmit  telescope  has  been  developed  to  study  atmospheric  transmission  of  key  High 
Energy  Laser  wavelengths.  The  ruggedized  system  has  been  field  deployed  and  tested 
for  propagation  distances  of  greater  than  1  km.  Initial  experiments  were  performed 
in  the  vicinity  of  molecular  oxygen  to  electronic  transition  lines  near  760 

nm.  The  potassium  version  of  the  Diode  Pumped  Alkali  Laser  (DPAL)  operates  in 
between  two  of  the  sharp  oxygen  rotational  features  in  the  PP  and  the  PQ  branches. 
By  scanning  across  many  laser  free  spectral  ranges  and  monitoring  the  laser  frequency 
with  a  very  precise  wavemeter,  the  full  structure  of  the  oxygen  molecular  feature  is 
observed.  The  device  can  also  be  used  to  observe  rotational  temperatures,  oxygen 
concentrations,  and  total  atmospheric  pressure. 

Introduction 

Tunable  Diode  Laser  Absorption  Spectroscopy  (TDLAS)  is  an  active  field  of  re¬ 
search  with  an  international  conference  [2]  and  a  very  broad  range  of  research  areas 
including  performance  analysis  [51],  combustion  [15],  in-situ  analysis  [25],  materials 
manufacturing  [56],  and  much  more.  The  background  of  TDLAS  in  field  use  is  out¬ 
lined  well  by  Werle,  [86]  who  has  also  authored  a  review  article  on  post-processing  of 
TDLAS  data  [87] .  There  is  also  a  review  article  of  tunable  lasers  with  some  example 
uses  in  the  infrared  region  [93].  A  majority  of  TDLAS  articles  are  generally  focused 
on  using  absorption  cells,  such  as  White  or  Herriot  cells,  to  contain  the  sample  gas 
with  a  folded  optical  path  while  also  controlling  parameters  like  temperature  and 
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pressure.  A  small  subsection  of  TDLAS  research  is  concerned  with  open  paths  in 
place  of  absorption  cells.  For  example,  research  areas  include  topics  such  as  sensing 
greenhouse  gas  emissions  from  agricultural  sites  [68],  sensing  effluents  from  military 
cargo  aircraft  co-requisitely  with  other  instruments  [18],  performing  TDLAS  exper¬ 
iments  along  the  leading  edge  of  an  aircraft  [30],  over  5  and  15  m  open  paths  with 
retro  reflectors  [81],  and  sensing  NO2  using  a  telescope  over  a  160  m  path  [75].  There 
is  even  an  article  that  focuses  on  the  use  of  TDLAS  with  O2  [6].  In  this  chapter,  a 
diode  laser  is  used  to  perform  multi-mode  laser  absorption  spectroscopy  (MUMAS) 
in  the  spectral  region  of  molecular  O2  A3E ~  to  electronic  transition  lines  near 
760  nm  over  a  10  m  absorption  cell.  Temperatures  over  the  range  of  300  to  500  K 
were  found  with  a  precision  of  ±  2%  are  found,  pressures  varying  from  200  to  760 
torr  with  an  experiment  uncertainty  of  ±  2%,  and  concentration  measurements  that 
agree  with  estimated  pressure  using  the  ideal  gas  law  and  pressure  sensors  in  the  ab¬ 
sorption  cell.  While  the  operation  of  MUMAS  is  slightly  different  from  the  technique 
described  in  this  paper,  many  of  the  fundamental  principles  are  the  same.  The  device 
described  here  is  an  open-path  TDLAS  system  that  can  be  used  to  investigate  the 
overall  atmospheric  transmission  within  the  tunable  area  of  the  laser. 

TDLAS  Apparatus 

Hardware 

A  complete  system  diagram  illustrating  optical  and  electrical  paths  is  provided 
in  Fig.  1.  A  New  Focus  Velocity  model  6312  laser  diode  with  a  tunable  range  of 
approximately  759-770  nm  was  used  as  the  source.  This  laser  typically  achieves  an 
output  power  of  7  to  15  111W  and  has  a  spectral  width  of  less  than  300  kHz.  This 
spectral  width  is  small  compared  to  the  pressure-broadened  widths  of  atmospheric 
absorption  features.  An  Optics  for  Research  10-5  optical  diode  is  placed  in  the  beam 
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Figure  1.  Diagram  of  TDLAS  system  demonstrating  overall  layout. 

path  to  prevent  back  reflections  from  entering  the  laser  cavity. 

The  source  is  split  by  Thor  Labs  PB108/PB145B1  2  m  thick  pellicle  beamsplit¬ 
ters,  first  using  a  coated  beamsplitter  to  a  High  Finesse  WSU-2  wavemeter.  The  free 
spectral  range  that  results  from  a  very  thin  etalon,  A Vfsr  =  c/(2nd),  is  on  the  order 
of  tens  to  hundreds  of  GHz  for  typical  pellicle  thicknesses.  This  is  wider  than  the 
tunable  area  under  investigation,  and  is  easily  removed  with  the  reference  detector 
signal.  The  Fizeau-type  wavemeter  is  fiber-coupled  to  fiber  switch  that  is  in  turn 
coupled  to  the  laser  source  and  a  SIOS  SL-03  stabilized  HeNe  laser  used  for  cali¬ 
bration  to  achieve  ±  10  MHz  0.00002  nm)  accuracy  in  the  spectral  region  being 
investigated.  It  should  be  noted  that  most  TDLAS  experiments  use  Fabry-Perot  or 
similar  devices  for  spectral  calibration  in  post  processing,  and  can  offer  a  variety  of 
free  spectral  ranges.  Generally,  Fabry-Perot  etalons  require  calibration  themselves. 
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The  wavemeter  requires  little  post  processing  of  the  wavenumber  axis. 

The  remaining  power  of  the  source  laser  is  then  amplitude  modulated  by  a  Stan¬ 
ford  Research  Systems  SR540  mechanical  chopper  at  frequencies  around  2  kHz  and 
then  coupled  via  a  single  mode  fiber  to  a  RC  Optical  Systems  12  military  grade 
Ritchey-Chretien  telescope.  A  Ritchey- Chretien  telescope  is  similar  in  design  to  a 
Cassegrain  telescope,  but  has  hyperbolic  mirror  surfaces.  Before  expansion,  another 
pellicle  beamsplitter  is  used  with  a  Thor  Labs  PDA  100A  silicon  detector /preampli¬ 
fier  as  the  laser  power  reference.  The  beam  is  expanded  to  nearly  fill  the  primary 
mirror  of  the  telescope.  Fine  focusing  of  the  collimated  laser  can  be  achieved  with 
movement  of  the  secondary  mirror  of  the  telescope.  The  New  Focus  laser,  now  with 
5  mW  of  power,  is  expanded  to  fill  the  telescope  aperture  and  can  be  seen  plainly  by 
the  naked  eye  hundreds  of  meters  away.  It  should  be  noted  that  once  the  laser  is  ex¬ 
panded,  the  expanded  beam  falls  well  within  eye-safe  limits  without  viewing  optics. 
Collimation  is  verified  by  direct  viewing  at  the  end  of  the  path  with  a  calibration 
target  board  and  an  infrared  camera.  The  camera  is  moved  across  the  center  of  the 
expanded  beam  and  the  expanded  width  is  measured.  Adjustments  are  made  at  the 
transmitting  telescope  until  the  expanded  beam  at  target  is  nearly  the  same  width  at 
the  exit  aperture.  A  6  inch  mirror  is  placed  some  distance  away  from  the  telescopes 
and  used  to  reflect  a  portion  of  the  collimated  light  back  to  a  receive  telescope  of 
the  same  make  and  model.  A  Thor  Labs  PDA  100A  silicon  detector/preamplifier  is 
coupled  to  the  receive  telescope,  and  the  receive  telescope  is  carefully  focused  so  that 
all  received  collimated  light  falls  on  the  detector.  To  aid  in  alignment,  an  eye  safe 
HeNe  laser  expanded  to  approximately  3  cm  is  attached  to  the  transmission  telescope 
is  co-aligned  with  the  expanded  laser  beam  and  used  as  a  reference  for  the  mirror 
alignment  process  to  redirect  the  tunable  laser  back  to  the  receive  telescope.  With 
experienced  personnel  and  hardware  in  place,  collimation  and  alignment  of  the  system 
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takes  about  10  minutes  over  typical  paths. 

The  detector  signals  are  analyzed  by  Stanford  Research  Systems  SR-850  dual- 
phase  lock-in  amplifiers.  The  reference  output  of  the  chopper  is  used  by  the  lock-in 
to  automatically  detect  the  chopper  frequency  and  the  phase  of  the  chopped  signal. 
The  time  constant  of  the  lock-in  is  adjusted  to  be  much  shorter  than  the  spectral 
scan  speed  of  the  laser,  and  the  phase  and  gain  of  each  lock-in  are  adjusted  to  give 
the  best  signal  over  the  dynamic  range  expected.  For  the  experiments  shown,  lock-in 
time  constants  were  typically  100  ms  with  a  filter  drop  off  of  about  6  dB/octave  and 
a  scan  speed  of  about  0.003  nm/sec.  An  iris  located  at  the  entrance  aperture  of  the 
reference  detector  is  adjusted  so  that  the  relative  signals  reported  by  the  lock-ins  for 
the  absorption  signal,  /,  and  the  reference  signal,  Jo,  are  of  similar  magnitude.  The 
dominant  noise  sources  from  (largest  to  smallest)  during  a  collection  are  wind  buffet¬ 
ing  of  telescopes  and  mirror,  turbulence  effects  over  the  path  being  investigated,  error 
in  tunable  laser  frequency  measurement,  and  response  variations  between  detectors. 
The  lock-in  magnitude  and  phase  along  with  laser  control  system  output,  and  the 
chopper  response  are  recorded  by  a  computer  using  a  National  Instruments  USB- 
6251  data  acquisition  device.  This  DAQ  is  also  capable  of  synchronous  analog  and 
digital  outputs,  and  can  be  used  to  drive  other  systems  in  the  experiment  from  laser 
warning  lights  to  laser  wavelength.  All  hardware  is  secured  in  a  frame  with  vibration 
shock  mounts.  The  optical  board  is  secured  on  the  top  for  easy  adjustments  to  optics 
in  the  held,  as  shown  in  Fig.  2. 

Software 

Instruments  are  controlled  via  Matlab®  with  the  Instrument  Control  Toolbox 
and  a  GPIB  interface  with  an  Agilent  82357A  USB-to-GPIB  converter.  The  National 
Instruments  USB-6251  DAQ  is  controlled  via  Matlab®  with  the  Data  Acquisition 
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Figure  2.  Photo  of  TDLAS  system  with  optical  board. 


Toolbox  through  a  USB  interface.  The  NI  USB-6251  DAQ  also  can  produce  analog 
and  digital  control  outputs,  and  drives  the  piezo  in  the  New  Focus  laser  cavity  for 
fine  wavelength  tuning.  The  laser  power  is  monitored  with  a  reference  detector  (Jo) 
just  before  the  open-path  to  be  investigated.  The  precise  wavemeter  used  in  this 
system  is  accessed  through  a  custom  Matlab®  interface  using  dynamically  linked 
libraries  provided  by  the  manufacturer.  Tuning  of  the  laser  is  achieved  with  a  rough 
motor  adjustment  via  the  GPIB  interface  and  has  a  repeatable  accuracy  of  about  0.1 
nm.  Fine  tuning  with  a  piezo  flexing  mechanism  scans  over  a  0.2  nm  window  via 
the  digital  to  analog  interface  on  the  DAQ.  The  full  spectrum  is  scanned  by  coarsely 
tuning  to  a  point  in  the  spectrum,  finely  tuning  through  a  0.2  nm  window,  then 
coarsely  moving  to  the  next  region  in  the  spectrum  until  the  full  tunable  range  of  the 
diode  has  been  explored.  The  wavemeter  is  calibrated  with  a  wavelength  stabilized 
SIOS  SL-03  HeNe  laser  where  calibration  takes  place  before  each  scan.  The  output 
of  the  wavemeter  is  collected  in  sync  with  the  DAQ  sampling  of  the  lock-in  outputs. 
The  wavemeter  checks  the  course  tuning  of  the  laser  so  that  spectral  areas  will  not  be 
missed  and  optimizes  scans  to  prevent  too  much  spectral  overlap.  After  each  piezo 
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scan,  the  Matlab®  interface  backs  up  data  to  the  hard  disk,  and  stores  individual 
device  statistics,  all  in  an  effort  to  diagnose  potential  problems  discovered  from  post 
processing  or  in  the  event  of  power  loss  during  an  experiment.  Uninterruptible  power 
supplies  are  used  prevent  damage  to  equipment  from  unreliable  power  sources  and 
prevent  loss  of  data  during  collection. 

Joining  FSR  scans 

Because  a  full  data  collection  is  achieved  using  many  overlapping  finely  tuned  piezo 
scans,  these  scans  must  be  joined  together  to  form  a  full  spectrum.  Initially,  the  spec¬ 
trum  itself  formed  a  complete  spectrum  using  a  crosscorrelation  between  adjoining 
scans  in  conjunction  with  known  line  positions  of  O2  absorption  features.  Spectral 
absorption  features  gave  a  clear  indication  of  how  a  set  of  piezo  scans  could  fit  to¬ 
gether.  However,  nonlinearity  in  each  piezo  scan  combined  with  spectral  areas  with 
no  absorption  features  made  the  merging  of  piezo  scans  a  difficult  task  to  carry  out 
accurately.  Even  when  piezo  scans  were  merged  together,  each  piezo  scans  nonlinear 
spacing  in  wavenumber  resulted  in  poor  information  about  the  lincshape  of  spectral 
features  and  the  integrated  area  for  each  spectral  feature.  A  piezo  scan  that  was  as¬ 
sumed  to  be  linear  in  wavenumber  could  be  in  error  by  as  much  as  half  a  wavenumber. 
A  precise  and  quickly  sampled  wavemeter  was  used  as  the  reference  throughout  each 
piezo  scan.  Typically  the  wavemeter  sampled  wavelength  on  the  order  of  hundreds  of 
Hz  for  a  scan  over  many  seconds.  For  each  piezo  scan,  wavemeter  output  was  used 
to  find  the  spectral  area  where  two  scans  overlapped  with  reliable  laser  output.  An 
example  of  the  overlapping  spectral  area  of  two  adjacent  piezo  scans  is  shown  in  Fig. 
3. 

Here,  the  reliable  overlap  between  spectral  scans  is  nearly  0.5  cm-1.  The  overall 
nonlincarly  from  sample  to  sample  can  also  be  seen  in  the  gentle  curvature  of  each 
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Figure  3.  Overlapping  spectral  region  for  two  adjoining  FSR  piezo  scans. 

scan,  and  the  similar  nonlinear  structure  is  observed  between  different  piezo  scans. 
Once  the  reliable  spectral  overlap  is  determined,  the  best  point  that  scans  can  be 
joined  is  determined  by  finding  the  overlapping  point  where  the  difference  between 
raw  signal  data  from  each  scan  is  the  smallest.  This  reduces  the  chance  for  anomalous 
jumps  in  data  and  preserves  the  lineshape  of  spectral  features  caught  between  scans. 
The  overlapping  area  of  two  scans  can  be  seen  in  Fig.  4  with  boundaries  where  the 
overlapping  area  is  considered  unreliable,  and  where  the  particular  switch  off  point 
between  adjacent  scans  lies.  The  residuals  shown  are  the  difference  between  the  final 
signals  of  each  scan  over  the  reliable  overlapping  area. 

This  process  is  performed  between  each  piezo  scan  to  form  the  full  spectrum.  The 
resulting  data  has  a  nonlinear  spacing  between  data  points  and  is  interpolated  to 
form  a  linearly  spaced  spectral  axis  with  the  number  of  points  required  by  the  final 
application.  Typically  a  full  scan  over  10  nm  is  generally  over  one  million  data  points 
after  piezo  scans  are  stitched  together.  Once  completed,  the  data  set  contains  stitched 
data  for  wavenumber,  voltage  applied  to  the  piezo,  the  rough  output  wavelength  from 
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Figure  4.  Demonstration  of  the  choice  for  spectral  area  overlap  of  ad¬ 
joining  piezo  scans. 


the  laser  as  a  function  of  voltage,  x  and  y  outputs  for  each  lock-in  amplifier,  and  a 
synchronized  absolute  time-of-sample  channel.  In  order  to  simulate  the  spectrum, 
have  expectations  about  transmission,  or  fit  to  collected  data,  parameters  about  the 
molecule (s)  of  interest  must  be  known.  Generally,  these  parameters  are  contained 
in  the  HITR.AN  database  in  the  forms  of  molecule,  isotope,  abundance,  Einstein 
coefficient,  broadening  rate,  temperature  shift,  and  the  like.  In  order  to  access  these 
parameters  in  MATLAB®  an  interface  was  created  to  perform  searches  for  atmospheric 
absorption  lines  over  particular  regions  of  the  spectrum.  This  interface  is  similar  to 
JavaHAWKS,  but  has  its  interface  in  Matlab®. 

Baseline  considerations 

The  usable  data  from  an  experiment  consists  of  the  signal  from  the  receive  tele¬ 
scope  detector  (/)  and  the  signal  from  the  laser  power  reference  detector  (/0)  both 
referenced  to  the  same  linearly  spaced  wavenumber  axis.  This  data  has  some  baseline, 
a  likely  result  of  nonlinearity  in  detectors  resulting  from  temperature  changes  over 
time  of  the  detection  and  inherent  in  their  manufactured  nature  (the  preamplifiers  of 
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Figure  5.  I,  Iq,  and  ln(///o)  with  an  instrument  baseline  visible  over  a 
102  m  path. 


PDAlOOAs  are  especially  sensitive  to  temperature).  This  baseline  must  be  factored 
out  before  fitting  to  the  data,  or  must  be  considered  during  fits.  An  example  of  the 
reference  signal  from  each  detector  and  the  result  of  the  using  /  and  1$  in  concert 
is  shown  in  Fig.  5. 

Typically,  the  baseline  is  mostly  linear,  with  some  minor  curvature  if  the  environ¬ 
ment  where  testing  takes  place  stays  at  a  relatively  constant  temperature.  For  the 
data  shown  above,  fits  to  the  baseline  with  a  simple  quadratic  generally  produce  a 
squared  term  4  orders  of  magnitude  smaller  than  the  slope  term.  For  a  higher  order 
fit,  a  cubic  term  is  8  orders  of  magnitude  smaller  than  the  slope  term  with  little 
change  in  the  resulting  standard  error  of  the  fit,  thus  demonstrating  that  there  is 
little  dependence  on  a  cubic  term  for  a  baseline  fit.  Etelon  effects  from  optics  before 
the  reference  detector  can  also  be  seen,  predominantly  in  green  curve  (J0)  but  also 
visible  in  the  blue  (/),  in  Fig.  5,  but  are  then  removed  when  the  two  signals  are  used 
together.  One  critical  point  of  discussion  is  that  removing  this  baseline  also  removes 
effects  from  broadband  aerosol  scattering,  an  effect  that  can  be  significant  depend- 
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Figure  6.  The  ln(///0)  compared  to  a  multiple  profiles  and  fit  residuals. 

ing  on  the  atmospheric  aerosols  present.  The  capability  to  retrieve  this  information 
cannot  currently  be  obtained. 

Spectral  Lineshapes 

Each  isolated  spectral  line  is  fit  with  an  individual  Voigt  lineshape,  defining  the 
line  center,  line  amplitude  (or  spectral  area  of  the  line),  the  Doppler  (Gaussian) 
temperature-broadened  width,  and  the  pressure-broadened  (Lorentzian)  width.  All 
four  parameters  are  allowed  to  vary  during  the  fitting  process.  The  least-squares 
fitting  routine  in  PeakFit®  was  iterated  until  the  fitting  process  converged  to  a 
solution  giving  the  final  result  of  the  fit  to  the  data  is  illustrated  in  Fig.  6  where  each 
individual  line  is  shown  to  be  well  isolated,  and  signal  to  noise  is  on  the  order  of  200 
for  the  larger  absorption  features. 

The  residuals  have  structure  for  each  absorption  feature  although  the  residuals  are 
very  small  and  are  typically  less  than  1%  of  the  absorption  feature.  The  least-squares 
fit  to  all  lines  shown  resulted  in  an  r2  of  0.999  and  a  fit  standard  error  of  0.0028.  The 
structured  residual  shown  in  Fig.  7  is  a  close-up  view  of  a  single  absorption  feature 
from  Fig.  6. 
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Sample  Point 


Figure  7.  Example  of  structure  in  residuals. 

The  direct  cause  of  this  residual  is  unknown  but  has  been  observed  in  other  lit¬ 
erature  [67].  Ritter  experienced  strikingly  similar  residuals  in  Voigt  fits  to  oxygen 
A-band  data.  The  W-shaped  residuals  Ritter  demonstrates  were  experienced  at  sev¬ 
eral  different  pressures  ranging  from  76  torr  to  600  torr  and  carry  the  same  structure 
throughout  pressure  variations.  It  should  be  noted  that  these  residuals  have  minimal 
effect  on  the  quality  of  data  obtained  and  speak  more  to  the  methods  used  for  fitting 
and  post  processing.  The  residual  depicted  in  Fig.  7  were  a  result  of  not  finding  the 
global  minimum  in  PeakFit®.  After  extensive  work,  the  structure  in  these  residu¬ 
als  were  eliminated  using  both  PeakFit®  and  non-linear  least  squares  methods  in 
MATLAB®,  the  results  of  which  are  demonstrated  in  Chap.  IV. 

Conclusions 

The  TDLAS  system  is  described  in  terms  of  hardware  and  software,  and  shown  to 
collect  information  for  molecular  absorption  over  atmospheric  paths  from  hundreds  to 
potentially  thousands  of  meters.  Full  spectra  can  be  collected  for  the  tunable  range 
of  the  laser  source  used,  and  the  laser  source  itself  can  easily  be  changed  to  produce 
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spectra  over  many  different  tunable  ranges.  Accuracy  in  wavenumber  is  limited  to 
the  spectral  reference  used  and  is  approximately  ±  10  MHz.  Noise  is  primarily  domi¬ 
nated  from  effects  of  telescope  and  mirror  buffeting.  Spectra  from  a  102  m  path  have 
been  presented  with  signal  to  noise  of  over  200.  The  high-quality  spectral  scans  from 
the  system  can  be  used  for  a  variety  of  tasks,  including  investigation  of  multi-species 
molecular  absorption,  estimation  of  temperature,  pressure,  and  concentration  of  at¬ 
mospheric  species  with  spectroscopic  methods,  high  energy  laser  (HEL)  atmospheric 
transmission,  and  validation  of  Monocular  Passive  Ranging  (MPR)  techniques. 


IV.  Investigation  of  atmospheric  O2  X3Efy  to  using 
open-path  tunable  diode  laser  absorption  spectroscopy 


Applied  A  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  device  fiber  cou¬ 
pled  to  a  pair  of  12.5”  telescopes  was  used  to  study  atmospheric  propagation  for  open- 
path  lengths  of  100  to  1,000  meters.  More  than  50  rotational  lines  in  the  molecular 
oxygen  A-band  AT3E“  to  b transition  near  760  nm  were  observed.  Temperatures 
were  determined  from  the  Boltzmann  rotational  distribution  to  within  1.3%  (less 
than  2  K).  Oxygen  concentration  was  obtained  from  the  integrated  spectral  area  of 
the  absorption  features  to  within  1.6%  (less  than  ±  0.04  x  1018  molecules  /  cm3). 
Pressure  was  determined  independently  from  the  pressure-broadened  Voigt  lineshapes 
to  within  10%.  A  Fourier  Transform  Interferometer  (FTIR)  was  also  used  to  observe 
the  absorption  spectra  at  1  cm-1  resolution.  The  TDLAS  approach  achieves  a  mini¬ 
mum  observable  absorbance  of  0.2%,  whereas  the  FTIR  instrument  is  almost  20  times 
less  sensitive.  Applications  include  atmospheric  characterization  for  high  energy  laser 
propagation  and  validation  of  monocular  passive  ranging. 

Introduction 

Atmospheric  transmission  in  the  vicinity  of  the  O2  X3S~  to  blYi+  absorption  fea¬ 
ture  is  important  for  remote  sensing  [25]  and  military  applications  [15].  Environmen¬ 
tal  applications  often  require  the  simultaneous  monitoring  of  multiple  atmospheric 
and  pollutant  concentration,  which  has  stimulated  the  development  of  multiplexing 
and  frequency  modulation  techniques.  [51]  In  contrast,  atmospheric  propagation  for 
high  energy  lasers  requires  high  spectral  resolution  and  low  detection  limits  over  long 
paths.  The  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  technique  is  well 
developed  and  desirable  for  such  applications.  [86] 

The  majority  of  TDLAS  experiments  are  performed  in  multi-path,  White,  or 
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Herriot  cells  where  the  temperature  and  pressure  are  readily  controlled.  A  few  appli¬ 
cations  for  open-path  instruments  include  monitoring  greenhouse  gas  concentrations 
from  agricultural  sites  with  10%  uncertainties  [69],  sensing  effluents  from  military 
cargo  aircraft  with  path  lengths  of  up  to  15  m  [18],  determining  atmospheric  con¬ 
stituents  onboard  commercial  aircraft  [30],  investigating  methane  and  ethane  over  5 
and  15  m  open  paths  [81],  and  sensing  NO2  with  telescopic  instruments  over  path 
lengths  exceeding  160  m  [75].  The  challenges  of  developing  TDLAS  instruments  for 
paths  of  1  km  or  greater  where  jitter  control,  turbulence  effects,  and  held  operations 
become  important,  are  largely  unaddressed. 

Recently  Arita,  et  al.  [7]  examined  the  O2  A"3E~  to  blYi+  transition  using  multi- 
mode  laser  absorption  spectroscopy  (MUMAS)  for  a  10  m  absorption  cell.  Spectra 
were  recorded  at  temperatures  of  300-500  K,  and  pressures  of  200-760  Torr.  Fitting 
simulated  spectra  to  the  observations  yields  uncertainties  in  extracted  temperatures 
of  8  K  and  pressures  of  less  than  14  Torr.  Concentration  was  reported  with  2% 
confidence  limits.  [7]  The  current  TDLAS  results  will  be  compared  to  this  previous 
MUMAS  work. 

The  oxygen  CUA"3!]"  to  51S+  transition  is  of  particular  interest  for  passive  ranging. 
Techniques  estimating  range  using  the  depth  of  spectral  absorption  features  when 
viewing  a  distant  bright  broadband  spectral  source  is  often  referred  to  as  monocular 
passive  ranging  (MPR).  [53]  The  O2  (X-b)  transition  is  desirable  for  passive  ranging 
for  several  reasons:  1)  it  has  a  reasonably  constant  atmospheric  concentration  ratio  for 
dry  air,  2)  is  characterized  well  as  a  function  of  altitude  and  meteorological  conditions, 
and  3)  is  spectrally  isolated  from  other  absorbing  species. 

A  demonstration  of  the  MPR  technique  for  a  static  rocket  motor  test  using  a 
Bomcm  MR-254  Fourier  transform  spectrometer  at  a  range  of  2.8  km  yielded  range 
estimates  with  an  accuracy  of  about  0.5%  (14  m).  [39]  More  recently,  an  intensified 
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CCD  array  coupled  to  variable  band  liquid  crystal  filter  was  deployed  for  a  ground 
test  of  a  static  jet  engine  in  afterburner  at  ranges  of  0.35  -  4.8  km,  establishing  a 
range  error  of  15%.  [4]  The  current  work  validates  the  MPR  FTIR  approach  with  the 
TDLAS  active  instrument.  A  second  application  for  the  TDLAS  instrument  involves 
characterizing  molecular  and  aerosol,  absorption,  and  scattering  for  high  energy  lasers. 
[57]  In  particular,  the  potassium  variant  of  the  Diode  Pumped  Alkali  Laser  (DPAL) 
operates  near  770.1  nm,  in  the  high  rotational  limit  of  the  O2  (X-b)  PQ  and  pP 
branches.  [105]  Studying  the  atmospheric  effects  of  high  irradiance  (kW/cm2)  laser 
propagation  to  distant  targets  will  be  greatly  aided  by  the  development  of  this  rugged 
TDLAS  instrument. 

The  TDLAS  Device 

The  tunable  diode  laser  absorption  spectroscopy  instrument  is  shown  schemati¬ 
cally  in  Fig.  8. 

A  detailed  description  of  the  device  configuration  and  method  for  collecting  spectra 
have  previously  been  discussed.  [66]  The  tunable  diode  is  a  New  Focus  Velocity  laser 
model  6312  with  a  100  mW  maximum  output,  a  10  nm  tunable  range,  and  less 
than  a  300  kHz  linewidth.  The  laser  source  near  760  nm  can  easily  be  changed 
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to  investigate  other  spectral  regions  of  interest  with  minimal  changes  to  the  overall 
system.  After  amplitude  modulation  at  a  frequency  of  2  kHz,  the  laser  is  fiber  coupled 
and  expanded  to  nearly  fill  a  military  grade  12.5  inch  RC  Optical  RitcheyChretien 
transmit  telescope.  The  laser  beam  is  directed  across  an  open-path  and  then  received 
by  a  second  identical  telescope.  Two  Thor  Labs  PDAlOOAs  are  used  for  the  reference 
and  signal  detectors  and  are  analyzed  by  Stanford  Research  Systems  SR850  dual¬ 
phase  lock-ins  and  recorded  by  a  National  Instruments  USB-6251.  Pellicle  beam 
splitters  are  used  to  minimize  etalon  effects.  The  reference  intensity  is  measured 
late  in  the  optical  train  (immediately  prior  to  the  telescope  beam  expansion)  and 
attenuated  to  balance  sent  and  received  signals.  A  High  Finesse  WSU-2  wavemeter 
is  used  to  determine  the  frequency  axis  and  is  calibrated  with  a  fiber-coupled  SIOS 
SL-03  frequency-stabilized  HeNe  laser  to  achieve  at  better  than  10  MHz  accuracy. 
All  hardware  is  connected  to  and  controlled  by  a  PC  using  MATLAB®.  During  an 
experiment,  the  laser  is  finely  tuned  by  driving  a  piezo  to  tune  over  an  approximately 
0.1  nm  range,  and  then  coarsely  tuned  to  the  next  spectral  region.  The  process  is 
repeated  for  95  free  spectral  ranges  to  cover  a  10  nm  spectral  region.  Calibration  of 
the  frequency  axis,  baseline  removal,  and  further  spectral  processing  have  recently 
been  reported.  [66] 

System  noise  is  generally  dominated  by  telescope  jitter  and  atmospheric  turbu¬ 
lence  over  the  path,  although  for  paths  of  100  meters,  these  effects  are  usually  small. 
Spectra  have  been  recorded  over  100  m  to  1  km  open-paths,  and  the  maximum  open- 
path  distance  is  expected  to  be  approximately  5  km,  a  limitation  of  the  tip  and 
tilt  pointing  accuracy  of  the  QuickSet  QPT-130  ruggedized  pan  and  tilt  mounts  used. 
Collection  geometries  consist  of  the  transmit  and  receive  telescopes  next  to  each  other 
separated  by  1  meter,  each  viewing  a  6  inch  mirror  placed  at  a  variety  of  distances 
giving  total  path  lengths  from  100-1,000  meters.  1  km  collections  took  place  outdoors 
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after  several  inches  of  snowfall  with  an  average  temperature  of  -8  C,  a  48%  relative 
humidity,  and  a  pressure  of  989  hPa  during  the  collection  time. 

The  performance  of  the  active  TDLAS  system  was  compared  to  a  passive  Fourier 
Transform  Spectrometer.  An  ABB  Bomern  MR-304  FTIR  spectrometer  coupled  to  a 
reflective  telescope  with  a  4.9  mrad  held  of  view  and  12  inch  aperture  was  placed  in 
between  the  TDLAS  transmit  and  receive  telescopes  and  viewed  a  700  Watt  halogen 
broadband  source  placed  above  the  FTIR  viewed  through  the  same  mirror  described 
above.  The  path  viewed  by  the  TDLAS  and  FTIR  are  nearly  the  same.  This  FTIR 
uses  a  silicon  (Si)  detector  with  a  spectral  range  of  8,000-15,000  cm-1  and  an  indium 
gallium  arsenide  (InGaAs)  detector  with  a  spectral  range  of  6,000-12,000  cm-1  and 
spectral  resolution  of  1  cm-1.  The  Si  detector  provided  superior  signal-to-noise  in  the 
760  nm  spectral  region.  The  FTIR  device  collected  interferograms  at  10  Hz  resulting 
in  spectra  with  1  cm-1  resolution.  Calibration  and  post  processing  techniques  have 
previously  been  described.  [39] 

Two  NIST-certified  Davis  Vantage  Pro2  weather  sensors  were  placed  near  the 
telescope  and  mirror.  Temperature,  pressure,  and  relative  humidity  measurements 
were  recorded  with  instrument  variations  of  1  C,  2  hPa,  and  6%  respectively. 

Results 

A  typical  TDLAS  spectrum  in  the  vicinity  of  the  molecular  oxygen  02V3S"  to 
(v’=0,  v”=0)  band  is  shown  in  Fig.  9a  with  approximately  95  free  spectral 
ranges  and  over  a  million  data  points  in  the  full  collected  spectrum  with  only  every 
25th  data  point  plotted  in  the  figure. 

The  rotational  spectrum  of  the  (%  (X-b)  absorption  is  assigned  as  AK  A  J(K") 
where  K  is  the  total  angular  momentum  without  electron  spin,  and  J  represents  the 
total  angular  momentum  (spin  +  rotation).  Fig.  10  shows  the  energy  levels  involved 
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Figure  9.  (a)  Undersampled  TDLAS  spectrum  of  O2  (X-b)  (0,0)  band 
with  rotational  assignments,  (b)  fit  of  Equation  (9)  using  a  series  of  Voigt 
lineshapes,  (c)  differences  between  observed  spectra  and  simulation  using 
weather  instrument  data,  and  (d)  difference  between  simulation  and  data 
using  TDLAS  derived  atmospheric  parameters. 
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Figure  10.  Example  energy  level  diagram  of  O2  (X-b)  for  the  pP( 3)  and 
PQ{ 3)  lines. 


for  the  pP(3)  and  PQ(3)  transitions. 

The  spectrum  of  Fig.  9  exhibits  two  of  the  four  rotational  branches, 
aaA J  —p  PpQrRrQ,  described  by  Hnnd’s  case  (b)  coupling  for  the  magnetic 
dipole  transition.  [8]  The  R-branch  lines  up  to  K"  =  27  are  observable  with  the  cur¬ 
rent  apparatus.  However,  the  current  analysis  is  limited  to  the  P-branch  where  the 
rotational  spacing  is  larger.  Only  the  odd  values  of  K  are  present  due  to  nuclear  spin 
statistics.  The  rotational  distribution  in  the  ground  state  peaks  at  K"  =  9,  consistent 
with  a  temperature  near  300  K.  The  long  path  length  of  100  m  and  lack  of  instru¬ 
mental  broadening  leads  to  a  large  peak  absorbance  of  A  =  2.36  ±  0.002,  despite  the 
long  radiative  lifetime  of  ~  11.4  seconds.  [73]  The  spectral  resolution  is  limited  by 
the  pressure  broadened  lineshapes  and  complete  spectral  isolation  is  achieved. 

The  absorbance,  A,  for  each  rotational  feature  is  described  by  the  Beer-Lambert 
law: 

Aak,aj,k"(v)  =  —  hr  (  =  (J akaj,k"(v)N(K" ,  J")L  (9) 

Vo  (v)J 

where  the  natural  logarithm  of  the  send,  Jo,  and  received,  It ,  signals  is  equal  to  the 
product  of  the  frequency  dependent  absorption  cross-section  of  the  specific  rational 
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feature,  cak.aj.k",  the  number  density  of  molecular  oxygen  in  the  specific  rotational 
level,  N(K",  J")  and  the  optical  path  length,  L.  I o  and  It  are  initially  balanced  using 
an  iris  on  the  Jo  detector  to  allow  similar  gain  settings  between  both  detectors  and 
their  respective  lock-ins.  The  absorption  cross-sections  are  specified  as, 

,  s  gj>  Aak,aj,k"Xak,aj,k"9v(v*^  v0(AK,  AJ,  K"),  Avl,  Avd) 
vak,aj,k"{v)  =  - - - — -  (10) 


g.j" 


8n 


where  the  degeneracies  for  the  ground  and  exited  states  are,  gj"  =  (2J"  +  1)  and  gj>  = 
(2 J'  +  1),  and  Av^  and  Av o  are  the  Lorentzian  and  Doppler  linewidths  (FWHM), 


AvL  =  (^)'!  (7 airip  —  Ps)  +  IselfPs) 
AvD  «  7.2xlO_7u0(A/i,  A  J,  K") 


(11) 


where  7 air  and  7 seif  are  the  air  and  self-broadening  coefficients,  p  and  ps  are  the  total 
pressure  and  partial  pressure,  and  n  is  the  coefficient  of  temperature  dependence 
on  the  air-broadened  linewidth,  all  defined  by  the  2008  HITRAN  database.  [73] 
M  is  the  species  mass  in  AMU;  T  is  the  temperature  of  the  atmosphere  in  Kelvin. 
The  state  specific  spontaneous  emission  coefficients,  Aaa'.aj./c",  and  line  positions 
uq(AK,  AJ,  K”)  are  also  referenced  from  HITRAN.  The  Voigt  lineshape  gy  is  area 
normalized  so  that  the  frequency  integrated  absorption  is, 


aAK,AJ,K"  —  /  aAK,AJ,K"(l')d1'  ~ 


gj,  A  ax,  A  J,K'AaK,AJ,K" 


gjn 


87T 


(12) 


The  Boltzmann  distribution  specifies  the  rotational  dependent  number  density, 


NIK"  J")  = 

1  ’  '  Q(T) 


(13) 


where  N  is  the  total  concentration  of  oxygen,  Q(T )  is  the  HITRAN  rotational  parti- 


36 


Figure  11.  A  fit  of  Equation  (9)  to  the  pP( 5)  rotational  line. 

tion  function  that  is  approximately  equal  to  hcB/kT  using  rotational  constant  B ,  / 
is  the  isotopic  abundance  of  1602,  and  is  the  rotational  energy  for  the  ground 

state. 

Fig.  9a  shows  a  fit  of  Eq.  (9)  to  each  line  in  the  observed  spectrum.  All  lines  in  the 
spectra  were  fit  simultaneously  with  Voigt  lincshapes  using  a  nonlinear  least  squares 
method.  Each  absorption  feature  was  described  using  four  parameters:  line  center, 
uo,  integrated  absorbance,  A0  =  f  A(is)dis,  Doppler  width,  Az/d,  and  the  Lorentzian 
width,  Avl-  Additionally,  a  small  baseline  (A  «  0.0019)  was  observed  and  a  cubic 
background  was  included.  The  RMS  fit  residual,  0.0022,  is  less  than  0.1%  of  the  peak 
absorbance,  as  shown  in  Fig.  9b.  This  detection  limit  is  comparable  to  that  achieved 
for  controlled  environment  multi-path  cells  without  multiplexed  detection.  The  fit 
residuals  are  unstructured  and  the  spectra  exhibit  a  signal-to-noise  ratio  of  about 
150.  Fig.  11  demonstrates  a  Voigt  fit  to  a  single  line,  the  pP( 5),  in  the  non-physical 
log-absorbance  to  demonstrate  the  quality  of  signal. 

Fig.  9c  shows  the  difference  between  the  data  and  a  simulation  using  the  Line- 
by-line  Radiative  Transfer  Model  (LBLRTM)  with  average  atmospheric  conditions 
reported  by  the  meteorological  instruments  ( T  —  25.1  ±  1  C,  P  —  0.977  ±  0.002 
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atm,  N  =  4.911  ±  0.16  x  1018  cm  3)  while  Fig.  9d  is  an  LBLRTM  simulation  using 
the  fitted  results  (T  =  30.2  C,  P  =  1.035  ±  0.053  atm,  N  =  4.825  ±  0.04  x  1018 
cm-3).  Modest  improvements  resulted  from  using  the  temperature,  pressure,  and 
concentration  derived  from  the  spectra  rather  than  the  meteorological  instruments, 
but  very  small  errors  in  wavelength  prevent  smaller  residuals  unless  line  centers  were 
included  as  a  fit  parameters.  The  fit  line  positions  differ  from  the  HITRAN  data 
base  by  0.00077  ±  0.00093  cm-1,  which  compares  favorably  with  the  wavemeter  best 
accuracy  of  0.00033  cm-1. 


Rotational  Temperature 

Temperature  is  readily  determined  from  the  rotational  intensity  distribution.  The 
integrated  absorbance  for  each  feature, 
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gj<  Aaa',aj,/<"^aa'aja'"  (2 J"  +  1)  ~ ek",j » 
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Qj"  8vr  Q(T) 


NLf 


(14) 


can  be  evaluated  to  achieve  a  linear  dependence  on  rotational  energy, 


4° 

aAK,AJ,K" 


(2J"  +  l)-^Aft',A  J,K"NaK,AJ,K'‘ 


NLf 
87 tQ(T) 


Independent  weighted  linear  fits  of  the  observed  absorbance  for  both  the  p P  and  PQ 
branches  are  are  illustrated  in  Fig.  12. 

The  variance  in  fit  absorbance,  A0  ,  as  defined  by  the  95%  confidence  bound  for 
the  Voigt  fit  of  each  absorption  feature,  were  employed  for  the  weights.  The  pP 
and  the  PQ  branches  resulted  in  temperatures  of  303.5  ±  2.0  K  and  303.1  ±  1.8 
K,  respectively.  The  uncertainty  in  temperature  is  defined  by  the  95%  confidence 
bounds  in  the  slope  parameter.  The  meteorological  instruments  used  during  the 
collection  recorded  temperature  every  minute,  and  the  average  temperature  over  the 
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Figure  12.  Estimations  of  temperature  using  the  intensity  distributions 
of  the  pP  and  PQ  branches. 

duration  of  collection  was  299.5  ±  1  K  with  confidence  derived  from  the  instrument 
performance  and  varied  by  3.6  K  during  the  experiment.  It  should  also  be  noted  that 
the  temperature  sensors  were  not  in-path,  but  were  located  near  exterior  doors  of  the 
building  within  5  meters  of  each  end  of  the  path.  The  average  outdoor  temperature 
was  about  290  ±  2  K,  and  may  have  contributed  to  the  lower  measured  meteorological 
temperature. 

Concentration 

Concentration  can  be  determined  from  the  intercept  of  Eq.  (15)  and  the  data  in 
Fig.  12  or  by  the  depth  of  any  rotational  feature.  Because  each  line  is  resolvable, 
concentration  will  be  calculated  for  each  line,  with  the  average  result  giving  the  final 
estimate  for  concentration  using  Eq.  (13).  The  resulting  concentrations  for  each 
transition  in  the  pP  and  PQ  branches  are  shown  in  Fig.  13. 

The  uncertainty  in  concentration  for  a  specific  line  is  0.1%.  The  standard  devia¬ 
tion  from  line  to  line  is  greater  than  0.8%.  The  best  estimate  for  oxygen  concentration 
is  4.82  ±  0.04  x  1018  molecules  /  cm3.  The  weather  station  provides  a  somewhat 
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Figure  13.  Estimation  of  concentration  from  each  rotational  line.  The 
lightly  dashed  line  gives  the  confidence  bounds  from  the  meteorological 
equipment,  the  heavily  dashed-dotted  line  gives  the  confidence  bounds 
from  the  concentration  estimate  using  the  y-intercept  method,  and  the 
heavily  dashed  line  gives  the  confidence  95  percent  bounds  for  concen¬ 
tration  using  the  fitted  areas  method. 


larger  estimate,  4.91  ±  0.16  x  1018  molecules  /  cm3,  with  a  larger  uncertainty.  Cal¬ 
culating  concentration  using  Eq.  (15)  yielded  results  that  agree  with  the  line  by  line 
method  giving  4.82  ±  0.05  x  1018  molecules  /  cm3. 

Pressure 

Determining  atmospheric  pressure  from  the  spectral  data  is  less  precise  than  con¬ 
centration  and  temperature.  The  high  resolution  spectra  exists  with  sufficient  signal- 
to-noise  to  distinguish  the  Lorentzian  component  of  the  Voigt  lineshape.  Fitting  a 
Voigt  lineshape  to  each  spectrally  isolated  line  provides  estimates  of  the  Lorentzian, 
Az/£,  and  the  Doppler,  Az/q,  full-width  half  maxima.  The  results  of  the  pressure 
broadened  lineshape,  with  fixed  and  unfixed  Doppler  widths,  along  with  the  expected 
linewidth  using  meteorological  data  and  the  H1TRAN  database  are  all  illustrated  in 
Fig.  14. 

There  is  a  strong  rotational  dependence  to  the  pressure  broadening  rates  due  to 
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Figure  14.  Estimation  of  pressure  from  meteorological  data  and  the 
HITRAN  database  (triangles),  from  the  fitted  Voigt  decomposition  by 
fixing  the  Doppler  width  (circles),  and  by  fitting  to  the  Doppler  and 
Lorentian  width  simultaneously  (diamonds). 


the  inelastic  energy  transfer.  [62]  Voigt  fits  to  individual  lines,  as  shown  in  Fig.  12, 
resulted  in  Lorentzian  widths  with  an  average  fit  uncertainty  of  ±  0.2%  for  each  line. 
Using  the  distribution  of  pressures  for  individual  features,  pressure  was  determined  as 
0.975  ±  0.04  atm  while  the  meteorological  instruments  reported  0.977  ±  0.001  atm. 
When  temperature  is  unconstrained,  the  Doppler  widths  are  reasonable  and  vary  by 
±  5.1%.  This  variation  is  demonstrated  in  Fig.  15. 

By  not  assuming  a  fixed  Doppler  width,  a  small  reduction  of  performance  in 
spectral  fitting  results  for  pressure  broadening,  and  decreased  fitted  temperature  con¬ 
fidence  boundaries  by  about  0.2  K.  For  additional  comparison,  pressure  broadening 
widths  calculated  from  data  collected  using  a  Bomem  DA-8  Fourier  Transform  Spec- 
tromieter  (FTS)  [62]  is  included  in  the  figure.  Spectra  in  Pope’s  data  was  collected 
at  0.032  cm-1  resolution  and  had  a  signal-to-noise  ratio  of  about  10  near  the  deepest 
rotational  lines.  The  results  compare  favorably  with  the  data  collected  here. 

Extending  collections  to  1  km  outdoor  open  paths  is  achieved  by  keeping  the 
same  general  geometry  but  moving  the  turning  mirror  farther  away  from  the  transmit 
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Figure  15.  Estimation  of  Doppler  widths  derived  from  NLLS  fitting 
(circles)  and  from  assuming  T  =  298.3  K  (triangles). 


and  receive  telescopes  and  realigning  until  the  best  signal  is  obtained.  1  km  paths 
produced  deeper  absorption  features  as  demonstrated  by  the  comparison  of  a  100  m 
path  to  a  1  km  path  in  Fig.  16. 

Implementation  of  longer  path  collections  make  the  system  more  sensitive  to  tele¬ 
scope  jitter  and  turbulence  shown  by  the  larger  standard  deviation  about  the  baseline 
of  A  A  fa  0.012,  more  than  seven  times  larger  than  for  the  100  meter  path  spectrum 
with  the  equivalent  system  settings.  The  absorbance  between  rotational  pairs  is 
clearly  nonzero  over  longer  paths  and  these  spectral  areas  must  be  excluded  from  the 
baseline  fitting  process.  Some  of  the  deeper  absorption  features  become  nearly  opaque 
over  the  200  kHz  width  of  the  tuning  laser,  and  the  measurement  of  absorbance  be¬ 
comes  limited  by  the  system  detection  noise  in  the  peak  of  a  long  path  absorbance 
feature. 

Monocular  Passive  Ranging 

The  process  of  viewing  bright  distant  sources  and  using  the  depth  of  absorption 
bands  to  give  an  estimate  of  range  to  the  source  is  referred  to  as  monocular  passive 
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Figure  16.  A  comparison  of  spectra  collected  at  100  m  and  1  km. 

ranging  (MPR)  and  can  be  used  with  a  variety  of  atmospheric  species.  Oxygen 
is  an  excellent  candidate  for  MPR,  being  spectrally  isolated  and  having  a  stable 
atmospheric  concentration.  [38]  The  TDLAS  instrument  was  deployed  simultaneously 
with  a  non-imaging  FTIR  spectrometer  to  validate  the  MPR  technique. 

An  example  FTIR  spectrum  of  the  oxygen  X to  transition  at  1  cm-1 

resolution  with  4,000  coadded  interferograms  is  demonstrated  in  Fig.  17a. 

The  FTIR  instrument  is  fast  framing  with  the  ability  to  produce  10  interferro- 
grams  per  second  at  its  maximum  resolution  of  1  cm-1.  The  spectral  region  with 
no  absorption  features  was  used  to  fit  a  cubic  baseline  through  the  X-b  transition  in 
order  to  give  the  spectrum  in  absorbance  as  shown  in  Fig.  17b.  The  lower  spectral 
resolution  of  the  FTIR  spectra  is  apparent  in  the  R-branch,  and  even  P-branch  ro¬ 
tational  lines  are  not  fully  resolved.  Baseline  noise  is  also  significantly  higher  than 
TDLAS  spectra. 

In  order  to  compare  performance  metrics  between  the  two  instruments,  statistics 
concerning  the  accuracy  of  integrated  absorbance  are  compared  as  a  function  of  col¬ 
lection  time.  The  uncertainty  in  absorbance  from  the  fully  integrated  spectra  as  a 
function  of  the  number  of  scans  is  shown  in  Fig.  18a  along  with  a  single  TDLAS  scan 
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Figure  17.  An  example  spectrum  produced  from  4,000  co-added  FTIR 
scans,  before  (a)  and  after  (b)  the  baseline  correction.  The  points  used 
for  fitting  a  baseline  to  raw  FTIR  data  are  shown  by  asterisks. 


having  an  absorbance  converges  to  5.915  with  a  standard  deviation  of  0.046. 

The  FTIR  spectra  exhibit  an  integrated  absorbance  of  5. 2-6. 2,  converging  to  5.784 
with  a  standard  deviation  of  0.252.  The  smaller  variance  in  TDLAS  absorbance  is 
shown  by  dashed  horizontal  lines  in  Fig.  18a.  The  standard  deviation  of  integrated 
absorbance  as  a  function  the  number  of  co-adds  in  Fig.  18b  give  the  confidence 
bounds  found  in  Fig.  18a  for  FTIR  data  points  while  overall  confidence  for  statistical 
performance  of  FITR  integrated  absorbance  are  shown  by  curved  dotted  lines.  It 
should  be  noted  that  the  data  shown  here  uses  both  the  P  and  R  branches,  but  the 
final  TDLAS  implementation  can  allow  for  a  single  branch  to  be  used  for  integrated 
absorbance  with  similar  performance  outcomes. 

The  time  required  to  obtain  a  single  FTIR  spectrum  (0.1  s)  is  shorter  than  a 
TDLAS  scan  (16  s).  The  standard  deviation  of  integrated  absorption  is  characterized 
in  terms  of  signal-to-noise  versus  collection  time  in  Fig.  19. 

The  raw  collected  baseline  defines  the  signal  and  the  noise  is  defined  by  the  RMS 
variation  about  that  baseline.  The  signal-to-noise  from  each  instrument  using  this  def¬ 
inition  results  in  similar  performance.  Collection  times  given  can  be  greatly  reduced 
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Figure  18.  Demonstration  of  the  statistical  standard  deviation  of  in¬ 
tegrated  absorption  (a)  as  a  function  of  the  number  collected  spectra. 
Statistics  from  the  FTIR  data  were  used  to  assign  confidence  bounds  for 
integrated  absorption  for  FTIR  collections  (b)  and  expected  performance 
of  FTIR  collections  as  a  function  of  scan  time.  Results  for  integrated  ab¬ 
sorption  are  shown  with  the  horizontal  line  with  one  standard  deviation 
of  confidence  lines  for  a  consecutively  long-running  TDLAS  collection 
are  also  shown  in  (a). 


Figure  19.  Demonstration  of  (dots)  FTIR  and  (open  circles)  TDLAS 
signal-to-noise  as  a  function  of  collection  time. 
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Figure  20.  The  location  of  the  potassium  D\  line  between  the  A'”  =31 
and  33  lines  in  the  molecular  oxygen  P-branch. 


for  the  TDLAS  instrument  if  the  lineshape  of  absorption  features  is  not  desired,  at  the 
expense  of  lost  information  regarding  pressure  and  possibly  temperature.  The  TD¬ 
LAS  instrument  achieves  the  same  performance  in  determining  the  total  absorbance 
for  the  same  acquisition  time,  while  adding  detailed  information  regarding  lineshape 
and  rotational  distribution. 

DPAL  Atmospheric  Transmission 

A  final  experiment  was  performed  to  demonstrate  the  application  of  this  TDLAS 
device  for  field  studies  of  laser  transmission.  The  potassium  variant  of  the  Diode 
Pumped  Alkali  Laser  (DPAL)  operates  on  the  D\  line  near  12985.185  cm-1,  [82]  in 
the  high  rotational  limit  of  the  O2  (X-b)  (0,0)  P  branch,  as  shown  in  Fig.  20. 

The  spectra  were  acquired  over  a  150  meter  outdoor  path  during  nighttime  con¬ 
ditions  during  weather  periods  with  low  wind.  During  the  data  collection,  a  Triad 
Technologies  potassium  vapor  cell  with  zero  pressure  (no  buffer  gas)  was  placed  in 
the  optical  path  just  before  the  detector  on  the  receive  telescope.  With  the  cell  in 
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place,  collections  were  conducted  with  cell  temperatures  of  65  C  giving  an  approxi¬ 
mate  vapor  pressure  of  2  x  10~5  mbar.  Temperatures  were  controlled  with  a  Watlow 
temperature  controller  and  the  cell  was  encased  in  a  custom  aluminum  heat  block 
with  thermocouples  placed  in  bored  holes  in  the  aluminum  block  within  4  mm  of  the 
vapor  cell.  Fig.  20  demonstrates  that  the  potassium  D\  absorption  line  lies  nearly 
midway  between  the  P  branch  K"  —  31  and  K"  —  33  (V’=0,  v'=0)  oxygen  lines. 

The  LBLRTM  simulation  using  the  TDLAS  extracted  atmospheric  conditions  is 
also  shown.  The  K  DPAL  system  operates  at  high  pressure,  1-20  atm,  depending 
on  the  spectral  width  of  the  diode  pump  lasers.  The  collision  induced  broadening 
and  lineshifts  can  shift  the  K  D\  line  by  up  to  0.5  cm-1.  Also  note,  at  longer  path 
lengths  the  O2  (X-b)  hot  bands  have  rotational  lines  closer  to  the  K  D\  line  and  could 
degrade  laser  transmission.  Further  characterization  of  atmospheric  transmission  for 
the  K,  Rb,  and  Cs  variants  of  the  DPAL  is  currently  in  progress  using  the  current 
TDLAS  apparatus. 

Typical  spectra  collected  using  the  TDLAS  technique  compares  favorably  com¬ 
pared  to  the  MUMAS  technique  outlined  in  Arita,  et  al.  [7].  While  TDLAS  collection 
times  are  longer,  the  resulting  raw  spectrum  is  detailed  enough  to  easily  extrapolate 
concentration,  temperature,  and  pressure  of  observed  atmospheric  species  without 
extensive  post-processing.  Sensitivity  of  MUMAS  is  reported  to  be  0.01%  for  128  av¬ 
eraged  spectra,  while  the  TDLAS  technique  can  detect  0.2%  over  a  single  scan.  Mea¬ 
sured  values  for  concentration,  temperature,  and  pressure  for  MUMAS  collections  of 
the  02 (X-b)  is  reported  as  ±  2%  for  each  measurement,  while  TDLAS  measurements 
are  ±  0.7%,  ±  0.8%,  and  ±  4%  respectively. 
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Conclusions 


A  system  to  conduct  open-path  TDLAS  measurements  has  been  developed  to 
investigate  the  molecular  oxygen  A-band  A13E“  to  b transition  lines  near  760 
nm.  The  resolution  of  the  system  was  limited  to  laser  lincwidth  of  less  than  300 
kHz  and  was  sufficient  to  separate  Voigt  lineshapes  into  their  Lorentzian  and  Doppler 
components.  Analysis  of  the  data  has  shown  that  accurate  estimates  for  temperature, 
concentration,  and  pressure  can  be  found.  Temperature  was  determined  to  within 
1.3%,  concentration  to  1.6%,  and  atmospheric  pressure  to  less  than  10%.  System 
noise  is  dominated  by  some  combination  of  turbulence  in  the  atmospheric  path  and 
telescope  jitter.  The  laser  source  can  easily  be  changed  to  investigate  other  spectral 
regions.  When  compared  to  collection  times  using  FTIR  instruments,  signal-to-noise 
is  similar  however  the  TDLAS  instrument  can  record  spectral  lineshapes  of  absorption 
features.  The  TDLAS  instrument  has  been  deployed  to  laser  test  ranges,  operated  at 
temperatures  as  low  as  -10  C  for  several  days  without  performance  degradation.  Path 
lengths  have  been  extended  to  1  km  with  minimum  detectable  changes  in  absorbance 
of  AA  =  0.06. 
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V.  Open-path  atmospheric  transmission  for  the  diode 

pumped  cesium  laser. 


A  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  device  was  used  to  prop¬ 
agate  a  laser  over  outdoor  atmospheric  open-paths.  Water  vapor  absorption  lines 
near  the  Cs  Diode  Pumped  Alkali  Laser  (DPAL)  wavelength  of  895  nm  were  inves¬ 
tigated.  Temperature,  pressure,  and  water  vapor  concentration  were  determined  for 
150  meter  and  1  km  open-paths  with  statistical  errors  of  0.2%.  Comparison  with 
meteorological  instruments  yields  agreement  for  the  1  km  path  to  within  0.6%  for 
temperature,  3.7%  for  pressure  and  2.4%  for  concentration. 

Introduction 

The  diode  pumped  alkali  laser  (DPAL)  was  first  proposed  in  2001  as  an  alternative 
to  high  power,  diode  pumped  solid-state  lasers  [44]  [47].  Diode  laser  excitation  on  the 
D2  2S1/2  -  2  P3/2  transition  and  collisional  energy  transfer  to  the  spin-orbit  split  2P\/2 
state,  yields  lasing  on  the  D\  2P\/2  -  2Si/2  transition  in  potassium,  rubidium,  or 
cesium  vapor.  The  DPAL  system  offers  excellent  thermal  control,  good  beam  quality, 
and  a  quantum  efficiency  of  95-99%.  A  high  power,  >1  kW,  system  with  48%  optical- 
to-optical  efficiency  has  recently  been  reported.  [14] 

The  cesium  version  of  the  DPAL  system  has  received  considerable  attention.  [28, 
45,  63,  94,  96,  97,  98,  99,  100,  101,  102]  Slope  efficiencies  exceeding  80%  have  been 
demonstrated  [96],  exciplex  assisted  pumping  schemes  [63],  and  transverse  excitation 
with  an  unstable  resonator  have  been  explored  [102],  The  larger  spin-orbit  splitting  in 
cesium  generally  requires  a  molecular  collision  partner  to  accelerate  the  fine  structure 
mixing.  [100]  Rare  gases  are  also  usually  employed  at  pressures  of  1-10  atmospheres  to 
broaden  the  absorption  profile  and  better  match  the  diode  laser  spectral  distribution. 
[60,  61] 
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Interest  in  the  atmospheric  transmission  characteristics  of  these  alkali  lasers  has 
led  us  to  develop  of  a  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  device  for 
long  open-path,  >1  km,  field  deployment.  [65]  The  system  was  originally  designed  to 
study  the  potassium  DPAL  that  operates  in  a  narrow  window  between  the  rotational 
features  of  the  O2  X-b  (0,0)  transition  near  760  nm.  In  the  current  work,  we  extend 
the  TDLAS  system’s  functionality  to  the  Cs  DPAL  at  wavelengths  near  894  nm. 
There  are  several  moderate  strength  water  lines  and  very  weak  CO2  and  OH  lines  in 
the  vicinity  of  the  Cs  DPAL  wavelength.  [74] 

TDLAS  has  previously  been  used  to  explore  atmospheric  water  vapor  absorption 
features  for  a  number  of  applications.  Water  vapor  detection  in  the  Martian  atmo¬ 
sphere  has  been  demonstrated  for  a  120  cm  open-path  and  a  10”4  minimum  detectable 
absorption  during  a  single  scan.  [9,  25]  The  measurement  of  water  vapor  with  an  un¬ 
certainty  of  about  1.14%,  was  demonstrated  with  Voigt  fits  having  residuals  only  1  - 
3  times  larger  than  the  noise.  Other  groups  use  diode  lasers  to  investigate  open-path 
water  vapor  on  airborne  LIAVs  [76]  and  manned  [21]  platforms.  For  example,  a  1.39 
m  diode  laser  over  an  open  folded  50  cm  absorption  path  placed  in  the  free-stream 
air  using  an  external  probe  with  a  signal-to-noise  ratio  of  3  and  a  minimum  mea- 
sureable  absorbance  of  10-5  Hz-1/2  was  demonstrated  on  a  NASA  P3B  aircraft.  [76] 
A  system  with  a  diode  operating  near  1.4  m  was  deployed  on  a  DC-8-72  between  a 
quartz  fuselage  window  and  a  retro-reflector  mounted  on  the  plane’s  outboard  engine 
nacelle.  [21] 

Several  balloon-borne  instruments  have  measured  water  vapor  concentrations  to 
within  10%  over  the  upper  troposphere  and  lower  stratosphere.  [27,  26,  24,  23,  37]  One 
device  employs  a  56  meter  Herriott  cell.  [23]  Atmospheric  sampling  for  altitudes  of  1 
-  31  km  yields  mixing  ratio  with  errors  of  less  than  10%.  [26]  Absorbance  detection 
limits  of  10~5  have  been  demonstrated.  [23]  A  more  recent  article  outlines  a  very 
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Figure  21.  Diagram  of  TDLAS  system. 


similar  device  with  an  updated  Herriott  absorption  cell  having  with  different  paths 
for  water  vapor  (36  meter  path)  and  methane  (74  meter  path).  [37] 

Finally,  diode  lasers  operating  near  1.37  m  are  used  to  determine  transpiration 
rates  of  plant  leaves  giving  water  vapor  concentrations  to  within  5%.  [42]  A  more 
recent  article  employs  similar  methods  investigating  stronger  vibrational  bands  near 
2.7  m  and  also  resolves  water  vapor  spatially  by  moving  the  absorption  path  using  a 
mechanical  scanning  platform.  [92]  Using  diode  lasers  over  open-paths  is  a  relatively 
low  risk  method  that  can  investigate  more  than  a  single  spatial  location  and  can  have 
detection  limits  in  the  neighborhood  of  10-5  in  absorbance. 

TDLAS  Device 

The  tunable  diode  laser  absorption  spectroscopy  instrument  is  shown  schemati¬ 
cally  in  Fig.  21. 

A  detailed  description  of  the  device  configuration  and  method  for  collecting  spec¬ 
tra  have  previously  been  discussed.  [65,  64]  The  tunable  diode  is  a  New  Focus  Velocity 
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laser  model  6318  with  a  10  mWatt  maximum  output,  a  890-910  nm  tunable  range,  and 
less  than  a  300  kHz  linewidth.  The  laser  source  near  900  nm  can  easily  be  changed 
to  investigate  other  spectral  regions  of  interest  with  minimal  changes  to  the  over¬ 
all  system.  After  amplitude  modulation  by  a  Thor  Labs  EO-AM-NR-C1  free  space 
modulator  at  102  kHz  (lock-in  maximum  useable  modulation  rate),  the  laser  is  fiber 
coupled  and  expanded  to  nearly  fill  a  military  grade  12.5  inch  RC  Optical  Ritchey- 
Chrtien  transmit  telescope.  The  laser  beam  is  directed  across  an  open-path  and  then 
received  by  a  second  identical  telescope.  Two  Thor  Labs  PDA100A  Si  detectors  are 
used  to  measure  the  reference  (or  sent)  and  transmitted  (or  received)  signals,  and  are 
analyzed  by  Stanford  Research  Systems  SR850  dual-phase  lock-ins  and  recorded  by 
a  National  Instruments  LTSB-6251  DAQ.  Pellicle  beam  splitters  of  2  nr  thickness  are 
used  to  minimize  etalon  effects.  The  reference  intensity  is  measured  late  in  the  optical 
train  (immediately  prior  to  the  telescope  beam  expansion)  and  can  be  attenuated  to 
balance  transmitted  and  received  signals.  A  High  Finesse  WSLI-2  wavemeter  is  used 
to  determine  the  frequency  axis  and  is  calibrated  with  a  2  MHz  fiber-coupled  SIOS 
SL-03  frequency-stabilized  HeNe  laser  to  achieve  a  better  than  10  MHz  accuracy.  All 
hardware  is  connected  to  and  controlled  by  a  PC  using  MATLAB®.  During  an  experi¬ 
ment,  the  laser  is  finely  tuned  over  an  approximately  0.1  nm  range  by  driving  a  piezo 
to  change  cavity  length,  and  then  coarsely  tuned  to  the  next  spectral  region.  The 
process  is  repeated  for  nearly  200  free  spectral  ranges  to  cover  the  available  20  nm 
spectral  scan.  Calibration  of  the  frequency  axis,  baseline  removal,  and  further  spec¬ 
tral  processing  have  also  been  reported.  [65,  64]  Noise  when  the  system  is  deployed 
outdoors  is  generally  dominated  by  telescope  jitter  and  atmospheric  turbulence  over 
the  path.  For  shorter  paths  near  100  meters,  these  effects  are  usually  small  and  less 
than  1%  of  the  signal.  The  maximum  open-path  distance  is  about  2  km,  a  limitation 
imposed  by  the  Quickset  QPT-130  ruggedized  pan  and  tilt  mounts. 
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Figure  22.  Data  collection  geometry  for  the  150  m  (top)  and  1  km 
(bottom)  paths  showing  the  position  of  the  TDLAS  transmit  and  receive 
telescopes  to  the  east  (A)  and  the  turning  mirror  west  side  of  the  open- 
path  (B).  Buildings  are  depicted  in  black  while  pavement  is  depicted  in 
grey.  Wind  prevailed  in  the  northerly  direction  for  both  cases. 

Collections  were  performed  over  a  short  150  meter  and  a  longer  1  km  open  path 
during  fair  weather  conditions  with  site  geometries  shown  in  Fig.  22. 

The  send  and  receive  telescopes  were  separated  by  1  meter,  each  viewing  a  6  inch 
flat  mirror  at  1  meter  above  the  ground.  All  collections  were  performed  at  night  to 
operate  with  minimal  turbulence.  Two  NIST-certified  Davis  Vantage  Pro2  weather 
sensors  were  placed  near  the  telescope  and  mirror  for  the  150  m  path.  Temperature, 
pressure,  and  relative  humidity  were  recorded  at  15  minute  intervals  with  95%  con¬ 
fidence  bounds  of  1  C,  2  hPa,  and  6%,  respectively.  The  150  m  total  path  length 
spectra  were  recorded  on  September  8th,  2011  over  cut  grass  with  average  meteoro¬ 
logical  conditions  of  12.3  C,  996.7  hPa,  and  87.2%  relative  humidity.  The  1  km  path 
length  spectra  were  collected  on  November  24th,  2011  over  a  paved  runway  surface 
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Figure  23.  An  example  spectrum  demonstrating  a)  the  /  and  I0  detector 
response  and  b)  the  A=  -log(/t/fo)  for  water  vapor  over  tunable  laser 
range  and  the  typical  baseline  curvature. 

with  average  conditions  of  4.5  C,  994.3  hPa,  and  90.7%  relative  humidity  recorded  by 
a  single  Davis  Vantage  Pro2  sensor  in  the  center  of  the  folded  path. 

Results 

A  typical  absorption  spectrum  from  11,000  to  11,250  cm-1  for  an  atmospheric 
path  length  of  100  m  is  shown  in  Fig.  23a. 

The  periodic  modulation  in  the  send  reference  intensity,  /0,  occurs  throughout  the 
171  free  spectral  range  scans  of  the  diode  laser  and  comes  from  etalon  effects  from 
the  instrument.  The  laser  scan  rate  is  0.002  nrn/s  and  signal  is  sampled  at  512  Hz. 
While  the  figure  demonstrates  a  single  sweep,  sweeps  can  be  averaged  to  improve  the 
signal-to-noise  ratio.  The  received  signal,  /t,  has  a  slightly  different  response,  so  that 
the  observed  spectral  absorbance, 

A(v)  =  -  In  b(v)  =  b{v)  Y  Vi{v)NiL  (16) 

includes  a  frequency  dependent  cubic  baseline,  6,  in  addition  to  the  absorption  cross- 
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Figure  24.  Example  absorption  spectra  for:  150  m  (black)  and  1  km 
(grey)  open-path  outdoor  collections  with  low  pressure  Cs  D\  hyperfine 
lines  located  in  path. 


section,  cr*(z>) ,  for  each  atmospheric  constituent,  i,  with  concentration,  Nt.  All  of 
the  spectral  features  in  Fig.  23b  are  assigned  to  water  vapor  and  exhibit  lineshapes 
limited  only  by  the  atmospheric  pressure  broadening.  Fig.  24  demonstrates  the 
spectra  from  the  150  meter  and  1  km  collections  in  the  vicinity  of  the  Cs  D\  2S\/2 
-2P\/2  line,  with  the  baseline  removed. 

All  atmospheric  absorption  lines  in  the  spectrum  are  a  result  of  water  vapor,  and 
are  assigned  using  the  vibration  quantum  numbers  of  the  upper  level  {vV ,  v2',  v3')  and 
the  ground  vibrational  (0,  0,  0)  level  is  common  for  all  the  features.  Rotational  states 
are  labeled  with  the  quantum  numbers  [J',  K'A,  K'C]  —  [J”,  K”  A,  K”C ]  for  the  upper 
and  lower  states,  respectively.  A  blended  pair  of  lines  (0,  0,  3)  [10, 0, 10]  —  [9,  0,  9]  and 
(0,  0,  3)  [10, 1, 10]  —  [9, 1,  9]  make  up  the  first  feature  near  11,174  cm-1  while  another 
three  lines  near  11,175  cm-1  create  a  broader  shape  just  to  the  right.  The  two 
features  on  either  side  of  the  Cs  Di  emission  line  are  the  (0,0,3)  [9,4,6]  —  [8,4,5] 
on  the  left  and  the  (0,0,3)  [7,2,5]  —  [6,2,4]  blended  with  the  (0,0,3)  [10,5,6]  — 
[9,  5,  5]  on  the  right.  There  are  additional  smaller  water  vapor  lines  in  this  spectral 
area,  but  they  are  not  immediately  visible  over  these  shorter  path  lengths.  The  150 
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m  collection  has  a  baseline  signal  RMS  from  11,184  to  11,185  cm-1  of  0.00169  in 
absorbance  with  a  peak  signal  located  on  the  blended  pair  of  water  lines  at  11,182.36 
cm-1  of  0.1238  giving  a  signal-to-noise  ratio  of  about  73.2.  The  1  km  collection 
exhibited  a  baseline  signal  RMS  of  0.0080  in  absorbance  with  the  same  water  vapor 
peak  absorbance  of  0.4437  giving  a  signal-to-noise  ratio  of  54.9.  As  path  lengths 
increase  noise  from  path  turbulence  and  system  pointing  jitter  increase  at  a  rate 
greater  than  the  absorption  signal  increases.  During  the  1  km  collection  a  Scintec 
BLS900  turbulence  measurement  device  was  setup  beside  the  TDLAS  device  with 
the  transmitter  placed  by  the  turning  mirror  and  the  receiver  between  the  transmit 
and  receive  telescopes  giving  a  500  m  path.  The  BLS900  device  is  limited  to  C2 
measurements  between  3  x  10~14  and  3  x  10~10  for  500  m  path  collections.  During 
the  experiment  values  for  C2  were  below  the  lower  threshold.  Even  assuming  a  high 
value  of  1  x  10~14  for  turbulence,  RMS  beam  wander  from  turbulence  over  a  1  km 
path  is  less  than  12  mm.  During  the  same  experiment  the  alignment  laser  was  seen 
to  move  on  a  target  board  by  over  4  cm  with  movements  that  correlated  in  time  with 
wind  gusts.  While  both  turbulence  and  platform  jitter  both  contribute  to  noise,  at  1 
km  paths  with  turbulence  under  1  x  10~14  platform  jitter  is  likely  to  be  the  primary 
cause  of  baseline  noise. 

A  low  pressure  (no  buffer  gas)  Cs  vapor  cell  inside  a  heater  block  at  30  C  was  placed 
in  the  path  between  the  receive  telescope  and  the  detector  to  provide  a  reference  for 
the  Cs  DPAL  wavelength.  The  strong,  narrow  absorption  features  near  11,178  cm-1 
in  Fig.  24  are  due  to  the  Di  2S\/2  -2P\/2  transition  in  cesium.  The  hyperhne  structure 
of  this  transition  is  easily  resolved.  The  hyperhne  splittings  are  approximately  9.192 
and  1.168  GHz  for  the  ground  and  excited  states,  respectively,  as  shown  in  Fig.  25. 

A  best  estimate  for  the  atmospheric  temperature,  pressure,  and  water  vapor  con¬ 
centration  were  obtained  using  the  Line-by-Line  Radiative  Transfer  Model  (LBLRTM) 
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Figure  25.  Structure  of  Cs  including  4  hyperfine  lines  that  make  up  the 
lasing  emission  on  the  D\  line. 


57 


[19]  in  conjunction  with  the  HITRAN  database  [74]  and  knowing  the  overall  path 
length,  a  non-linear  least  squares  fitting  process  was  used  by  simulating  an  atmo¬ 
spheric  spectra  defined  for  a  homogenous  atmosphere  with  a  single  temperature, 
pressure,  and  water  vapor  concentration  as  fit  parameters.  The  cubic  baseline  was 
removed  before  the  fitting  process  using  only  data  points  that  lie  outside  of  the  Cs 
hyperhne  lines  and  the  water  vapor  lines  present  in  the  spectrum.  The  iterative  fitting 
process  also  includes  the  95%  confidence  interval  for  each  fitted  parameter  based  on 
results  from  each  of  the  iterative  fits,  giving  the  final  estimates  for  the  average  path 
temperature,  pressure,  and  water  vapor  concentration.  Fig.  26a  illustrates  the  fit 
quality  for  the  150  nr  path  with  residuals  in  6b  resulting  in  an  estimated  temperature 
of  15.5  C  (288.6±0.15  K),  a  pressure  of  993.8±0.9  hPa,  and  a  water  vapor  concentra¬ 
tion  of  3.096  ±  0.004  x  1 0 1 '  molecules  per  cm3.  For  comparison,  the  meteorological 
conditions  during  the  test  are  reported  in  Figs.  27a,  27b,  and  27c. 

The  agreement  is  poorer  than  the  statistical  error  from  the  spectral  fitting  and 
agree  within  2.4%.  However,  the  weather  station  data  indicates  variations  within 
test  of  up  to  5%.  The  fit  results  for  the  1  km  path  yield  a  temperature  of  6.1  C 
(279.2±0.20  K),  a  pressure  of  958.3±1.5  hPa,  and  a  water  vapor  concentration  of 
1.940±0.004  x  101'  molecules  per  cm3,  all  agreeing  with  the  meteorological  instru¬ 
ments  to  within  3.7%  or  less  with  the  fit  shown  in  Fig.  26c  with  residuals  in  Fig.  26d. 
The  meteorological  data  collected  during  the  1  km  experiment  is  shown  in  Figs.  27d, 
27e  and  27f.  The  particularly  poor  estimate  for  pressure  for  the  long  path  is  easy 
to  understand  considering  the  fit  to  pressure  is  especially  sensitive  to  the  Lorentzian 
(wing)  portion  of  the  Voigt,  making  it  a  parameter  especially  sensitive  to  the  overall 
signal-to-noise  ratio.  Structured  residuals  are  evident  resulting  from  minor  wavemeter 
errors  expected  because  of  the  ridged  line  center  constraints  used  in  the  fit.  Includ¬ 
ing  linecenters  in  the  fit  removes  structure  in  the  residuals  leaving  only  noise  and 
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Figure  26.  Demonstration  of  a)  150  m  spectrum  with  NLLS  fit  and  b) 
residuals,  c)  1  km  spectrum  with  NLLS  fit  and  d)  residuals. 


Figure  27.  Variations  in  temperature,  relative  humidity,  and  pressure 
during  experimental  collections  for  the  150  m  path  (a,b,c)  with  a  sensor 
located  near  the  transmit  telescope  (black)  and  another  near  the  turning 
mirror  (grey)  and  1  km  path  lengths  (d,e,f)  with  a  single  sensor  placed 
in  the  center  of  the  path. 
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reducing  the  residual  RMS  by  about  20%.  Adding  linecenters  as  a  fitting  parameter 
increases  fitting  complexity  without  improving  overall  fitting  results,  and  results  from 
this  fitting  method  are  not  used. 

Discussion 

Having  anchored  LBLRTM  simulations  to  observed  TDLAS  spectra  over  exam¬ 
ple  paths,  atmospheric  simulations  can  be  compared  to  Cs  laser  lineshapes.  DPAL 
systems  require  1-10  atmospheres  of  added  buffer  gases  for  efficient  operation.  High 
power  diode  bars  and  stacks  exhibit  broad  spectral  lineshapes  relative  to  atomic,  gas 
phase  absorption  features.  Pressure  broadening  of  the  alkali  D\  and  D2  lines  occur 
at  rates  of  «  20  MHz/Torr.  For  example,  argon  broadens  the  Cs  D2  line  at  18.31 
MHz/Torr.  [60]  Furthermore,  a  buffer  gas  is  required  to  induce  transitions  from  the 
pumped  2P3/2  level  to  the  2Pi/2  level.  The  foie  structure  splitting  in  Cs  is  relatively 
large  (10.4  GHz),  as  shown  in  Fig.  25.  Molecular  collision  partners  such  as  methane 
or  ethane  are  typically  employed  to  increase  the  rate  of  collisional  energy  transfer.  [14] 
Several  Cs  lasers  have  been  demonstrated  with  powers  of  up  to  49  W  [99]  with  either 
helium  or  argon  and  ethane  buffer  gases  at  total  pressures  of  about  1  atmosphere. 
Fig.  28  shows  a  demonstration  of  a  1  km  homogenous  atmosphere  at  0  C  (273.15  K), 
1  atm  (1013.25  hPa),  and  50%  RH  (8.10  x  1017  molecules  per  cm3)  at  ground  level 
alongside  three  examples  of  Cs  gain  profles,  one  with  only  100  Torr  C2H6  and  0  Torr 
He,  one  with  100  Torr  C2H6  and  760  Torr  He,  and  one  with  100  Torr  C2H6  and  7,600 
Torr  He.  [61,  59] 

The  simulation  with  only  100  Torr  C2H6  shows  some  hyperfine  structure.  Not  all 
four  hyperfine  lines  are  resolvable,  but  each  pair  of  hyperfine  lines  are  distinguishable 
and  only  the  far  wings  of  the  gain  profile  lie  over  any  water  vapor  absorption  lines. 
The  Voigt  lineshape  has  a  Doppler  component  of  width  0.032  cm-1  and  a  Lorentzian 
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Figure  28.  Simulation  of  1  km  horizontal  atmospheric  path  (grey)  with 
the  simulated  gain  profiles  of  a  Cs  vapor  gain  cell  with  0  atm  He  (solid 
black),  1  atm  He  (dashed  black)  and  10  atm  He  (short  dash  back),  each 
with  a  buffer  gas  mixture  including  100  Torr  C2H6.  Gain  profiles  are 
normalized  to  one. 

width  of  0.194  cm-1  for  this  low  pressure  simulation.  With  1  atm  of  He,  the  hyperhne 
lines  are  completely  washed  out  and  the  Cs  absorption  profile  begins  to  encompass 
the  (0,0,3)  [9,4,6]  -  [8,4,5],  (0,0,3)  [7,2,5]  -  [6,2,4],  and  (0,0,3)  [10,5,6]  -  [9,5,5] 
water  vapor  lines.  With  10  atm  of  added  to  the  100  Torr  C2H6,  the  gain  profile 
completely  surrounds  many  major  water  lines  and  is  visibly  bine  shifted.  If  the 
same  amount  of  laser  power  is  distributed  across  each  Cs  absorption  profile,  the 
emission  line  from  a  high  pressure  Cs  gain  cell  would  encounter  significant  water 
vapor  absorption  features.  This  absorption  would  also  contribute  to  higher  thermal 
blooming  for  high  energy  lasers.  While  a  lower  pressure  (1  atmosphere)  variant  of  the 
Cs  DPAL  would  be  preferable  for  atmospheric  propagation,  the  demand  for  narrow 
banding  the  pump  diode  bars  requires  emerging  technologies  such  as  external  cavity 
volume  Bragg  gratings  (VBG).  Narrow  banding  of  high  power  diode  bars  to  pump  the 
Rb  DPAL  system  have  recently  been  developed  with  linewidths  of  18  GHz,  enabling 
the  1  atmosphere  variant  of  the  DPAL  laser.  [36] 

It  should  also  be  noted  that  observed  spectra  of  highly  pressure  broadened  Cs 
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Figure  29.  Simulation  of  1  km  horizontal  atmospheric  path  (thin  black) 
to  simulated  Cs  cell  with  100  Torr  C2H6  (grey),  with  3  atm  of  Ar  (thick 
black)  and  a  laboratory  Cs  cell  with  3  atm  of  Ar  (x). 

exhibit  asymmetric  line  shapes  are  not  modeled  here.  The  theoretical  D\  lineshape 
for  a  Cs  gain  cell  with  3  atm  of  argon  is  depicted  in  Fig.  29  alongside  a  symmetric  Cs 
absorption  spectrum  [52],  On  the  blue  side  a  bump  in  the  absorption  profile  is  seen 
experimentally  but  not  modeled  here. 

ft  should  be  noted  that  this  spectral  asymmetry  could  result  in  Cs  DPAL  with  a 
spectrally  wider  emission  line  than  indicated  in  the  figures  above.  The  DPAL  gain 
is  high  with  a  low  cavity  Q,  unstable  resonators  are  likely  required  for  high  power 
operation.  Rather  than  consider  line  narrowing  associated  with  the  laser  resonator, 
we  use  the  small  signal  gain  lineshape  as  a  conservative  estimate  of  the  laser  linewidth. 

A  scenario  was  created  with  a  slant  range  of  100  km  with  both  starting  and  ending 
points  at  10  km  altitudes,  a  geometry  having  a  low  point  altitude  in  the  center  of 
about  9.8  km  clue  to  earth  curvature.  LBLRTM  was  used  to  compute  optical  depth 
and  configured  to  calculate  Voigt  lineshapes  to  within  25  cm-1  of  their  respective 
line  centers,  to  include  all  continuum,  and  using  an  atmospheric  profile  created  by 
LBLRTM  with  the  1976  US  Standard  model.  No  aerosol  effects  were  included  in 
the  simulation.  Fig.  30  shows  a  close-up  highlighting  the  100  km  high  altitude  path 
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Figure  30.  Simulations  of  a  10  km  altitude  atmospheric  path  with  a  100 
km  slant  range  (black  line),  a  1  km  horizontal  atmospheric  path  (dotted 
line),  and  a  100  Torr  C2H6  gain  cell  (grey). 

compared  to  the  1  km  low  altitude  atmospheric  path  and  low  pressure  Cs  gain  profile 
-  both  shown  also  in  Fig.  29. 

Molecular  scattering  over  the  long  path  is  obvious  in  the  baseline,  and  with  a 
low  pressure  of  268.5  hPa  (~  0.26  ATM)  the  lines  are  visibly  narrower  than  the  low 
altitude  case. 

Conclusion 

The  TDLAS  system  is  shown  to  collect  information  for  molecular  absorption  over 
atmospheric  paths  from  hundreds  to  thousands  of  meters.  Accuracy  in  the  wavenum¬ 
ber  axis  is  limited  to  the  spectral  reference  used  and  is  approximately  +/-  10  MHz. 
Noise  is  primarily  dominated  from  effects  of  telescope  and  mirror  jitter  with  RMS 
jitter  for  150  m  and  1  km  paths  of  0.00169  and  0.0080  giving  a  signal-to-noise  ratio 
of  73.2  and  54.9,  respectively,  in  absorbance.  The  high-quality  spectral  scans  from 
the  system  can  be  used  for  a  variety  of  tasks,  including  investigation  of  multi-species 
molecular  absorption,  estimation  of  temperature,  pressure,  and  concentration  for  wa¬ 
ter  vapor  with  results  to  within  2.4%  or  less  from  150  m  paths  and  3.7%  for  1  km 
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paths  of  meteorological  data,  all  using  a  nonlinear  least  squares  fit  with  LBLRTM. 
Overall  Cs  DPAL  atmospheric  transmission  will  be  excellent  over  atmospheres  with 
very  little  water  vapor.  Cs  DPALs  with  lower  pressure  gain  cell  configurations  will 
have  the  best  atmospheric  transmission  over  well-hydrated  atmospheric  paths  while 
Cs  DPALs  with  spectral  broad  emission  lines  will  suffer  from  water  vapor  absorption 
which  may  in  turn  cause  molecular  absorption  induced  thermal  blooming.  Asym¬ 
metry  of  highly  pressure  broadened  Cs  absorption  features  may  result  in  spectrally 
broader  Cs  DPALs. 
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VI.  Molecular  absorption  of  diode  pumped  alkali  lasers 


Molecular  oxygen  absorption  lines  near  the  potassium  (K),  and  water  vapor  ab¬ 
sorption  lines  near  the  rubidium  (Rb),  and  cesium  (Cs)  Diode  Pumped  Alkali  Lasers 
(DPAL)  at  wavelengths  near  770  nm,  795  nm,  and  895  nm,  respectively,  were  inves¬ 
tigated  using  the  Line-by-Line  Radiative  Transfer  Model  (LBLRTM)  with  the  High 
Energy  Laser  End-to-End  Simulation  (HELEEOS).  A  tunable  diode  laser  absorption 
spectroscopy  (TDLAS)  device  was  used  to  anchor  simulations  to  outdoor  150  m  at¬ 
mospheric  open-path  collections.  The  Cs  DPAL  operates  in  the  wing  of  water  vapor 
lines  and  exhibits  absorbances  similar  to  the  Chemical  Oxygen  Iodine  Laser  (COIL). 
The  Rb  DPAL  absorbance  is  lower  for  many  engagement  scenarios  with  A  =  0.002 

-  0.2.  The  K  DPAL  operates  between  the  rotational  lines  of  O2  to  blYi+  with 

absorbances  as  high  as  0.07  for  tactical  air-to-ground  missions. 

Introduction 

The  diode  pumped  alkali  laser  (DPAL)  was  first  proposed  in  2001  as  an  alternative 
to  high  power,  diode  pumped  solid-state  lasers.  [44]  Laser  excitation  on  the  D2  2<S'i/2 

-  2  P3/2  (pump)  transition  and  collisional  energy  transfer  to  the  spin-orbit  split  2P\/2 
state,  yields  lasing  on  the  Di  2Pi/2  -  2Si/2  (emission)  transition  in  potassium  (K), 
rubidium  (Rb),  or  cesium  (Cs)  vapor.  The  energy  level  structure  used  in  DPAL 
operation  for  K,  Rb,  and  Cs  are  similar  in  nature,  shown  in  Figs.  31,  32,  and  25, 
respectively,  giving  pump  and  laser  emission  locations  along  with  emission  hyperffiie 
lines. 

The  three  level  structure  they  each  share  results  in  high  efficiency  with  losses 
mostly  incurred  from  the  spin-orbit  relaxation  from  the  2P3//2  -  2Pi/2-  The  spectral 
width  of  the  pump  and  emission  lines  are  primarily  dependent  on  pressure  and  buffer 
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Figure  31.  Structure  of  potassium  (K)  DPAL  operation  from  pumping 
the  D2  line  near  766.7  nm,  with  spin-orbit  relaxation  from  the  42P3/2  to 
the  A2*Pi/2  levels  resulting  in  lasing  emission  from  the  D\  line  near  770.1 
nm  made  up  of  4  hyperfine  lines  closely  spaced. 
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Figure  32.  Structure  of  rubidium  (Rb)  DPAL  operation  from  pumping 
the  D2  line  near  780.2  nm,  with  spin-orbit  relaxation  from  the  52P3/2  to 
the  52P ! j2  levels  resulting  in  lasing  emission  from  the  D\  line  near  795.0 
nm  made  up  of  4  hyperfine  lines  with  moderate  spacing. 


gas.  First  demonstrated  with  Rubidium  (Rb)  [47],  DPAL  systems  offer  excellent 
thermal  control,  good  beam  quality,  and  a  quantum  efficiency  of  95-99%,  generating 
considerable  interest  within  the  laser  community.  [107] 

The  potassium  DPAL  lases  at  the  shortest  wavelength  at  770.1  nm  near  the  O2 
(X-b)  transition.  [95]  The  rubidium  DPAL  has  also  shown  successful  implementation 
[46,  80,  106,  103]  and  the  limitations  imposed  by  atom  cycle  rates  of  rubidium  have 
been  discussed.  [54]  The  cesium  version  of  the  DPAL  system  has  received  considerable 
attention.  [28,  96,  94,  98,  63,  99,  100,  102]  Recently,  a  high  power  flowing  system 
greater  than  1  kW  with  a  48%  optical-to-optical  efficiency  has  been  reported.  [14] 
Lithium  and  sodium  have  a  much  smaller  separation  between  the  D\  and  D-2  lines, 
making  them  less  feasible  to  pump  by  the  broad  spectral  emission  of  diode  stacks. 

Interest  in  the  atmospheric  transmission  characteristics  of  these  alkali  lasers  led 
to  the  development  of  a  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  device 
for  long  open-path,  >1  km,  field  deployment.  The  system  was  originally  designed  to 
study  the  potassium  DPAL  that  operates  in  a  narrow  window  between  the  rotational 
features  of  the  O2  X-b  (0,0)  transition  near  760  nm.  In  the  current  work  we  extend 
the  TDLAS  systems  functionality  to  the  Rb  DPAL  at  wavelengths  near  795  nm  and 
review  previous  work  for  K  [65]  and  Cs  [17].  There  are  several  moderate  strength 
water  lines  in  the  vicinity  of  the  Rb  and  Cs  DPAL  emission  lines  and  molecular 
oxygen  transition  near  the  K  DPAL  emission  line.  [73] 

Performance  differences  in  molecular  absorption  for  each  DPAL  variant  for  a  va¬ 
riety  of  weapons  engagement  scenarios  are  considered.  Simulations  are  validated 
with  open-path  TDLAS  spectra.  The  High  Energy  Laser  End-to-End  Operational 
Simulation  (HELEEOS)  is  used  to  compare  the  K,  Rb,  and  Cs  DPAL  systems  to 
the  Chemical  Oxygen  Iodine  Laser  (COIL).  Path  integrated  absorbance,  total  sys¬ 
tem  Strehl  ratio,  and  required  dwell  time  are  evaluated  for  long  range  “air-to-air”, 
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“air-to-ground”,  “ground-to-air”,  “low  altitude”  ocean,  and  “low  power”  jamming 
scenarios. 

Simulation  of  DPAL  Transmission 

High  energy  laser  simulations  like  HELEEOS  [29],  High  Energy  Laser  Consol¬ 
idated  Modeling  and  Engagement  Simulation  (HELCOMES)  [48],  and  HELSEEM 
(High-Energy  Laser  System  End-to-End  Model)  [70],  treat  laser  emission  using  a 
monochromatic  model.  The  wings  of  atmospheric  lineshapes  are  usually  considered, 
but  only  for  the  final  laser  wavelength  in  question.  For  the  present  work,  the  spec¬ 
tral  overlap  of  the  absorption  lineshapes  with  the  laser  small  signal  gain  is  used  to 
evaluate  the  path  integrated  absorbance.  The  Line-by-Line  Radiative  Transfer  Model 
(LBLRTM)  [19]  was  implemented  within  HELEEOS  3.0  BETA  to  estimate  molecu¬ 
lar  absorption  for  a  variety  of  scenarios.  LBLRTM  can  also  compute  a  variety  of 
other  effects,  like  continuum,  molecular  scattering,  aerosol  effects,  and  the  like,  but 
for  these  simulations  LBLRTM  was  compiled  and  driven  with  inputs  to  only  consider 
molecular  absorption.  Aerosol  absorption,  aerosol  scattering,  molecular  scattering, 
and  continuum  were  all  computed  using  HELEEOS  routines  while  LBLRTM  was 
passed  inputs  on  molecular  concentration,  temperature,  and  pressure  and  returned 
molecular  absorption  with  very  high  resolution  back  to  HELEEOS.  The  high  resolu¬ 
tion  molecular  absorption  profile  was  then  used  with  the  DPAL  laser  gain  profile  to 
determine  overall  molecular  absorption: 

OO 

f  A(p)G(y)dv 

A  =  L -  (17) 

/  G(v)dD 
o 

where  the  laser  lineshape,  G(u),  is  integrated  over  the  spectral  area  where  the  laser 
gain  profile,  typically  10  times  the  FWHM.  The  atmospheric  path  in  length  ab- 
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sorbance  is  defined  as, 


Aiy)  =  £  /  p(-)’  T(-))  •  Ni(z)dR  (18) 

i 

where  z(R)  is  the  altitude  as  a  fraction  of  slant  range,  R.  The  absorption  cross- 
section,  <7j,  for  each  atmospheric  constituent,  i,  is  dependent  on  frequency,  z>,  and 
pressure,  P,  and  temperature,  T.  The  absorption  depends  on  the  species  specific 
concentration,  Nt}  as  a  function  of  altitude,  z,  and  path  length,  R.  Overall  molecular 
absorbance,  A,  is  used  in  HELEEOS  as  an  input  to  compute  thermal  blooming. 
Predicted  atmospheric  absorbance,  and  DPAL  small  signal  gain  lineshapes  with  buffer 
gas  pressures  of  1  atm,  10  atm,  and  20  atm  for  each  DPAL  variant  are  provided  in 
Figs.  33,  34,  and  35.  A  similar  spectra  for  a  COIL  device  is  provided  in  Fig.  36. 
When  calculating  molecular  absorption,  HELEEOS  breaks  up  the  slant  path  into 
many  segments,  typically  1,000,  where  section  averages  for  temperature,  pressure,  and 
atmospheric  concentration  are  passed  to  LBLRTM,  which  in  turn  calculates  spectral 
molecular  absorbance.  Figs.  33,  34,  35,  and  36  each  represent  a  single  section  from 
the  1,000  segments  using  the  “ground-to-air”  scenario  described  in  Table  1.  The  laser 
lineshapes  are  used  in  conjunction  with  the  LBLRTM  (km)-1  spectral  absorbance  for 
each  slant  path  segment.  HELEEOS  uses  volume  absorption  coefficient,  a  =  Ecpn;, 
at  each  altitude,  z,  along  the  slant  path,  R. 

The  K  DPAL  emission  line  is  hemmed  in  by  the  molecular  oxygen  A"3E-  to  blTip 
(u”=0,  v'=0)  pP{ 33)  and  PQ{ 33)  lines  on  the  left  and  the  pP(31)  and  PQ(31)  lines 
on  the  right  in  Fig.  33.  The  four  closest  lines  are  assigned  with  absorbances  of  0.04 
or  greater. 

For  a  narrow  line  K  DPAL,  most  of  the  molecular  absorption  occurs  from  the  hot 
bands  (v”—l,  u’=l)  of  molecular  oxygen,  each  nearly  a  magnitude  less  in  absorbance 
than  the  (n”=0,  n’=0)  lines. 
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Figure  33.  Structure  of  atmospheric  absorption  by  molecular  oxygen 
near  the  potassium  (K)  DPAL  D\  line  near  770.1  nm.  Also  shown  is  the 
small  signal  gain  lineshape  with  buffer  gas  pressures  of  1,  10  and  20  atm. 

The  Rb  DPAL  is  very  different  from  both  the  K  and  Cs  variant  -  different  from 
the  K  because  of  water  vapor  instead  of  molecular  oxygen  -  and  different  from  the 
Cs  because  the  Rb  has  a  water  vapor  line  nearly  on  resonance  of  its  laser  emission. 
All  atmospheric  absorption  lines  seen  in  the  spectrum  are  a  result  of  water  vapor, 
and  are  assigned  using  the  vibration  quantum  numbers  of  the  upper  level  (iq’,  v2\ 
t>3 ’)  with  the  ground  vibrational  (0,0,0)  level  common  for  all  the  features.  Rotational 
states  are  labeled  with  the  quantum  numbers  [J\  A' A,  K'c\-[J'\  A'” a,  K” c\  f°r  the 
upper  and  lower  states,  respectively.  The  water  vapor  line  on  resonance  of  the  Rb 
DPAL  is  the  (0,1,3)  [2,1,1]-[2,1,2]  seen  in  Fig.  34. 

CO2  transitions  are  defined  in  HITRAN  as  a  linear  triatomic  with  large  Fermi 
resonance  using  the  definition,  [iq’,  v2\  l2\  V3,r’]-  [iq”,  v2\  l2n ,  v 3” yr” ]  ( Branch ,  J”, 
Symmetry ” ,  A”),  where  v  is  the  quantum  number  associated  with  the  normal  mode  of 
vibration  j,  and  lj  is  the  vibrational  angular  momentum  quantum  number  associated 
with  the  degenerate  bending  mode  j,  and  r  is  the  level  in  the  Fermi  resonating 
group.  A  P-branch  of  CO2  runs  through  the  Rb  emission  line  and  the  CO2  [1,1, 1,5, 2]- 
[0,1, 1,0,1]  (P,l,4,f)  and  (P,l,3,e)  are  on  each  side  of  the  emission  line,  which  also  has 
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Figure  34.  Structure  of  atmospheric  absorption  by  water  vapor  near  the 
rubidium  (Rb)  DPAL  D\  line  near  795.0  nm.  Also  shown  is  the  small 
signal  gain  lineshape  with  buffer  gas  pressures  of  1,  10  and  20  atm. 

some  OH  lines  nearby.  OH,  a  diatomic  with  a  doublet-7r  electronic  state,  is  defined 
in  HITRAN  as  [X\  i\  X” ,  i” ,  Vi”]  ( Branch ,  J”,  Symmetry" ,  F”),  X  being 
the  electronic  state,  and  specifically  for  OH,  the  branch  definition  for  the  lower-state 
quanta  also  accommodates  the  total  orbital  angular  momentum  N  as  well  as  J.  The 
OH  lines  [X, 3/2, 4]-[X, 1/2,0]  (OP,7.5,f,f)  and  (PP,11.5,f,f)  also  lie  on  each  side  of  the 
emission  line.  The  line  strengths  associated  with  these  CO2  and  OH  transitions  along 
with  their  relatively  low  natural  concentrations  result  in  absorption  lines  that  are 
very  small  compared  to  the  water  vapor  lines  nearby. 

The  atmospheric  structure  around  the  Cs  is  very  similar  to  Rb,  dominated  by 
water  vapor  with  only  a  few  OH  and  CO2  lines  of  low  line  strength  and  low  relative 
concentration  to  be  seen  in  Fig.  35. 

It  is  also  clearly  seen  that  the  water  vapor  absorbance  near  the  Cs  DPAL  is  higher 
for  the  same  relative  path  conditions,  but  the  Cs  water  vapor  lines  are  on  each  side, 
surrounded  in  even  more  closely  than  the  K  DPAL  is  by  molecular  oxygen. 

The  COIL  line  seen  in  Fig.  36  demonstrates  the  narrowness  of  the  emission  line, 
a  result  of  cooler  temperatures  and  significantly  lower  operating  gain  cell  pressures 
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Figure  35.  Structure  of  atmospheric  absorption  by  water  vapor  near  the 
cesium  (Cs)  DPAL  D\  line  near  894.6  nm.  Also  shown  is  the  small  signal 
gain  lineshape  with  buffer  gas  pressures  of  1,  10,  and  20  atm. 

because  of  the  fundamental  difference  in  laser  design  between  a  DPAL  and  COIL. 
Again,  only  water  vapor  lines  are  seen  in  the  spectrum,  but  lines  do  exist  in  the 
neighborhood  of  the  COIL  emission  line  from  NO,  CO2,  and  CH4.  The  lines  do 
not  have  the  abundance  or  line  strength  to  contribute  to  absorption,  and  are  also 
relatively  far  away  from  the  narrow  emission  line  with  the  exception  of  the  CO2 
[4, 1, 1, 1,3]- [0,1 ,1 ,0,1]  (R,3,3,e)  that  sits  nearly  at  the  center  of  COIL  emission. 

Five  scenarios  were  designed  to  represent  typical  laser  engagements  shown  in  Fig. 
37.  The  Airborne  Laser  (ABL)  was  megawatt  class  chemical  oxygen-iodine  laser 
(COIL)  on  a  Boeing  747-400F  platform  designed  for  defense  against  intercontinental 
ballistic  missile  (ICBM)  attacks  by  destroying  the  missile  in  boost  phase.  An  “air-to- 
air”  scenario  to  represent  a  typical  ABL  engagement  has  a  high  overall  altitude,  very 
long  slant  path,  and  low  water  content.  The  Advanced  Tactical  Laser  (ATL)  was  a 
COIL  device  installed  on  board  a  C-130  to  support  special  operations.  The  “air-to- 
ground”  scenario  designed  to  represent  it  has  a  moderate  slant  range  propagating  from 
the  aircraft  at  lower  altitude  to  a  ground  level  target.  The  Tactical  High  Energy  Laser 
(THEL)  was  a  ground-based  laser  designed  through  a  joint  effort  between  the  United 
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Figure  36.  Structure  of  atmospheric  absorption  by  water  vapor  near 
the  COIL  iodine  line  near  1.315  /an  is  shown  with  a  typical  emission 
linewidth,  significantly  narrower  than  DPAL  emission  lines. 


States  and  Israel  to  defend  ground  battle  spaces  against  incoming  artillery.  The 
“ground-to-air”  scenario  depicting  its  use  has  a  slant  range  propagating  from  ground 
level  to  an  airborne  target.  The  Navy  has  been  interested  in  using  lasers  for  ship 
defense,  and  while  the  higher  powers  for  ship  defense  have  focused  on  Free  Electron 
Lasers  (FEL),  there  has  also  been  interest  in  smaller  high  energy  laser  to  compliment 
the  Phalanx  close-in  weapon  system  (CIWS).  A  “low  altitude”  scenario  having  a  short 
path  length  with  very  high  overall  water  vapor  content  was  designed  to  represent 
this  case.  Large  Aircraft  Infrared  Countermeasures  (LAIRCM)  systems  such  as  the 
AN/AAQ-24(V)  are  also  considered  with  a  low-power  laser  “jamming”  scenario  with 
low  aircraft  lasing  a  low-altitude  incoming  missile.  Multiple  scenarios  share  many  of 
the  same  input  parameters  and  settings  to  simplify  the  analysis.  Only  the  “air-to-air” 
scenario  uses  adaptive  optics.  All  scenarios  are  set  to  use  the  nearest  Extreme  and 
Percentile  Environmental  Reference  Tables  (ExPERT)  location.  [71,  77]  The  ExPERT 
database  is  a  compilation  of  data  from  573  ground  sites  that  have  collected  specific 
atmospheric  information  for  each  location.  The  final  implementation  of  the  database 
contains  information  that  describes  the  localized  and  upper  air  weather.  Generally, 
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the  database  was  designed  to  describe  the  weather  extremes  experienced  by  each 
ground  site,  thus  information  is  stored  as  a  function  of  percentile  humidity,  time  of 
day,  etc.  [77]  Critical  differences  from  scenario  to  scenario  are  atmospheric  constituent 
content,  overall  engagement  geometry,  laser  wavelength,  laser  power,  laser  spectral 
lineshape.  Table  1  provides  the  primary  differences  for  each  scenario,  including  the 
ratio  of  laser  power  (LP)  to  target  damage  threshold  (TDT),  if  the  aero-optic  model 
includes  a  conformal  aperture,  and  other  general  scenario  settings.  A  visual  overview 
is  included  in  Fig.  37. 


Table  1.  HEL  Engagement  Scenarios 


Scenario  Name 

“air-to-air” 

“air-to-ground” 

“ground-to-air” 

“low  altitude” 

“jamming” 

Start  Alt  (km) 

12 

1.5 

0.01 

0.02 

1 

End  Alt  (km) 

15 

0.001 

2.5 

0.005 

1 

Slant  Range  (km) 

100 

10 

5 

5 

5 

Aerosol  Model 

Continental  Avg. 

Urban 

Desert 

ANAM 

Desert 

Turbulence  Model 

HV  5/7 

HV  5/7 

Climatological 

NSLOT 

HV  5/7 

Approx.  Lat 

39.03 

29.33 

31.87 

22.95 

34.55 

Approx.  Lon 

125.75 

47.52 

35.22 

120.2 

69.22 

ExPERT  Climate  Type 

Maritime 

Desert 

Desert 

Maritime 

Desert 

Aero-Optic  Model 

Conf.  Ap. 

Conf.  Ap. 

None 

None 

Conf.  Ap. 

Exit  Beam  Dia.  (m) 

1.5 

0.5 

0.7 

0.5 

0.02 

Circular  Target  Diameter  (m) 

0.5 

0.25 

0.25 

0.25 

0.05 

LP/TDT  (cm2  •  s-1) 

200 

10 

50 

50 

1 

TDLAS  Device 

The  tunable  diode  laser  absorption  spectroscopy  instrument  has  been  used  to 
validate  simulations  of  atmospheric  paths  and  detailed  descriptions  of  the  device 
configuration  and  methods  for  collecting  spectra  have  previously  been  discussed.  [17, 
65,  64]  The  tunable  diode  lasers  compatible  are  the  New  Focus  Velocity  laser  series, 
here  models  6312,  6314,  and  6318  represent  the  K,  Rb,  and  Cs,  respectively.  Model 
6324  can  be  used  to  represent  a  COIL  laser.  Typical  output  powers  are  about  10 
mW,  and  each  has  a  10  -  20  nm  tunable  range  depending  on  the  model,  all  with  less 
than  a  300  kHz  linewidth.  The  laser  source  can  easily  be  changed  to  investigate  other 
spectral  regions  of  interest  with  minimal  changes  to  the  overall  system.  The  laser 
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Figure  37.  Engagement  geometry  for  a)  “air-to-air”  b)  “air-to-ground” 
c)  “ground-to-air”  d)  “low  altitude”  e)  “jamming”  scenarios. 

output  is  amplitude  modulated  by  a  Thor  Labs  E0-AM-NR-C1  free  space  modulator 
at  102  kHz  (lock-in  maximum  useable  modulation  rate),  and  the  laser  is  fiber  coupled 
and  expanded  to  nearly  fill  a  military  grade  12.5  inch  RC  Optical  RitcheyChretien 
transmit  telescope.  The  laser  can  also  be  split  before  the  fiber  in  order  to  include 
a  reference  cell,  White  cell,  etc.  for  more  detailed  investigations.  The  expanded 
laser  beam  is  directed  across  an  open-path  and  then  received  by  a  second  identical 
telescope.  Thor  Labs  PDA100A  Si  detectors  are  used  to  measure  the  initial  and 
attenuated  signals,  and  are  filtered  by  Stanford  Research  Systems  SR850  dual-phase 
lock-ins.  The  dual-phase  lock-in  signal  is  recorded  by  a  National  Instruments  USB- 
6251  DAQ.  2  m  thick  Thor  Labs  pellicle  beam  splitters  are  used  in  the  optical  chain 
to  minimize  etalon  effects.  The  intensity  of  the  reference  signal  is  measured  late  in 
the  optical  train  immediately  prior  to  the  telescope  beam  expansion.  A  High  Finesse 
WSU-2  wavemeter  records  the  frequency  axis.  The  wavemeter  is  calibrated  with  a 
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2  MHz  fiber-coupled  SIOS  SL-03  frequency-stabilized  HeNe  laser  to  achieve  a  better 
than  10  MHz  accuracy.  All  hardware  in  the  TDLAS  device  is  controlled  by  a  PC 
using  MATLAB®.  Calibration  of  the  frequency  axis,  baseline  removal,  and  further 
spectral  processing  have  also  been  reported.  [17,  65,  64]  Baseline  curvatures  have  a 
typical  0.2  magnitude  in  absorbance  and  are  removed  with  quadratic  fits.  The  baseline 
removal  process  obfuscates  any  wavelength  dependent  scattering,  which  simulations 
suggest  change  by  0.003  in  absorbance  over  the  tunable  range  of  these  lasers  for  a  1 
km  boundary  layer  path.  In  this  configuration,  the  maximum  open-path  distance  is 
limited  by  the  Quickset  QPT-130  ruggedized  pan  and  tilt  mounts  used  and  is  about 
2  km. 

Collections  were  performed  over  a  short  150  m  open-path  during  fair  weather 
conditions.  The  send  and  receive  telescopes  were  placed  1  meter  above  the  ground 
and  separated  by  1  meter,  each  viewing  a  6  inch  flat  mirror  at  the  far  end  of  the  path 
being  investigated.  All  collections  were  performed  at  night  reducing  the  effects  of 
turbulence  during  collection.  Weather  data  was  collected  by  a  NIST-certified  Davis 
Vantage  Pro2  weather  sensor  placed  near  the  turning  mirror. 

Results 

A  best  estimate  for  the  atmospheric  temperature,  pressure,  and  molecular  oxygen 
(K)  or  water  vapor  (Rb  and  Cs)  concentration  were  obtained  using  the  Line-by-Line 
Radiative  Transfer  Model  (LBLRTM)  [19] ,  which  in  turn  uses  the  HITR AN  database 
[72],  By  knowing  the  open-path  distance,  a  non-linear  least  squares  fitting  process  was 
used  by  simulating  an  atmospheric  spectra  defined  for  a  homogenous  atmosphere  with 
a  single  temperature,  pressure,  and  water  vapor  concentration  as  fit  parameters.  This 
process  has  been  described  in  more  detail.  [17]  For  each  collection,  a  low  pressure  (no 
buffer  gas)  alkali  vapor  cell  inside  a  heater  block  was  placed  in  the  path  between  the 
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receive  telescope  and  signal  detector  to  provide  a  reference  for  the  DPAL  wavelength 
being  investigated.  The  spectra  are  illustrated  in  Figs.  38  -  40.  The  strong,  narrow 
absorption  features  near  from  the  D\  2P\/2  -  2*S'i/2  transition  in  potassium,  rubidium, 
and  cesium  are  clearly  seen  near  12,985.19  cm-1,  12,578.95  cm-1,  11,178.27  cm-1, 
respectively,  and  are  excluded  from  the  NLLS  fitting  model.  The  hyperhne  structure 
of  the  Di  transition  is  easily  resolved  at  low  pressures  for  rubidium  and  cesium,  but 
the  potassium  hyperhne  lines  are  spaced  more  closely  than  their  Doppler  widths.  For 
each  spectrum,  a  cubic  baseline  was  removed  before  the  fitting  process  using  only  data 
points  that  lie  outside  of  the  hyperhne  lines  and  the  molecular  oxygen  or  water  vapor 
lines  present  in  the  spectrum.  The  iterative  fitting  process  includes  a  95%  confidence 
interval  for  each  htted  parameter  (temperature,  pressure,  and  concentration)  based 
on  results  from  each  of  the  iterative  hts,  giving  the  final  estimates  for  the  average 
path  temperature,  pressure,  and  constituent  concentration. 

Fig.  38a  shows  the  ht  quality  for  a  150  m  path  with  residuals  in  Fig.  38b  in  the 
neighborhood  of  potassium,  resulting  in  an  estimated  temperature  of  15.5  C  (272.5 
±  3  K),  a  pressure  of  829.1  ±  1  hPa,  and  a  molecular  oxygen  concentration  of  7.011 
±  0.05  x  1018  molecules  per  cm3.  During  the  experiment,  the  weather  station  data 
indicates  variations  during  the  test  of  3.5  K  in  temperature,  10  percentage  points  in 
relative  humidity,  and  0.6  hPa  in  station  pressure.  The  agreement  between  the  ht 
and  meteorology  data  is  poorer  than  the  statistical  error  from  the  spectral  fitting. 
Considering  the  variation  of  weather  parameters  during  the  experiment  combined 
with  the  relatively  low  signal-to-noise  ratio  over  the  higher  rotational  lines,  the  large 
confidence  bounds  from  ht  parameters  is  not  surprising.  Fitting  to  the  entire  P-branch 
would  yield  more  accurate  results,  as  would  including  the  R-branch,  but  would  require 
significantly  longer  collection  times.  [19]  During  the  collection,  an  object  passed 
through  the  open  path  collection  and  is  seen  near  12,995.67  cm-1  and  should  not  be 
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Figure  38.  TDLAS  data  collected  near  the  K  DPAL  emission  line  a)  over 
a  150  m  path  with  NLLS  fit  and  b)  residuals. 

attributed  to  a  spectroscopic  function. 

Fig.  39a  similarly  demonstrates  the  NLLS  fit  for  the  150  meter  path  with  residuals 
in  Fig.  39b  near  the  rubidium  emission  line  resulting  in  an  estimated  temperature 
of  27.8  C  (301.0  ±  10.8  K),  a  pressure  of  815.5  ±  33.4  hPa,  and  a  water  vapor 
concentration  of  4.2141  ±  0.013  x  1 0 1 '  molecules  per  cm3.  This  150  m  outdoor  path 
was  collected  at  night  using  two  weather  stations  at  the  east  side  of  the  path,  and 
using  the  temporal  average  of  data  from  both  instruments  for  the  time  of  collection 
the  temperature  was  21.5  C  (294.65  K),  78.3  %  relative  humidity,  and  985.6  hPa.  The 
agreement  between  the  fit  and  meteorology  data  is  poorer  than  the  statistical  error 
from  the  spectral  fitting,  just  as  with  the  potassium  case.  However,  absorbance  of 
the  water  lines  nearby  rubidium  are  nearly  an  order  of  magnitude  smaller  than  those 
near  the  cesium  DPAL  emission  lines,  as  can  be  seen  in  the  fit.  As  expected  fitting  to 
smaller  absorption  lines  result  in  poorer  fitting  statistics,  and  demonstrate  the  lower 
overall  molecular  absorbance  near  the  Rb  DPAL. 

Fig.  40a  illustrates  the  fit  quality  for  the  150  meter  path  with  residuals  in  Fig. 
40b  resulting  in  an  estimated  temperature  of  15.5  C  (288.6  ±  0.15  K),  a  pressure  of 
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Figure  39.  TDLAS  data  collected  near  the  Rb  DPAL  emission  line  a) 
over  a  150  m  path  with  NLLS  fit  and  b)  residuals. 


993.8  ±  0.9  hPa,  and  a  water  vapor  concentration  of  3.096  ±  0.004  x  101'  molecules 
per  cm3.  Weather  conditions  during  the  experiment  collection  were  fairly  constant, 
varying  less  than  1  degree  K,  less  than  4  percentage  points  in  relative  humidity,  and 
less  than  0.4  hPa  in  pressure.  The  agreement  between  the  measurement  and  the 
meteorology  data  is  poorer  than  the  statistical  error  from  the  spectral  fitting  and 
but  fitted  temperature  still  agrees  with  meteorological  temperature  to  within  2.1%, 
excellent  agreement. 

A  typical  simulated  absorption  spectrum  for  a  1  km  atmospheric  path  used  to 
calculate  integrated  absorption  coefficient  for  a  Cs  DPAL  is  shown  in  Fig.  41a  along 
with  the  product  of  the  small  signal  gain  and  absorption  coefficient  before  integration 
in  Fig.  41b. 

Using  Eq.  (17)  the  final  integrated  absorption  of  each  laser  type  for  each  scenario 
is  summarized  in  Fig.  42. 

The  “air-to-air”  scenario  of  in  Fig.  42a  demonstrates  the  disadvantage  a  K  DPAL 
for  over  extreme  path  lengths.  Only  by  making  a  K  DPAL  with  a  narrow  emission 
line  results  in  molecular  absorbance  that  rivals  the  Rb  variant  with  similar  gain  cell 
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pressure.  Each  flavor  of  DPAL  performs  better  than  COIL  in  regards  to  molecular 
absorption,  with  the  exception  of  a  20  atm  K  and  Cs  getting  close  to  the  COIL 
molecular  absorption  value,  K  because  of  its  long  path  length  and  Cs  because  of  the 
stronger  water  lines  nearby.  Surprisingly,  the  Rb  DPAL  benefits  from  the  higher  buffer 
gas  pressures  because  the  water  vapor  absorption  line  that  lies  nearly  on  top  of  the 
Rb  D\  line  -  while  the  K  and  Cs  DPALs  are  hemmed  in  by  their  respective  absorption 
features.  The  Rb  escapes  absorption  by  spreading  its  power  over  a  broader  spectral 
area.  Fig.  42b  shows  the  results  from  a  “air-to-ground”  scenario,  demonstrating  how 
a  shorter  path  begins  to  favor  the  K  variant,  putting  its  performance  on  par  with  the 
Rb  because  of  lower  molecular  oxygen  in  the  overall  path.  For  this  case,  both  the  Cs 
and  COIL  lasers  have  about  the  same  molecular  absorption  for  higher  Cs  gain  cell 
pressures.  Fig.  42c  shows  how  a  “ground-to-air”  scenario  also  gives  similar  results 
when  comparing  a  highly  pressure-broadened  Cs  DPAL  to  a  COIL  device  while  the 
K  and  Rb  DPAL  performances  are  similar.  Given  that  its  path  length  is  twice  as 
long,  it  is  no  surprise  the  molecular  absorption  is  approximately  twice  that  of  the 
“air-to-ground”.  The  over-water  “low  altitude”  scenario  in  Fig.  42d  demonstrates 
the  advantage  a  K  DPAL  would  have  over  these  shorter  paths  with  very  high  water 
vapor  while  the  Rb  and  Cs  DPALs  and  COIL  are  all  affected  by  the  water  vapor 
content  near  the  sea  surface.  Finally,  Fig.  42e  shows  similar  trend  for  the  scenario 
describing  laser  jamming  of  IR  seeking  missiles.  The  worst  case  K  and  Rb  DPAL 
configurations  perform  similarly  to  the  best  case  Cs  configuration.  A  worst  case  Cs 
experiences  similar  molecular  absorption  as  the  COIL,  a  trend  shared  in  all  plots. 

Molecular  absorption  is  not  the  only  metric  of  interest,  and  all  inputs  and  outputs 
from  HELEEOS  are  available  to  compare  overall  laser  performance.  Scattering  at 
shorter  wavelengths  starting  in  the  visible  regime  and  shorter  by  air  molecules  can 
be  significant.  The  scattering  cross-section  of  air  molecules  is  approximately  pro- 
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portional  to  the  Rayleigh  regime,  of  which  a  COIL  device  lies  at  the  border  of,  or 
arguably  beyond.  While  scattering  reduces  fluence  on  target,  it  also  reduces  thermal 
blooming,  a  trade  off  which  is  computed  in  HELEEOS  simulations.  Thus,  extinction 
due  to  scattering  can  be  nearly  8  times  greater  for  a  K  DPAL  compared  to  a  COIL 
device.  For  example,  the  “air-to-air”  scenario  calculated  total  aerosol  and  molecular 
scattering  cross  section  to  be  0.47,  0.41,  0.26,  and  0.06  for  K,  Rb,  Cs,  and  COIL 
respectively.  Losses  clue  to  diffraction  have  the  opposite  effect,  with  the  diffraction- 
limited  beam  diameter  on  target  approximately  proportional  to  wavelength,  giving 
a  smaller  minimum  focusable  spot  size  for  a  DPAL  laser  and  increasing  irradiance. 
For  example,  the  “air-to-ground”  scenario  has  a  diffraction  limited  beam  width  of 
0.019,  0.021,  0.024,  and  0.034  for  K,  Rb,  Cs,  and  COIL  respectively.  Strchl  ratio,  a 
numerical  value  ranging  from  zero  to  one,  represents  many  effects  that  increase  the 
effective  spot  size,  and  reduces  the  peak  at  beam  center  to  a  lower  value  than  the 
peak  for  the  diffraction  limited  case.  Strehl  ratio,  S,  is  defined  here  as  the  ratio  of 
transmitted  to  diffraction-limited  irradiance,  /,  at  beam  center, 


I(x  =  0  ,y  =  0) 
ID(x  =  0,2/  =  0) 


(19) 


If  the  Strchl  ratio  is  unity,  the  system  performance  is  considered  “diffraction  limited” . 
As  beams  are  spread  by  atmospheric  effects  like  turbulence  and  thermal  blooming, 
the  Strchl  ratio  is  lowered  as  the  irradiance  of  the  beam  is  spread  outward.  Fig.  43 
demonstrates  Strehl  ratio  for  the  various  DPAL  variants  and  scenarios  Fig.  42. 

Just  as  with  molecular  absorption,  performance  of  Strehl  ratios  for  Rb  DPAL 
gain  cell  buffer  gas  pressures  are  typically  inverted  compared  to  K  and  Cs,  except  for 
the  low  power  “jamming”  scenario  where  all  DPAL  cases  are  very  close  to  identical. 
Overall,  it  is  also  seen  that  molecular  absorption  has  an  effect  on  overall  strehls,  with 
a  typical  Strehl  ratio  difference  from  the  minimum  pressure  to  the  maximum  pressure 
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of  about  15%,  12%,  and  22%  for  K,  Rb,  and  Cs  respectively,  further  demonstrating 
how  the  Cs  DPAL  is  more  dependent  on  gain  cell  pressure  than  K  or  Rb.  Interestingly, 
beam  spread  is  greater  for  all  DPALs  than  for  COIL,  giving  COIL  the  highest  overall 
Strehl  ratio  in  nearly  every  case  with  the  exception  of  the  water-rich  “low  altitude” 
case  over  the  water.  Separating  the  results  for  overall  Strehl  ratio  for  each  individual 
run  into  components  cause  is  difficult  because  of  the  internal  design  of  HELEEOS 
and  its  dependence  on  the  High  energy  laser  And  Relay  Engagements  (SHaRE),  an 
internal  toolkit  provided  by  MZA  Associates. 

Dwell  time,  D^,  was  calculated  using  a  90%  confidence  in  the  probability  of  kill 
and  using  target  damage  thresholds  20%  of  the  overall  target  damage  threshold.  As 
expected,  dwell  time  is  inversely  related  to  beam  Strehl  ratio,  with  very  high  Strehl 
ratios  exhibiting  the  lowest  dwell  times.  The  “air-to-air”  scenario  resulted  in  dwell 
times  with  little  difference  resulting  from  molecular  absorption  for  different  DPAL 
gain  cell  pressures,  and  with  COIL  performing  best  over  very  long  paths  from  its 
superior  molecular  scattering  performance,  and  coupled  with  low  absorption  from 
high  altitudes  and  low  pressures,  there  is  little  thermal  blooming  in  this  scenario.  The 
“air-to-ground”  and  “ground-to-air”  had  similar  trends,  with  K  and  Rb  performing 
similarly  for  their  best  case  pressures,  and  COIL  having  the  longest  required  dwell 
time.  The  “low  altitude”  scenario  demonstrated  a  high  dependence  on  DPAL  gain 
cell  pressures  for  Rb  and  Cs,  with  K  easily  outperforming  all  other  devices,  a  likely 
result  given  the  high  water  vapor  content  over  water. 

For  reference,  a  plot  of  diffraction- limited  spot  area  compared  to  scattering  cross- 
section  and  turbulence  relative  variance,  each  normalized  to  one,  is  shown  in  Fig. 
45.  Higher  relative  values  indicate  poorer  beam  quality  on  target.  The  area  of  the 
diffraction- limited  beam  size  is  proportional  to  A2,  while  scattering  cross-section  is 
proportional  to  A-4,  and  turbulence  is  proportional  to  A~7/6,  with  scattering  having 
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Figure  44.  Dwell  time  for  a)  “air-to-air”  b)  “air-to-ground”  c)  “ground- 
to-air”  d)  “low  altitude”  e)  “jamming”  scenarios  for  DPAL  lasers  with 
gain  cells  containing  buffer  gasses  at  1  atm  (slant  right),  10  atm  (dotted), 
and  20  atm  (slant  left)  of  pressure,  all  compared  to  a  COIL  laser  (white). 
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Figure  45.  Diffraction  Limited  Spot  Area  (-),  scattering  cross-section  (:), 
and  turbulence  relative  variance  (-  -),  each  normalized  to  one  to  showing 
each  effect’s  spectral  dependence. 

the  largest  relative  change  in  overall  effect  over  the  wavelength  range  because  of  its 
A-4  dependence.  Because  of  scattering’s  strong  degradation  of  shorter  wavelength 
beams,  laser  devices  with  shorter  wavelengths  will  be  greatly  affected  by  scattering 
over  long  slant  ranges,  a  trend  seen  in  Fig.  44a. 

With  Rb,  Cs,  and  COIL  devices  being  dependent  on  water  vapor  concentration, 
the  distribution  of  relative  humidity  recordings  at  each  site  using  data  from  the  nearest 
ExPERT  site  was  used  to  set  relative  humidity  for  a  one-percentile  day,  a  fiftieth- 
percentile  day,  and  for  a  ninety-ninth-percentile  day,  where  relative  humidity  on  a 
one-percentile  day  would  be  very  low,  for  a  fiftieth-percentile  day  would  have  median 
humidity  found  at  the  site  used,  and  for  a  ninety-ninth-percentile  day  would  exhibit 
very  high  relative  humidity.  The  same  series  of  simulations  for  the  five  scenarios  dis¬ 
cussed  above  were  carried  out  with  a  10  atm  gain  cell  for  each  DPAL  and  the  typical 
COIL  low  pressure  gain  cell.  As  expected,  molecular  absorption  for  the  “air-to-air” 
case  is  nearly  identical  for  each  percentile  relative  humidity  due  to  the  low  water 
vapor  content  and  minimal  broadening  of  those  lines  at  high  altitudes.  A  close  look 
at  the  “ground-to-air”  case  reveals  the  dependence  Rb,  Cs,  and  COIL  devices  have 
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(b)  (c) 

Figure  46.  Variation  in  percentile  relative  humidity  for  a)  integrated 
absorbance,  b)  strehl  ratio,  and  c)  dwell  time  for  the  “ground-to-air” 
scenarios  for  DPAL  a  laser  with  a  10  atm  gain  cell  buffer  gas  pressure 
each  with  a  1%  (horizontal  lines),  a  50%  (slanted  left-slanted  right),  and 
a  99%  (horizontal-vertical)  percentile  relative  humidity  days. 

on  water  vapor  content.  Fig.  46a  demonstrated  the  vast  differences  in  molecular 
absorbance  from  a  one-percentile  to  a  ninety-ninth-percentile  relative  humidity  day, 
where  all  water  dependent  devices  show  the  change  in  integrated  molecular  absorp¬ 
tion  as  relative  humidity  increases.  For  example,  the  COIL  device  has  over  5  times 
higher  molecular  absorbance  for  the  “ground-to-air”  between  the  one-percentile  to 
a  ninety-ninth-percentile  relative  humidity  day.  This  higher  overall  molecular  ab¬ 
sorbance  contributes  to  overall  beam  Strehl  ratio,  seen  in  Fig.  46b,  where  the  COIL 
still  has  the  best  Strehl  ratio  performance.  Overall  performance  for  dwell  times  in 
Fig.  46c  demonstrates  how  very  high  water  content  adds  nearly  50%  more  dwell  time 
to  heavily  water  dependent  devices.  As  expected,  the  K  and  Rb  DPAL  devices  are 
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not  as  susceptible  to  seasonal  and  daily  changes  in  relative  humidity,  and  can  been 
seen  in  this  example  in  variation  of  relative  humidity  for  the  “ground-to-air”  scenario. 

Conclusions 

The  TDLAS  system  is  shown  to  collect  information  for  molecular  absorption  over 
atmospheric  paths  accurately  estimating  atmospheric  content  validating  atmospheric 
models.  K  DPAL  atmospheric  transmission  will  be  excellent  over  atmospheres  with 
high  water  vapor  and  shorter  path  lengths  with  lower  laser  linewidths.  Rb  DPAL 
atmospheric  transmission  will  be  excellent  over  atmospheres  with  lower  water  vapor 
concentrations  and  will  perform  relatively  better  than  the  Cs  DPAL  variant.  Rb 
DPAL  lasers  will  benefit  from  broader  laser  linewidths  by  spreading  laser  emission 
away  from  the  water  line  located  near  resonance.  Cs  DPAL  atmospheric  transmission 
will  be  the  poorest  compared  to  Rb  and  K  variants,  and  performs  best  with  very  little 
water  vapor.  Cs  DPALs  will  benefit  from  lower  pressure  gain  cell  configurations  over 
hydrated  atmospheric  paths.  All  DPAL  systems  suffer  from  higher  molecular  scatter¬ 
ing  compared  to  COIL  devices,  and  DPAL  Strehl  ratios,  and  their  associated  dwell 
times  are  less  desirable  than  those  from  equivalent  power  COIL  devices.  Seasonal 
distribution  of  relative  humidity  can  result  in  dwell  times  50%  longer  for  Cs  DPAL 
and  COIL  devices  while  having  little  effect  on  K  and  Rb  DPALs  and  little  effect  on 
all  systems  at  high  altitudes.  Performance  will  depend  strongly  on  details  of  the  laser 
lineshape. 


91 


VII.  Conclusion 


A  Tunable  Diode  Laser  Absorption  Spectroscopy  (TDLAS)  apparatus  has  been 
designed,  built,  characterized,  and  deployed  for  High  Energy  Laser  (HEL)  weapons 
applications.  Long,  open-path  atmospheric  characterization  with  path  lengths  ex¬ 
ceeding  1  km  are  studied  and  documented.  The  prior  state-of-the  art  limited  path 
lengths  to  150  m.  The  instrument  has  been  ruggedized  and  deployed  in  the  held, 
operating  at  temperatures  of  -9  C  to  25  C.  A  total  17  indoor  and  outdoor  campaigns 
were  executed  collecting  a  total  of  43  sets  of  data  totaling  over  500  scans  of  the  potas¬ 
sium,  rubidium,  and  cesium  D\  emission  wavelengths.  The  TDLAS  system  is  being 
prepared  for  deployment  to  DoD  HEL  test  ranges,  including  the  October  2012  test 
scheduled  at  the  Naval  Surface  Warfare  Center  (NSWC)  directed  energy  test  facility 
in  Dahlgren,  VA. 

The  high  spectral  resolution  (300  kHz)  applied  over  the  broad  spectral  coverage 
(over  6  THz)  of  the  instrument,  and  multiple  wavelength  regions  diode  laser  wave¬ 
lengths  allow  for  detection  of  more  than  72  O2  spectral  lines  and  over  100  H2O  lines. 
High  sampling  rates  yield  over  a  million  samples  over  the  full  tunable  diode  range, 
1,500  samples  per  Doppler  limited  lineshape,  and  15,000  per  atmospheric  broadened 
linshape  enabling  robust  Voigt  profile  characterization  and  extraction  of  both  tem¬ 
perature  and  pressure.  Absolute  wavelength  accuracy  of  10  MHz  was  achieved  using 
a  wavemeter  limited  by  a  stabilized  HeNe  calibration  source  and  was  validated  by 
both  alkali  frequency  standards  and  by  pressure  shifted  and  broadened  atmospheric 
lines  from  the  HITRAN  database. 

The  minimum  detectable  peak  absorbance  in  a  laboratory  environment  is  0.01. 
System  jitter  and  atmospheric  turbulence  introduce  additional  noise  for  field  deploy¬ 
ments,  yielding  a  detectable  absorbance  of  0.015  over  a  150  m  paths  and  0.04  over 
1  km  paths.  Most  collections  took  place  at  times  of  the  day  when  turbulence  is 
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low,  pointing  to  telescope  jitter  as  the  dominant  source  of  noise.  Even  for  the  very 
weak  O2  to  &1E+  magnetic  dipole  transitions,  the  signal-to-noise  ratio  under 

field  conditions  exceeds  1,000  with  a  peak  signal  of  2.51  in  absorbance  with  a  0.0025 
noise  standard  deviation.  The  12  inch  diameter  collimated  beam  transmitted  over 
1  km  maintains  excellent  collimation,  and  at  over  4  km  distances  maintains  an  12 
inch  diameter  but  because  of  platform  jitter  and  turbulence  the  overall  beam  center 
moves  up  two  2  inches  from  its  center  over  1  km  paths  making  4  km  paths  impossible 
without  larger  turning  optics. 

The  instrument  is  supported  by  instrument  control  software  based  in  MATLAB®  us¬ 
ing  the  Instrument  Control  Toolbox  and  the  Data  Acquisition  Toolbox,  as  well  as 
manufacturer-supplied  shared  libraries  customized  for  use  with  the  TDLAS  appa¬ 
ratus.  Conversion  of  observed  signals  to  calibrated  absorbance  as  function  of  cal¬ 
ibrated  wavelength  is  also  achieved  by  post-processing  methods  also  written  and 
implemented  in  MATLAB®,  and  spectral  interpretation  to  yield  atmospheric  tem¬ 
peratures,  pressures  and  species  concentrations  were  based  on  spectroscopic  methods 
and  non-linear  least-squares  fitting  methods  using  the  Line-by-Line  Radiative  Trans¬ 
fer  Model  (LBLRTM).  The  new  TDLAS  instrument  represents  a  significant  advance 
of  existing  instrument,  investigating  broad  spectral  areas  at  ultra-high  spectral  reso¬ 
lution  over  long  outdoor  paths. 

The  O2  A f3E“  to  51S+  (0,0)  rotational  spectrum  in  the  vicinity  of  762  nm  has 
been  characterized  by  the  TDLAS  instrument.  Temperature  was  determined  to  within 
0.7%,  concentration  to  0.8%,  and  atmospheric  pressure  to  less  than  4%.  The  TDLAS 
spectra  agree  with  the  meteorological  data  for  temperature  to  within  1.3%,  for  con¬ 
centration  to  within  1.6%,  and  to  atmospheric  pressure  within  less  than  10%.  Hot 
bands  can  be  observed  over  3  km  path  lengths  in  multi-pass  cells.  Broad  spectral  band 
width  allows  for  up  to  72  rotational  features  to  be  observed  in  a  single  unified  scan  and 
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direct  observation  of  Boltzman  statistical  distribution  to  extract  the  J  dependence 
of  pressure  broadening  rates  observed  from  Voigt  lineshapes.  The  results  suggest  a 
possible  improvement  to  HITRAN  database  as  the  held  data  agrees  more  closely  with 
Pope’s  [62]  results.  Validation  of  Monocular  Passive  Ranging  (MPR)  using  Fourier 
Transform  Interferometer  (FTIR)  device  was  achieved,  demonstrating  that  the  FTIR 
has  a  variance  in  measured  absorbance  5  times  greater  than  the  TDLAS  device  mea¬ 
surement.  Spectral  observations  are  more  sensitive  to  the  state  of  atmosphere  than 
most  modest  weather  stations  predict.  Observed  baseline  changes  slowly  over  tuned 
wavelength  but  can  be  up  to  10%  of  the  maximum  signal  observed  and  attributing 
any  spectral  dependence  of  the  scattering  is  problematic.  Total  scattering  extinction 
measurements  may  require  off-axis  methods. 

Atmospheric  absorbance  for  all  three  DPAL  variants  was  characterized  in  the 
held  with  path  lengths  up  to  1  km.  A  high  spectral  resolution  is  required  for  the 
characterization  of  HEL  propagation  due  to  the  narrow  linewidths  of  laser  emission, 
making  techniques  like  FTIR  with  lower  0.25  cm-1  resolution  inadequate.  Transmis¬ 
sion  spectra  were  recorded  with  atmospheric  features  co-registered  with  low  pressure 
alkali  cells,  yielding  wavelength  calibration  to  within  10  MHz,  and  nearly  2  MHz 
in  close  spectral  proximity  to  alkali  features.  Field  data  was  interpreted  to  validate 
LBLRTM  simulations  using  the  HITRAN  database.  With  improved  atmospheric  con¬ 
stituent  estimates,  no  signihcant  variation  between  simulations  and  observations  can 
be  observed  within  the  statistical  limits  of  the  instrument.  However,  improved  mete¬ 
orological  estimates  for  the  state  of  atmosphere  are  required  to  adequately  simulate 
slant  paths,  and  the  need  for  high  fidelity  weather  instrumentation  is  apparent.  The 
02  X3E~  to  b 1 E+  is  critical  to  the  potassium  DPAL,  especially  the  p P  (31,33)  and 
PQ  (31,33)  features  and  nearly  hot  bands.  Water  vapor  is  important  to  both  the 
rubidium  and  cesium  DPALs,  with  cesium  having  water  vapor  lines  nearly  an  order 


94 


of  magnitude  greater  than  those  near  rubidium  DPAL  emission.  Water  vapor  was 
used  to  measure  temperature,  pressure  and  concentration  to  within  2.4%  or  less  from 
150  m  path  collections  and  3.7%  for  1  km  path  collections  when  compared  to  mete¬ 
orological  data.  OH  and  CO2  lines  are  scattered  near  DPAL  emission  lines,  however 
the  atmospheric  abundance  of  these  species  are  so  low  absorption  lines  do  not  im¬ 
pact  DPAL  atmospheric  transmission.  Potassium  DPAL  devices  will  be  insensitive 
to  changes  in  relative  humidity,  and  rubidium  and  cesium  DPALs  will  have  the  least 
molecular  absorption  at  high  altitude  and  dry  slant  paths. 

The  impact  of  atmospheric  transmission  on  High  Energy  Laser  weapon  perfor¬ 
mance  for  the  emerging  DPAL  systems  was  assessed.  LBLRTM  was  integrated  into 
the  High  Energy  Laser  End-to-End  Operational  Simulation  (HELEEOS)  to  compute 
molecular  absorption  and  five  scenarios  were  considered  in  this  broad  application 
study:  a  100  km  slant  path  at  high  altitude  “air-to-air”,  a  10  km  “air-to-ground” 
slant  path,  a  5  km  “ground-to-air”  slant  path,  an  over-ocean  “low  altitude”  slant 
path,  and  a  low  power  “jamming”  scenario.  When  considering  only  molecular  ab¬ 
sorption,  the  potassium  and  rubidium  DPALs  seems  to  be  the  clear  winners,  however 
when  considering  the  system  engineering  metric  strehl  ratio,  the  differences  between 
DPAL  variants  are  less  exaggerated,  with  the  only  clear  scenario-specific  standout  be¬ 
ing  potassium  for  very  well  hydrated  slant  paths.  The  warfighter  metric,  dwell  time, 
also  highlights  the  advantage  of  a  laser  system  that  operates  independently  of  water 
vapor  absorption.  It  also  clearly  demonstrates  the  disadvantage  shorter  wavelengths 
have  in  scattering  compared  to  existing  technologies  like  the  Chemical  Oxygen  Iodine 
Laser  (COIL),  making  a  cesium  DPAL  an  excellent  choice  for  the  dry  and  long  path 
in  a  high  altitude  ” air-to-air”  engagement. 
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Recommendations 


As  DPAL  systems  become  more  mature,  additional  opportunities  will  be  avail¬ 
able.  Using  higher  source  laser  powers  could  be  used  to  study  side  scattering,  and 
with  careful  experiment  geometry  the  angular  dependence  of  scattering  could  also 
be  studied.  Performing  experiments  in  heavy  aerosol  environments  will  also  improve 
chances  of  viewing  and  quantifying  side  scatter  and  could  be  validated  with  particle 
counters  and  a  temporal  variation.  Additional  validation  could  be  achieved  by  col¬ 
lecting  information  about  the  size  distribution  of  scatterers.  Specialized  balanced  and 
calibrated  photo  detectors  could  reduce  baseline  variation  in  a  single  scan.  Variations 
in  noise  baseline  could  also  be  captured  using  shorter  lock-in  time  constants  on  the 
order  of  10  msec.  Collecting  turbulence  data  over  the  same  path  simultaneously  us¬ 
ing  the  TDLAS  system  and  an  off  the  shelf  turbulence  measurement  device  like  the 
BLS900  may  demonstrate  the  capability. 

The  maximum  effective  open-path  length  the  of  the  TDLAS  system  could  be 
extended  by  upgrading  telescope  mounts,  altering  the  optical  system  design,  and  in¬ 
creasing  laser  power.  Signal-to-noise  ratio  could  be  improved  by  increasing  noise  sup¬ 
pression  with  faster  modulation  frequency  and  dynamic  range  could  be  improved  with 
upgraded  lock-ins  operating  at  significantly  higher  modulation  frequencies.  Whole- 
system  reliability  and  ease  of  use  could  be  improved  by  transitioning  all  optical  ele¬ 
ments  in  the  optical  board  to  fiber  coupled  devices  and  enclosing  sensitive  devices  in 
thermally  controlled  enclosures.  Large  hollow  retro-reflective  arrays  could  be  used  to 
eliminate  the  receive  telescope  and  make  the  alignment  process  significantly  easier, 
but  may  exhibit  etelon  effects. 

Improvements  to  customized  hardware  control  could  be  implemented  increasing 
the  probability  of  successful  data  collection.  Closed  absorption  cells  could  be  in¬ 
cluded  during  collections  to  improve  signal-to-noise  for  detailed  studies  of  very  small 
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absorption  features.  Different  modulation  techniques  could  be  used  to  track  single 
absorption  features  over  long  time  periods  to  characterize  temporal  variation  along 
paths.  The  TDLAS  system  could  be  deployed  in  a  variety  of  different  atmospheric 
environments  or  place  in  geometries  similar  to  those  where  a  DPAL  may  be  deployed. 
Stable  and  unstable  DPAL  resonators  could  be  simulated  and  experimentally  tested 
to  determine  possible  emission  lineshapes  which  could  then  be  applied  to  make  at¬ 
mospheric  transmission  estimates. 

The  device  could  also  be  used  to  monitor  species  not  investigated  here,  like  CO2, 
OH,  and  virtually  any  species  within  the  range  of  tunable  diode  lasers  available. 

Comparisons  between  LBLRTM  and  the  Laser  Environmental  Effects  Definition 
and  Reference  (LEEDR),  and  LEEDR  can  be  used  to  upgrade  internal  HELEEOS 
molecular  absorption  calculations.  More  detailed  HELEEOS  studies  can  be  carried 
out,  including  different  scenarios,  comparisons  to  other  modern  HEL  systems,  and 
more  studies  of  the  effects  diurnal  and  seasonal  variations  have  on  DPAL  propaga¬ 
tion. 
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Appendix  A.  Physical  Constants 


Table  2.  Physical  constants^ 


Symbol 

Description 

Value/Units 

Oifs 

Fine-structure  constant 

1/137.035  999  679(94) 

c 

Speed  of  light  in  vacuum 

299  792  458  ms-1 

h 

Planck  constant 

6.626  068  96(33)  x  10"34Js 
4.135  66733(10)  x  10~15eVs 

h 

Planck  constant  /  2n 

1.054  571628(53)  x  10~34Js 
6.582118  99(16)  x  KT16eVs 

ks 

Boltzmann  constant 

1.380  650  4(24)  x  10"23  JK"1 
8.617343(15)  x  lO^eVK”1 

eo 

Permittivity  (Electric  constant) 

8.854187817  x  lO-^Fm"1 

IM) 

Permeability  (Magnetic  constant) 

12.566  370  614  x  10-7NA“2 

e 

Elementary  charge 

1.602176  487(40)  x  10~19C 

me 

Electron  mass 

9.109  38215(45)  x  10"31  kg 

rrip 

Proton  mass 

1.672  621637(83)  x  10~27kg 

mu 

Atomic  mass  constant 

1.660  538  782(83)  x  10~27kg 

Patm 

Standard  atmosphere 

101  325  Pa  =  760  Torr 

t  NIST  2006  CODATA  recommended  values 


Appendix  B.  Alkali  Data 


Potassium  Properties 


Table  3.  Potassium  (K)  physical  properties 


Property 

Symbol 

Value 

Atomic  number' 

Z 

19 

Total  nucleons* 

Z  +  N 

58 

Ground  state'*' 

— 

4  s2Sl/2 

Relative  natural  abundance* 

rj(39K ) 

93.2581(44)% 

Atomic  mass'* 

m 

39.0983  u 

Melting  point*’ 

Tm 

63.38  °C 

Boiling  point* 

Tb 

759  °C 

Density  (solid)* 

P 

856  kg  m-3 

Nuclear  spin* 

I 

3/2 

Magnetic  moment 

P 

+0.391465  8 

Ionization  energy* 

Ej 

4.340  663  45  eV 

*  Steck,  Daniel  A.  “Rubidium  87  D  Line  Data.”  2001 


Table  4.  Potassium  D\  (42S1/2  — ►  42P]y2)  transition  optical  properties [82] 


Property 

Symbol 

Value 

Frequency 

389.286  058  716(62)  THz 

Energy 

his0 

1.609  957  707  081 91  eV 

Wavelength 

A 

770.108  385  049(123)  nm 

Lifetime 

r 

26.37(5)  ns 

Natural  Line  Width  (FWHM) 

6.03(1)  MHz 

Rubidium  Properties 


Cesium  Properties 
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Table  5.  Potassium  Z)2  (42S!/2  — ►  42P3/2)  transition  optical  properties [82] 


Property 

Symbol 

Value 

Frequency 

D) 

391.016  170  03(12)  THz 

Energy 

hv 0 

1.617112  872  240  59  eV 

Wavelength 

A 

766.700  921 822(24)  nm 

Lifetime 

T 

26.37(5)  ns 

Natural  Line  Width  (FWHM) 

6.035(11)  MHz 

Table  6.  Rubidium  (Rb)  physical  properties 


Property 

Symbol 

Value 

Atomic  number' 

Z 

37 

Total  nucleons* 

Z  +  N 

87 

Ground  state! 

— 

5s2Ai/2 

Relative  natural  abundance! 

r](87Rb ) 

27.83% 

Atomic  mass' 

m 

86.909180  520(15)  u 

Melting  point' 

Tm 

39.31  °C 

Boiling  point' 

Tb 

688  °C 

Density  (solid)  1 

p 

1532  kg  m-3 

Nuclear  spin' 

I 

3/2 

Magnetic  moment 

p 

+2.75124 

Ionization  energy! 

Ej 

3.893  90(2)  eV 

!  Tiecke,  T.G.  “Properties  of  Potassium.”  2010 


Table  7.  Rubidium  D\  (52S!/2  — ►  52Pi/2)  transition  optical  properties [78] 


Property 

Symbol 

Value 

Frequency 

D) 

377.107463  5(4)  THz 

Energy 

hv 0 

1.559  590  99(6)  eV 

Wavelength 

A 

794.978  850  9(8)  nm 

Lifetime 

T 

27.70(4)  ns 

Decay  Rate 

r 

36.10(5)  x  HRs”1 

Natural  Line  Width  (FWHM) 

5.746(8)  MHz 

Absorption  Oscillator  Strength 

/ 

0.3420(14) 

100 


Table  8.  Rubidium  D 2  (52S!/2  — »•  52P3/2)  transition  optical  properties[78] 


Property 

Symbol 

Value 

Frequency 

D) 

384.230  484 468  5(62)  THz 

Energy 

hv 0 

1.589  049439(58)  eV 

Wavelength 

A 

780.241  209  686(13)  nm 

Lifetime 

T 

26.24(4)  ns 

Decay  Rate 

r 

38.11(6)  x  lOS"1 

Natural  Line  Width  (FWHM) 

6.065(9)  MHz 

Absorption  Oscillator  Strength 

/ 

0.6956(15) 

Table  9.  Cesium  (Cs)  physical  properties 

Property 

Symbol 

Value 

Atomic  number^ 

Z 

55 

Total  nucleons! 

Z  +  N 

133 

Ground  state! 

— 

6s2£i/2 

Relative  natural  abundance! 

rj(133Cs) 

100% 

Atomic  mass! 

m 

132.905452(1)  u 

Melting  point! 

Tm 

28.44  °C 

Boiling  point! 

Tb 

671.0  °C 

Density  (solid)! 

P 

1879  kgm-3 

Nuclear  spin! 

I 

7/2 

Magnetic  moment 

P 

+2.579 

Ionization  energy! 

Ej 

3.893  90(2)  eV 

^  Steck,  Daniel  A.  “Cesium  D  Line  Data.”  1998 

Table  10.  Cesium  D\  (62S!/2  — »  6" 

*P !/2)  transition  optical  properties[79] 

Property 

Symbol 

Value 

Frequency 

D) 

335.116  048  807(41)  THz 

Energy 

hv o 

1.385  928  475(50)  eV 

Wavelength 

A 

894.592  959  86(11)  nm 

Lifetime 

r 

34.894(44)  ns 

Decay  Rate 

r 

28.659(36)  x  lOS”1 

Natural  Line  Width  (FWHM) 

4.5612(57)  MHz 

Absorption  Oscillator  Strength 

/ 

0.3438(13) 

101 


Table  11.  Cesium  D 2  (62Si/2  — >  62P3/2)  transition  optical  properties [79] 


Property  Symbol  Value 


Frequency 
Energy 
Wavelength 
Lifetime 
Decay  Rate 

Natural  Line  Width  (FWHM) 
Absorption  Oscillator  Strength 


D) 

351.725  718  50(11)  THz 

hv 0 

1.454  620  542(53)  eV 

A 

852.347275  82(27)  11111 

T 

30.473(39)  ns 

r 

32.815(41)  x  lOS"1 

5.2227(66)  MHz 

/ 

0.7148(13) 
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Appendix  C.  Matlab®  Code 


Importing  PeakFit  Fitting  Data 

This  code  was  created  to  take  PEAKFlT®fitting  results  and  import  them  into 
Matlab®.  The  code  is  specifically  designed  to  find  a  specific  type  of  Voigt  fit  in 
PeakFit,  the  Voigt  G/L,  but  can  easily  be  altered  to  extract  other  types  of  lines 
available  in  PeakFit®.  The  algorithm  runs  through  a  standard  PEAKFlT®numeric 
export  files  and  outputs  a  structure  containing  the  data  on  peaks,  line  center,  ampli¬ 
tude,  Lorentzian  and  Doppler  widths,  and  the  associated  confidence  boundaries  for 
each.  A  similar  routine  is  available  to  extract  the  fitting  data  from  PEAKFlT®data 
export  files. 


function  PeakFitData  =  peakFitNumericParse; 

%% 

%  Created  -  2010-12-14  -  CAR  -  Chris  Rice  -  crice@afit.edu 
%  Adapted  -  2010-12-14  -  CAR  -  From  pollOutputs . m 
%  Adapted  txt  file  parser  for  peakfit  needs. 

%  m-File  that  loads  a  specialized  PeakFit  Data  Output  Text  File. 

%% 

%  Cleanup  Environment: 

%  clear  all; 

%  close  all; 

%  clc; 

%  tic; 


%% 

%  Heavy  Lifting: 

%-+-+-+-+-+-+-+-+-+-+-+ 

%Choose  File  From  A  Directory: 

[txtFileName,  txtPathName]  =  uigetf ile ({' *Numeric* . txt Please  Pick'  ... 

'a  Peakfit  Numeric  File  (Must  have  Numeric  In  Name'); 
disp  (  [txtFileName] ) ; 

%Strings  To  Find: 

strl  =  'Voigt  Area  G/L'; 
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str2  =  'rA2  Coef  Det  DF  Adj  r~2  Fit  Std  Err  F-value ' ; 
%0pen  the  File: 

FID  =  fopen ( [txtPathName, txtFileName] ,  ' rt ' ) ; 

%Pull  in  a  line  from  the  file: 
theLine  =  fgetl (FID) ; 

%Placeholder  to  find  matches: 
numStrMatches  =  0; 

%4  Bundles  of  Data  in  a  PeakFit  File: 
currentDataBundle  =  0; 

%Scan  file,  determine  how  many  peaks  have  been  found: 
while  ( 1 ) 

%Get  Next  Line : 
theLine  =  fgetl (FID) ; 

%Break  at  end  of  file: 
if  (-dschar (theLine) ) 

%Never  Found  The  String: 
break 

end 

%Find  A  String: 

K  =  strfind (theLine, strl )  ; 

%Check  For  A  Match: 
if  (-dsempty  (K)  ) 

%Found  a  Match: 

numStrMatches  =  numStrMatches  +  1 ; 

end 

%Find  A  String: 

K  =  strfind (theLine, str2 )  ; 

%Check  For  A  Match: 
if  (-dsempty  (K)  ) 

%Get  the  next  line: 
theLine  =  fgetl  (FID) ; 

%Found  a  Match: 

fitParams  =  str2num (theLine ) ; 

end 

end 

%Currently  End  of  File,  Rewind: 
f rewind (FID) ; 

%Disp  num  of  peaks: 
numPeaks  =  numStrMatches/4 ; 

disp ( [num2str (numPeaks) , '  Peaks  Found:  Parsing. ..']); 

%Pre-form  arrays: 
bundlel  =  zeros (numPeaks,  4)  ; 
bundle2  =  zeros (numPeaks, 6) ; 
bundle3  =  zeros (numPeaks, 6) ; 
bundle4  =  zeros (numPeaks, 5*4 ) ; 
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while  ( 1 ) 

%Get  Next  Line : 
theLine  =  fgetl (FID) ; 

%Break  at  end  of  file: 
if  (-lischar  (theLine)  ) 

%Never  Found  The  String: 
break 

end 

%Find  A  String: 

K  =  strfind (theLine, strl )  ; 

%Check  for  a  Match: 
if  (-lisempty  (K)  ) 

%Rotate  to  next  bundle  of  data: 
currentDataBundle  =  currentDataBundle+1 ; 
if (currentDataBundle==l ) 

bundlel(l,:)  =  str2num (theLine (24 : end) ) ; 
for (m  =  2:numPeaks) 

theLine  =  fgetl (FID) ; 

bundlel (m, : )  =  str2num (theLine (24 : end) ) ; 

end 

elseif (currentDataBundle==2 ) 

bundle2(l,:)  =  str2num (theLine (24 : end) ) ; 
for (m  =  2:numPeaks) 

theLine  =  fgetl (FID) ; 

bundle2 (m,  : )  =  str2num (theLine (24 : end) )  ; 

end 

elseif  (currentDataBundle==3) 

bundle3(l,:)  =  str2num (theLine (24 : end) ) ; 
for (m  =  2:numPeaks) 

theLine  =  fgetl (FID) ; 

bundle3 (m, : )  =  str2num (theLine (24 : end) ) ; 

end 

elseif (currentDataBundle==4 ) 

theLine  =  fgetl (FID) ; theLine  =  fgetl (FID) ; 
bundle4 ( 1 , 1 : 5)  =  str2num (theLine ( 8 : end) ) ; 
theLine  =  fgetl (FID) ; 

bundle4 ( 1 , 6 : 10 )  =  str2num (theLine ( 8 : end) ) ; 
theLine  =  fgetl (FID) ; 

%Check  if  Doppler  Width  has  been  fixed: 
if ( length ( str2num (theLine ( 8 : end) ) ) ==1 ) 

bundle4 (1, 11 )  =  str2num (theLine ( 8 : end) ) ; 

else 

bundle4 ( 1 , 11 : 15 )  =  str2num (theLine ( 8 : end) ) ; 

end 

theLine  =  fgetl (FID) ; 

bundle4 ( 1 , 1 6 : 20 )  =  str2num (theLine ( 8 : end) ) ; 
for (m  =  2:numPeaks) 

theLine  =  fgetl (FID) ; theLine  =  fgetl (FID) ; theLine 
theLine  =  fgetl (FID) ; 

bundle4 (m, 1 : 5 )  =  str2num (theLine ( 8 : end) ) ; 
theLine  =  fgetl (FID) ; 


fgetl (FID) ; 
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bundle4 (m,  6 : 1 0 )  =  str2num (theLine ( 8 : end) ) ; 
theLine  =  fgetl (FID) ; 

%Check  if  Doppler  Width  has  been  fixed: 
if ( length ( str2num (theLine ( 8 : end) ) ) ==1 ) 

bundle4 (m, 11 )  =  str2num (theLine ( 8 : end) ) ; 

else 

bundle4 (m, 11 : 15)  =  str2num (theLine ( 8 : end) ) ; 

end 

theLine  =  fgetl (FID) ; 

bundle4 (m, 16 : 20 )  =  str2num (theLine ( 8 : end) ) ; 

end 

else 

%Error : 

error ('Error  In  File  Formatting'); 

end 

end 

end 

%Close  the  File: 
f close (FID) ; 

%% 

%  Reorganize  Data: 

%Overall  Fit  Params : 

PeakFitData . rSq  =  f itParams ( 1 ) ; 

PeakFitData . stdErr  =  f itParams ( 3 ) ; 

%0rder  the  File  Output: 

PeakFitData . B1  =  f lipud (bundlel ) ; 

PeakFitData . B2  =  f lipud (bundle2 ) ; 

PeakFitData . B3  =  f lipud (bundle3 ) ; 

PeakFitData . B4  =  f lipud (bundle4 ) ; 

%File  Format: 

PeakFitData . BIFormat  =  { ' aO ' , ' al ' , ' a2 ' , ' a3 ' } ; 

PeakFitData . B2Format  =  {' Amplitude ',' Center ',' FWHM ' ,  ' Asym50 ' ,  ' FW  Base '  ,  ' Asyml 0 ' } ; 
PeakFitData . B3Format  =  { 'Anlytc  Area', '%  Area','Int  Area', '%  Area ',' Cent roid ',.. . 
' Moment 2 ' }  ; 

PeakFitData . B4Format  =  { 'Area' , 'Area  Std . Err ',' Area  t-val', 'Area  95%',... 

'Area  105% ',' Center ',' Center  Std . Err ',' Center  t-val'  ,  'Center  95%',... 

'Center  105%',  'Width  1',  'Width  1  Std . Err ',' Width  1  t-val ',' Width  1  95%',... 
'Width  1  105%', 'Width  2', 'Width  2  Std . Err ',' Width  2  t-val ',' Width  2  95%',... 
'Width  2  105%'}; 

%Save  MatFile: 

save ( [txtFileName (1:19),  ' PeakFit . mat ' ] ,  ' PeakFitData ' ) ; 

%Save  Excel  Spreadsheet: 
warning  off  MATLAB : xlswrite : AddSheet ; 
fName  =  [ txtFileName ( 1 : end-4 ),' -PFD . xls '] ; 
delete (fName) ; 

xlswrite (fName, PeakFitData .BIFormat, ' DataSetl ' , ' Al ' ) ; 
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A1  ' 


xls write ( fName, PeakFitData . B1 ,  ' DataSet 1 '  ,  ' A3 ' )  ; 
xlswrite ( fName, PeakFitData . B2Format , 'DataSet2;, ' 
xls write ( fName,  PeakFitData . B2  ,  ' DataSet 2 ' ,  ' A3 ' ) ; 
xlswrite ( fName, PeakFitData . B3Format , 'DataSet3:, "A1 
xlswrite ( fName, PeakFitData . B3 ,  ' DataSet 3 ' ,  ' A3 ' )  ; 
xlswrite ( fName, PeakFitData . B4 Format , ' DataSet 4 ’ , ' A1 
xlswrite ( fName,  PeakFitData . B4  ,  ' DataSet 4 '  ,  ' A3 ' ) ; 
warning  on  MATLAB : xlswrite : AddSheet ; 

%Trash  Dumb  Default  Sheets: 
fPath  =  pwd;  %  Current  working  directory. 
sheetName  =  'Sheet';  %  EN :  Sheet,  DE :  Tabelle,  etc 
%  Open  Excel  file. 

objExcel  =  actxserver (' Excel . Application ') ; 
objExcel .Workbooks .Open (fullfile (fPath,  fName));  % 
%  Delete  sheets, 
try 

%  Throws  an  error  if  the  sheets  do  not  exist. 
objExcel . ActiveWorkbook . Worksheets . Item ( [ sheetName 
objExcel .ActiveWorkbook .Worksheets . Item ( [sheetName 
objExcel .ActiveWorkbook . Worksheets . Item ( [ sheetName 
catch 

%  Do  nothing, 
end 

%  Save,  close  and  clean  up. 
objExcel .ActiveWorkbook . Save; 
objExcel .ActiveWorkbook . Close; 
objExcel . Quit ; 
objExcel . delete; 

%Done : 

disp ( 'Parse  Complete'); 
return 


) ; 

) ; 

(Lang,  dependent) 

Full  path  is  necessary! 

'1  ] )  . Delete; 

'2  ] )  . Delete; 

'3  ] ) . Delete; 


107 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


Model  for  K  D1  Absorbance 


The  following  code  calculated  the  pressure-broadened  absorption  for  the  K  Z^hyperfine 
lines.  This  code  was  based  off  existing  code  demonstrated  with  cesium.  [5]  A  line  cen¬ 
ter  shift  was  added,  pressure  broadening  values  were  changed,  as  well  as  parameters 
related  to  the  hyperfine  lines  of  K. 


function  [abs]  =  fnAbsCoef fKDl ( f req, T, PN2 , PHe, PCH4 , PC2H6, P3He) 

%% 

%  Adapted  -  2012-04-12  -  CAR  -  Chris  Rice  -  crice@afit.edu  -  From  Monte 
%  Anderson's  Rb  Simulation,  fnAbsCoef fRbDl  and  D2. 

%  Updated  -  2012-04-12  -  CAR  -  Additional  updates,  finding  refs  for  numbers 
%  Updated  -  2012-04-19  -  CAR  -  Finished  rate  updates  &  final  checks. 

%  Calculates  Potassium  D1  absorption  coefficent  (1/cm)  at  a  given  frequency 
%  at  temperature  T  (K)  and  buffer  gas  pressures  (Torr) . 


%%  CONSTANTS 

global  c  eo  h  hbar  kB  AMU; 

global  AmK39  MK39  Inam39  vo39  Iso39  tau; 

global  AmK41  MK41  Inam41  vo41  Iso41; 

global  TransM  numTrans; 

global  dv39S12  dv41S12  vHF  SFF; 


%  Physical  Constants  (NIST) 
c  =  2.997  92458e8 ; 
eo  =  8 . 85418781762e-12; 
h  =  6. 62606896e-34; 
hbar  =  1 . 054571628e-34; 
kB  =  1 ,3806504e-23; 

AMU  =  1 . 660538782e-27; 


%  vacuum  speed  of  light  (m/s) 
%  vacuum  permittivity  (F/m) 

%  Planck  const  (Js) 

%  Planck  const/2pi  (Js) 

%  Boltzmann  const  (J/K) 

%  Atomic  Mass  Unit  (kg) 


%  Potassium  Atomic  Data 
AmK3  9  =  38.96370668; 
AmK4 1  =  40.96182576; 
MK39  =  AmK39*AMU; 

MK4 1  =  AmK41*AMU; 

Inam39  =  3/2; 

Inam41  =  3/2; 

Iso39  =  39; 

Iso4 1  =  41; 


%  Atomic  mass  (AMU)  of  K39 
%  Atomic  mass  (AMU)  of  K41 
%  Mass  of  K39  atom  (kg) 

%  Mass  of  K41  atom  (kg) 

%  Total  nuclear  angular  momentum  K39 
%  Total  nuclear  angular  momentum  K41 
%  Isotope  number  for  K39 
%  Isotope  number  for  K41 


%  K  D1  Data 

vo39  =  389 . 28 60587 1 6el2 ;  %  Frequency  (Hz)  K39 

vo41  =  389 . 286294205el2;  %  Frequency  (Hz)  K41 
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tau  =  26.37e-9; 


%  Lifetime  (s)  K39  and  K41 


%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http :/ /staff . science . uva . nl/-itgtiecke/PotassiumProperties . pdf 


%  Transition  Matrix,  contains  allowed  transitions 


numbers  in  the  format  (Inam,  Jupper, 


%  NOTE: 
TransM 


Jlower  is 
=  [Inam39 
Inam3  9 
Inam3  9 
Inam3  9 
Inam4 1 
Inam4 1 
Inam4 1 
Inam4 1 


1/2  (ground  state) . 
1/2  1  2  I s o 3 9 ; 

1/2  1  2  I s o 3 9 ; 

1/2  2  1  I s o 3 9 ; 

1/2  2  2  I s o 3 9 ; 

1/2  1  1  I so4 1 ; 

1/2  1  2  I so4 1 ; 

1/2  2  1  I so4 1 ; 

1/2  2  2  Iso41] ; 


Flower, 


and  pertenent 
Fupper ) 


numTrans  =  length (TransM) ; 


quantum 


%%  PLACEHOLDER  -  EVERYONE  ABOVE  THIS  LINE  IS  GOOD  TO  GO  FOR  K! 


%%  Hyperfine  Structure  -  Transition  Frequencies  (from  TransM) 
vHF  =  zeros ( 1 , numTrans ) ; 
for  j  =  1: numTrans 

vHF ( j )  =  (hf sEnergyK (TransM ( j , 1 ) , 1 , TransM (j,2) , TransM (j, 4) , TransM (j, 5) ) /h)  -  . . . 
(hf sEnergyK (TransM (j,l), 0,1/2, TransM ( j , 3 ) , TransM ( j  ,  5) ) /h)  +  dvlso (TransM ( j , 5 ) ) ; 

end 

dv39S12  =  hfsEnergyK (Inam39, 0, 1/2, 3, Iso39) /h  -  hfsEnergyK (Inam39, 0, 1/2, 2, Iso39) /h; 
dv41S12  =  hfsEnergyK (Inam41, 0, 1/2, 2, Iso41) /h  -  hfsEnergyK (Inam41, 0, 1/2, 1, Iso41) /h; 

%%  PLACEHOLDER  -  EVERYONE  ABOVE  THIS  LINE  IS  GOOD  TO  GO  FOR  K! 

%%  Line  Strength  Factors,  SFF 
SFF  =  zeros ( 1 , numTrans ) ; 
for (n  =  l:numTrans) 

SFF (n)  =  S_fun(TransM(n,3),TransM(n,4),0.5, TransM (n,  2 )  ,  TransM (n,  1 ) )  ; 

end 

%%  main  function  output 

abs  =  NumDensityK (T) *CrossSectionK (freq, T, PN2,  PHe,  PCH4,  PC2H6,  P3He)  ; 

end  %end  of  f nAbsCoef fKDl  main  function  - 


%%  f nAbsCoef fKDl  SUBFUNCTIONS 


%%  HYPERFINE  ENERGY  LEVELS  - 

function  [dEhfs]  =  hfsEnergyK ( Inam, L, J, F, IsoNum) 
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%  Calculates  hyperfine  structure  energy  for  Cs  D1  and  D2 . 

o, 

o 

%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http :/ /staff . science . uva . nl/-itgtiecke/PotassiumProperties . pdf 
global  h  Iso39  Iso41; 

%Hyperfine  Structure  Constants  for  K39: 


A39S12  =  h*2  30 . 8598601e6;  %(Hz) 

A39P12  =  h*27 . 775e6;  %(Hz) 

A39P32  =  h*6.093e6;  %(Hz) 

B39P32  =  h*2 . 786e6;  %(Hz) 

%Hyperfine  Structure  Constants  for  K41: 
A41S12  =  h*127 . 0069352e6;  %(Hz) 

A41P12  =  h*15 . 245e6;  %(Hz) 

A41P32  =  h*3 . 363e6;  %(Hz) 

B41P32  =  h*3 . 351e6;  %(Hz) 


K  =  F*(F  +  1)  -  Inam* ( Inam  +1)  -  J*(J  +  1); 

%  tempStr  =  fprintf ( ' lnam=%0 . 5g, L=%0 . 5g, J=%0 . 5g, F=%0 . 5g, IsoNum=%0 . 5g ' , Inam, 

%L, J, F, IsoNum) ; 

%  disp (tempStr ) ; 

%  keyboard; 

if (IsoNum  ==  Iso39)  %  K39  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5*A39S12*K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*A39P12*K; 
elseif(L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A39P32*K)  +  B39P32* ( ( 1 . 5*K* (K+l )  -  2*Inam* ( Inam+1 ) * J* . . . 
(J+l)  ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) ) ; 

else 

tempStr  =  fprintf (' lnam=%0 . 5g, L=%0 . 5g, J=%0 . 5g, F=%0 . 5g, IsoNum=%0 . 5g ' ,  ... 

Inam, L, J, F, IsoNum)  ; 
disp (tempStr) ; 

tempStr  =  fprintf (' IMPROPER  INPUTS  -  Hyperfine  energy  level'  ... 

' determination ' )  ; 
disp (tempStr) ; 
dEhfs  =  0; 

end 

elseif ( IsoNum  ==  Iso41)  %  K41  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5  *A4 1S12  *K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5  *A4 IP 12  *K; 
elseif (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A41P32*K)  +  B4 1P32* ( ( 1 . 5*K* (K+l )  -  2*Inam* ( Inam+1 ) *J* .. . 
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(J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) ) ; 

else 

tempStr  =  f print f ( ’ lnam=%0 . 5g, L=%0 .  5g, J=%0 .  5g, F=%0 .  5g, IsoNum=%0 . 5g ' ,  . . . 

Inam, L, J , F, IsoNum) ; 
disp (tempStr) ; 

tempStr  =  fprintf (’ IMPROPER  INPUTS  -  Hyperfine  energy  level  '  ... 

1  determination  1 )  ; 
disp (tempStr) ; 
dEhfs  =  0; 

end 

else 

tempStr  =  fprintf ( 1 lnam=%0 . 5g, L=%0 . 5g, J=%0 . 5g, F=%0 . 5g, IsoNum=%0 . 5g ’ ,  . . . 

Inam, L, J, F, IsoNum) ; 
disp  (tempStr) ; 

tempStr  =  fprintf (' IMPROPER  INPUTS  -  Hyperfine  energy  level  determination'); 
disp  (tempStr) ; 
dEhfs  =  0; 

end 

end  %end  of  function 

%Hyperfine  subfunction 
function  dvisoOut  =  dvlso ( IsoNum) 

%  Provides  Isotope  shift  for  K39  -  all  frequencies  relative  to  K41. 
global  vo39  vo41  Iso39  Iso41; 

dvo  =  vo41  -  vo39;  %  Isotope's  hyperfine  spectrum  rel .  to  same 

%  central  frquency  vo  so  that 
%  voK39  +  dviso  =  voK41 

if (IsoNum  ==  Iso39) 
dvisoOut  -  dvo; 
elseif ( IsoNum  ==  Iso41) 
dvisoOut  =  0; 

else 

fprintf (' \nERROR  —  Isotope  Number  is  outside  allowed  range . \n\n ') ; 

end 

end 


%%  LINE  STRENGTH  FACTORS  - 

function  [S]  =  S_f un (F, Fprime, J, Jprime, I ) 
SI  =  Wigner6j (J, Jprime, 1, Fprime, F, I) ; 

S  -  (2*Fprime+l) * (2*J+1) * (Sl~2)  ; 

end  %end  of  function 


%%  WIGNER  6 j  FACTORS - 

function  Wigner  =  Wigner6 j ( jl, j2, j3, Jl, J2, J3) 

%  Wigner  6j-symbol  calculator.  Written  by  Amita  B  Deb,  Clarendon  Lab.  2007. 
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%  Improved  by  Richard  A.  Holt,  Univ.  of  Western  Ontario,  2009. 

%  Further  improved  to  deal  with  large  input  arguments  by  Lee  Harper, 

%  CRL  Oxford,  2009. 

%  Builds  on  the  m-file  Wigner 6 jcoef f  posted  by  Richard  A.  Holt. 

%  Calculates  {  jl,  j2  , j 3 }  using  Racah  formula.  See:  Sobelman:  Atomic 
%Spectra  and  Radiative  Transitions. 

%  Jl  J2  J3 

%  Finding  Triangular  coefficients 

tril  =  triangle.coef f ( jl, j2, j3) ; 
tri2  =  triangle.coef f ( jl, J2, J3) ; 
tri3  =  triangle.coef f ( Jl, j2, J3) ; 
tri4  =  triangle.coef f ( Jl, J2, j3) ; 


if  (tril  =  =  0  1 1  tri2  =  =  0  | 
Wigner=0 ; 
return 

end 


tri3  =  =  0||tri4==0) 


Finding  the  range  of  summation  in  the  Racah  formula. 


a  ( 1 ) 
a  (2 ) 
a  (3) 
a  ( 4 ) 


ji 

ji 

ji 

ji 


j2 

J2 

j2 

J2 


j3; 
J3 ; 
J3 ; 
j3; 


rangei  =  max (a); 


k ( 1 )  =  jl  +  j2  +  Jl  +  J2 ; 

k(2)  =  j2  +  j3  +  J2  +  J3; 

k ( 3 )  =  j3  +  jl  +  J3  +  Jl; 

rangef  =  min(k); 

Wigner  =  0; 

for  t=rangei : rangef 

Wigner  =  Wigner  +  (  (-1 )  ~t ) *exp (gammaln (t  +  2 ) 


fung  (t,  jl,  j2,  j3,  Jl,  J2,  J3)  )  ; 


end 

Wigner  =  (tril*tri2*tri3*tri4) " (0 . 5) *Wigner; 


end 

%  -  Wigner6j  subfunctions  - 

function  r  =  fung(t,  jl,  j2,  j3, Jl, J2, J3) 

%  Calculating  the  logarithm  of  the  denominator  in  Racah  Formula,  using 
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%  the  gamma  function  in  place  of  the  factorial. 

r  =  sum (gamma In ( [ (t-jl-j2-j3) ;  (t- jl-J2-J3) ;  (t-Jl- j2-J3) ;  (t-Jl-J2- j  3 ) ;  .  .  . 

( jl  +  j2+Jl+J2-t) ;  ( j2  + j3+J2+J3-t) ;  ( j3  +  j 1+ J3+ Jl-t ) ]  +  1)  )  ; 


end 

o, _ 

o 

function  tri  =  triangle_coef f (a,  b,  c) 

%  Calculates  triangle  coefficients  for  angular  momenta. 

%  This  version  returns  0  if  the  triangle  inequalities  are  violated.  (RAH) 

if  (a<0  |  |  b<0  |  |  c<0) 

tri=0 ; 
return 

end 

for  xa  =  abs (a-b) : 1 : a+b 
if  c==xa 

tri  =  factorial (a+b-c) *factorial (a-b+c) ^factorial (-a+b+c) /  (  ... 

factorial (a+b+c+1) ) ; 
return 

end 

end 

tri=0 ; 

end  %end  of  function 


%%  CROSS  SECTION  AT  LINE  CENTER  - 

function  [sO]  =  sigma_0K (SFF, vD, x) 

%  Calculates  cross  section  at  line  center  for  each  hyperfine  transition 
global  c  vo39  vo41  tau  TransM  Iso39  Iso41; 

if  (TransM (x,5)==Iso39) 
vo  =  vo39; 

el seif (TransM (x, 5) ==Iso41) 
vo  =  vo41; 

else 

fprintf ( ' \nERR0R  in  sigma_0K  -  -  Inam  is  outside  allowed  range . \n\n ') ; 
return ; 

end 

dvnat  =  1/  (2*pi*tau) ; 

sO  =  SFF (x) * (1/tau) * ( ( (c*100) ~2) / (8*pi) ) * (1/ ( (vo+vD (x) ) ~2) )*. . . 

( (2*TransM (x, 2) +1) / (2*. 5  +  1) ) * (2/ (pi* dvnat ) )  ; 


end 


%%  NUMBER  DENSITY 
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function  [ndensity]  =  NumDensityK (T) 

%  calculates  number  density  (crrU-3)  given  Temperature  in  Kelvin 
global  kB; 

ndensity  =  (KVaporPressure (T) )/ (kB*T) *101325* (le-6) ; 
end 


%%  VAPOR  PRESSURE  - 

function  [Pv]  =  KVaporPressure (T) 

%  Provides  K  vapor  pressure  in  atmospheres  given  T  in  Kelvin 

o, 

o 

%  REFERENCE: 

%  C.B.Alcock,  V.P.Itkin,  and  M.K.Horrigan 

%  Vapour  pressure  equations  for  the  metallic  elements:  298-2500K 
%  Canadian  Metallurgical  Quarterly,  1984,  23,  309-313 

o, 

o 

%  "Precise"  equation  for  K  vapour  pressure. 

%  Uses  melting  point  for  K  at  336. 8K  (63.65  degC) 

%  Valid  for  temperatures  from  melting  point  to  600K. 

o, 

o 

%  For  best  precision,  use  temperatures  between  melting  point  and  600K. 

if (T  >  336.8  &&  T  <  600) 
a  =  8.233; 
b  =  -4693; 
c  =  -1.2403; 

Pv  =  10. "(a  +  b*(T~-l)  +  c*loglO (T) )  ; 
elseif  (T  >  298  &&  T  <  336.8) 
a  =  4.961; 
b  =  -4646; 

Pv  =  10 . ~  (a  +  b* ( T  ~  — 1 ) ) ; 
elseif  (T  <  298  | |  T  >  600) 

fprint f ( ' \nlemperature  is  outside  of  allowed  range. \n\n') 
return 

end 

end 


%%  CROSS  SECTION  - 

function  [cs]  =  CrossSectionK (x, T, PN2 , PHe, PCH4 , PC2H6 , P3He) 

%  Calculates  cross  section  as  function  of  frequency,  x 
global  tau  sO  vo39  vo41  TransM  vHF  numTrans  Iso39  Iso41  c  SFF; 

dvnat  =  1/ (2*pi*tau) ;  %  natural  linewidth,  a  nu  natural 

csHF  =  zeros (1, numTrans) ;  %  Cross  section  for  an  individual  HF  transition 

temp  =  0; 

for  n  =  1: numTrans 

Flwr  =  TransM(n,3); 

IsoNum  =  TransM(n,5); 
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[ gamma ,  shift, dvD, dvL]  =  Buf ferGasK (T, IsoNum,  PN2, PHe,  PCH4,  PC2H6,  P3He)  ; 
if(IsoNum  ==  Iso39)  %  Set  central  frequency  based  on  isotope  (IsoNum) 

vo  =  vo3  9; 

elseif (IsoNum  ==  Iso41) 
vo  =  vo41; 

else 

fprint f ( ' \nERROR  in  CrossSect ionK-Isotope  Number  is  outside  allowed  range. \n\n' 
return; 

end 

sO  =  SFF (n) * (1/tau) * ( ( (c*100) ~2) / (4*pi) ) * (1/ ( (vo+vHF (n) + shift) *2) ) * . . . 

( (2*TransM (n, 2 ) +1 ) / (2*. 5+1) ) ; 

csHF  =  KNat Abundance (IsoNum)  *  BoltzRatioK (IsoNum, Flwr, T)  *  sO  *  ... 

( (x  +  vHF (n)  +  shift  +  vo)  /  (vHF (n)  +  shift  +  vo) )  *  . . . 
gVoigt ( (x  -  vHF (n) -shift ), dvL, dvD) ; 
temp  =  temp  +  csHF; 

end 

cs  =  temp; 
end 

o.  o. 
o  o 

%  Need  Total  Nuclear  Ang.  Momentum 
function  eta  =  KNatAbundance ( IsoNum) 

%  Return  Potassium  natural  abundance  based  on  Total  nuclear  angular  momentum: 
global  Iso39  Iso41; 

%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http : //staff . science . uva . nl/-itgtiecke/PotassiumProperties . pdf 

%Natural  Abundance  Ratios: 

Eta39  =  0.932581;  %Nuclear  Spin  =  3/2 
Eta41  =  0.067302;  %Nuclear  Spin  =  3/2 

%Note:  Eta40  =  0.000117;  %Nuclear  Spin  =  4  and  is  ignored  here  b/c  of  low  natural 
% abundance . 

if  (IsoNum  ==  Iso39) 
eta  =  Eta39; 
return ; 

elseif  (IsoNum  ==  Iso41) 
eta  =  Eta41; 
return ; 

else 

fprintf ( ' \nIsotope  Number  is  outside  of  allowed  range  for  K39  or  K41.\n\n'); 
return ; 

end 
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end 


%%  BUFFER  GAS  EFFECTS  - 

function  [gamma, shift , dvD, dvL]  =  Buf ferGasK (T, IsoNum, PN2 , PHe, PCH4 , PC2H6, P3He) 
%  Rb  buffer  gas  effects  -  pressure  broadening  and  frequency  shift 

o, 

o 

%  REFERENCE: 

%  Pitz,  G.  A. 

%  Pressure  broadening  and  shift  of  the  potassium  D1  transition  by  the 
%  noble  gasses  and  N2,  H2,  HD,  D2,  CH4 ,  C2H6,  C38,  and  n-C4H10 
%  JQSRT,  113,  387-395  2012. 

global  tau  c  kB  Iso39  Iso41  MK39  MK41  vo39  vo41; 


%Calculate  Broadening: 


gamma  = 

(17 . 78e6) * ( (T/328) ~  (1/2) ; 

1  *PN2 

o. 

.  o 

MHz/Torr 

(Pitz) 

+ 

(13 . 08e6) * ( (T/328) ~  (1/2) ; 

1  *PHe 

o, 

,  o 

MHz/Torr 

(Pitz) 

+ 

(29 . 35e6) * ( (T/328) ~  (1/2) ; 

1  *PCH4  .  .  . 

o_ 

.  o 

MHz/Torr 

(Pitz) 

+ 

(26. 63e6) * ( (T/328) ~  (1/2) ; 

1  *PC2H6 .  .  . 

o, 

,  o 

MHz/Torr 

(Pitz) 

+ 

(17 . 46e6) * ( (T/328) ~  (1/2) ; 

1 *P3He; 

o_ 

o 

MHz/Torr 

(Pitz) 

%  gamma 

=  0;  %for  test  purposes 

to  check 

the  buffergas  shift 

%Calculate  Shift: 

shift  = 

(-6 . 80e6) * ( (T/328) ~  (1/2) ; 

1  *PN2 

o. 

.  o 

MHz/Torr 

(Pitz) 

+ 

(  1 . 63e6) *( (T/328) ~  (1/2) ; 

1  *PHe 

o. 

.  o 

MHz/Torr 

(Pitz) 

+ 

(-7 . 41e6) * ( (T/328) ~  (1/2) ; 

1 *PCH4  . . . 

o. 

.  o 

MHz/Torr 

(Pitz) 

+ 

(-8 . 32e6) * ( (T/328)  ~  (1/2) ; 

1 *PC2H6 . . . 

o, 

,  o 

MHz/Torr 

(Pitz) 

+ 

(  6.82e6)*( (T/328)  "  (1/2) ; 

1 *P3He; 

o, 

o 

MHz/Torr 

(Pitz) 

%  Set  central  frequency  and  mass  based  on  isotope 
if  (IsoNum  ==  Iso39) 
vo  =  vo39; 

MK  =  MK39 ; 

elseif ( IsoNum  ==  Iso41) 
vo  =  vo41; 

MK  =  MK4 1 ; 

else 

fprintf ( ' \nERR0R  in  Buf ferGasK-Iso  is  outside  range . \n\n ') ; 
return ; 

end 

%Calculate  Doppler  and  Lorentz: 

dvD  =  (vo  +  shift) * ( (8*kB*T*log (2) ) / (MK* (c~2) ) ) "  (1/2) ; 
dvL  =  (1/ (2*pi*tau) )  +  gamma; 

end 

%%  BOLTZMAN  RATIO  - 

function  [BR]  =  BoltzRatioK ( IsoNum, F, T) 
global  kB  h  dv39S12  dv41S12  Iso39  Iso41; 

if (IsoNum  ==  Iso39) 


116 


456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 


if (F  ==  1) 

BR  =  5-/(5  +  7* exp (- (h*dv39S12)  . /  (kB . *T) ) ) ; 
elseif (F  ==  2) 

BR  =  (7*exp (- (h*dv39S12)  ./ (kB.*T) ) )  ./ (5  +  7 *exp (- (h*dv3 9S12 )  . / (kB . *T) ) ) ; 
elseif (F  #  1  &&  F  #  2) 

fprintf ( 1 \nIMPROPER  INPUTS-4 55\n\n ' ) 

BR  =  0; 

end 

elseif  ( IsoNum  ==  Iso41) 
if (F  ==  1) 

BR  =  3 . / (3  +  5* exp (- (h*dv41S12 )  . /  (kB . *T) ) ) ; 
elseif (F  ==  2) 

BR  -  (5*exp (- (h*dv41S12)  . / (kB. *T) ) )  . / (3  +  5*exp (- (h*dv4 1S12 )  . / (kB . *T) ) ) ; 
elseif (F  #  1  &&  F  #  2) 

fprintf ( ' \nIMPROPER  INPUTS-4  64\n\n ' ) 

BR  =  0 ; 

end 

else 

fprintf ( ' \nIMPROPER  INPUTS-4 6 8\n\n ' ) 

BR  =  0 ; 

end 

end 


%%  VOIGT  LINESHAPE  - 

function  [y]  =  gVoigt (v, dvL, dvD) 

%  A  useful  approximation  for  the  Voigt  convolution  is  given  by  Thompson 

o, 

o 

%  REFERENCE: 

%  William  J.  Thompson 

%  Numerous  neat  algorithms  for  the  Voigt  profile  function. 

%  Computers  in  Physics,  Vol  7,  No.  6,  Nov/Dec  1993. 

o. 

o 

%  gV(u,a)  =  Re [Exp (z~2) erfc (z) ] 

%  z  =  a  +  iu 

%  a  =  sqrt(ln(2))  vLorentzian/vDoppler 

%  u  =  2  sqrt(ln(2))  (v-vo) /vDoppler 

%  NOTE:  Not  valid  for  Doppler  free  spectra  where  vDoppler  =  0 
%  NOTE:  erfc(z)  =  1  -  erf(z)  (related  to  the  Faddeeva  function) 

a  =  sqrt (log (2) ) * (dvL/dvD) ; 
u  =  2*sqrt (log (2) ) * (v/dvD)  ; 
z  =  a  +  li*u; 

y  =  2*sqrt ( log (2 ) /pi )*( 1/dvD) *real ( faddeeva ( li*z ))  ; 
end 

%  -  gVoigt  subfunctions  - 

function  w  =  f addeeva ( z ,  N) 

%  FADDEEVA  Faddeeva  function 
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o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 
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o, 
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o 
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o 


W  =  FADDEEVA(Z)  is  the  Faddeeva  function,  aka  the  plasma  dispersion 
function,  for  each  element  of  Z.  The  Faddeeva  function  is  defined  as: 

w(z)  =  exp(-z~2)  *  erfc(-j*z) 

where  erfc(x)  is  the  complex  complementary  error  function. 

W  =  FADDEEVA (Z, N)  can  be  used  to  explicitly  specify  the  number  of  terms 
to  truncate  the  expansion  (see  (13)  in  [1]) .  N  =  16  is  used  as  default. 

Example : 

x  =  linspace (-10, 10, 1001) ;  [X,Y]  =  meshgrid (x, x) ; 

W  =  faddeeva (complex (X, Y) ) ; 
figure; 

subplot (121)  ;  imagesc (x, x, real (W) ) ;  axis  xy  square;  caxis ( [-1  1]); 
title (' re (faddeeva ( z ) ) ' ) ;  xlabel ( ' re ( z ) ' ) ;  ylabel ( ' im (z ) ' ) ; 
subplot ( 122 ) ;  imagesc (x, x, imag (W) ) ;  axis  xy  square;  caxis ( [-1  1]); 
title ( ' im (faddeeva ( z ) ) ' ) ;  xlabel ('re(z) ' ) ;  ylabel ( ' im ( z ) ' ) ; 

Reference : 

[1]  J.A.C.  Weideman,  "Computation  of  the  Complex  Error  Function,"  SIAM 
J.  Numerical  Analysis,  pp.  1497-1518,  No.  5,  Vol .  31,  Oct.,  1994 
Available  Online:  http://www.jstor.org/stable/2158232 


if  nargin<2,  N  =  [];  end 
if  isempty(N),  N  =  16;  end 


w  =  zeros (size (z) ) ;  %  initialize  output 


0,0,0,0  0, 
o  o  o  o  o 

%  for  purely  imaginary-valued  inputs,  use  erf  as  is  if  z  is  real 
idx  =  real(z)==0;  % 

w(idx)  =  exp (-z (idx) . ~2)  . *erfc (imag (z (idx) ) ) ; 

if  all (idx),  return;  end 
idx  =  -iidx; 


o,  o,  o,  o,  o, 
o  o  o  o  o 

%  for  complex-valued  inputs 

%  make  sure  all  points  are  in  the  upper  half-plane  (positive  imag.  values) 
idxl  =  idx  &  imag(z)<0; 
z (idxl)  =  con j (z (idxl) ) ; 

M  =  2  *N; 

M2  =  2  *M; 

k  =  (-M+1:1:M-1) ';  %  M2  =  no .  of  sampling  points. 

L  =  sqrt (N/sqrt ( 2 ) ) ;  %  Optimal  choice  of  L. 

theta  =  k*pi/M; 

t  =  L*tan (theta/2 ) ;  %  Variables  theta  and  t. 
f  =  exp  (-t .  "2)  . * (L~2+t . ~ 2) ; 
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f  =  [0;  f] ;  %  Function  to  be  transformed. 

a  =  real (f ft (f ftshift (f ) ) ) /M2;  %  Coefficients  of  transform, 
a  =  f lipud (a (2 : N+l ) ) ;  %  Reorder  coefficients. 

Z  =  (L+li*z (idx) )  . / (L-li*z (idx) )  ; 
p  =  polyval (a, Z ) ;  %  Polynomial  evaluation. 

w(idx)  =  2 *p . / (L-li* z ( idx) ) . ~ 2  +  (1/sqrt (pi) ) . / (L-li*z (idx) ) ; 

%  convert  the  upper  half-plane  results  to  the  lower  half-plane 
w(idxl)  =  con j (2*exp (-z (idxl ) . "2 )  -  w(idxl)); 

end 


Evaluate  w ( z ) . 
if  necessary 
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Model  for  K  D2  Absorbance 


The  following  code  calculated  the  pressure-broadened  absorption  for  the  K  D2 
hyperfine  lines.  This  code  was  based  off  existing  code  demonstrated  with  cesium.  [5] 
A  line  center  shift  was  added,  pressure  broadening  values  were  changed,  as  well  as 
parameters  related  to  the  hyperfine  lines  of  K. 


function  [abs]  =  fnAbsCoef fKD2 ( f req, T, PN2 , PHe, PCH4 , PC2H6, P3He) 

%% 

%  Adapted  -  2012-04-12  -  CAR  -  Chris  Rice  -  crice@afit.edu  -  From  Monte 
%  Anderson's  Rb  Simulation,  fnAbsCoef fRbDl  and  D2. 

%  Updated  -  2012-04-12  -  CAR  -  Additional  updates,  finding  refs  for  numbers 
%  Updated  -  2012-04-19  -  CAR  -  Finished  rate  updates  &  final  checks. 

%  Calculates  Potassium  D2  absorption  coefficent  (1/cm)  at  a  given  frequency 
%  at  temperature  T  (K)  and  buffer  gas  pressures  (Torr) . 


%%  CONSTANTS 

global  c  eo  h  hbar  kB  AMU; 

global  AmK39  MK39  Inam39  vo39  Iso39  tau; 

global  AmK41  MK41  Inam41  vo41  Iso41; 

global  TransM  numTrans; 

global  dv39S12  dv41S12  vHF  SFF; 


%  Physical  Constants  (NIST) 
c  =  2.997  92458e8 ; 
eo  =  8 . 85418781762e-12; 
h  =  6. 62606896e-34; 
hbar  =  1 . 054571628e-34; 
kB  =  1 ,3806504e-23; 

AMU  =  1 . 660538782e-27; 


%  vacuum  speed  of  light  (m/s) 
%  vacuum  permittivity  (F/m) 

%  Planck  const  (Js) 

%  Planck  const/2pi  (Js) 

%  Boltzmann  const  (J/K) 

%  Atomic  Mass  Unit  (kg) 


%  Potassium  Atomic  Data 
AmK3  9  =  38.96370668; 
AmK4 1  =  40.96182576; 
MK39  =  AmK39*AMU; 

MK4 1  =  AmK41*AMU; 

Inam39  =  3/2; 

Inam41  =  3/2; 

Iso39  =  39; 

Iso4 1  =  41; 


%  Atomic  mass  (AMU)  of  K39 
%  Atomic  mass  (AMU)  of  K41 
%  Mass  of  K39  atom  (kg) 

%  Mass  of  K41  atom  (kg) 

%  Total  nuclear  angular  momentum  K39 
%  Total  nuclear  angular  momentum  K41 
%  Isotope  number  for  K39 
%  Isotope  number  for  K41 


%  K  D2  Data 

vo39  =  391.01617003el2; 
vo4 1  =  391.01640621el2; 


%  Frequency  (Hz)  K39 
%  Frequency  (Hz)  K41 
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tau  =  26.37e-9; 


%  Lifetime  (s)  K39  and  K41 


%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http  :  /  /staff  .  science  .  uva  .  nl/-itgtiecke/PotassiumProperties  .  pdf 


Transition  Matrix,  contains  allowed  transitions 
numbers  in  the  format  (Inam,  Jupper,  Flower, 
NOTE: 


Jlower  is 

i  1/2 

(ground 

[Inam39 

3/2 

1 

0 

Iso39 

Inam3  9 

3/2 

1 

1 

Iso39 

Inam3  9 

3/2 

1 

2 

Iso39 

Inam3  9 

3/2 

2 

1 

Iso39 

Inam3  9 

3/2 

2 

2 

Iso39 

Inam3  9 

3/2 

2 

3 

Iso39 

Inam4 1 

3/2 

1 

0 

Iso41 

Inam4 1 

3/2 

1 

1 

Iso41 

Inam4 1 

3/2 

1 

2 

Iso41 

Inam4 1 

3/2 

2 

1 

Iso41 

Inam4 1 

3/2 

2 

2 

Iso41 

Inam4 1 

3/2 

2 

3 

Iso41 

numTrans  =  length (TransM); 


and  pertenent 
Fupper,  IsoNum) 


quantum 


%%  Hyperfine  Structure  -  Transition  Frequencies  (from  TransM) 
vHF  =  zeros ( 1 , numTrans ) ; 
for  j  =  1: numTrans 

vHF ( j )  =  (hf sEnergyK (TransM ( j , 1 ) , 1 , TransM ( j , 2 ) , TransM ( j ,  4 )  ,  TransM (j,5)  )  /h)-. . . 
(hf sEnergyK (TransM (j,l), 0,1/2, TransM ( j , 3) , TransM ( j , 5) ) /h)  +  dvlso (TransM ( j , 5 ) ) ; 
end 

dv39S12  =  hfsEnergyK (Inam39, 0, 1/2, 3, Iso39) /h  -  hfsEnergyK (Inam39, 0, 1/2, 2, Iso39) /h; 
dv41S12  =  hfsEnergyK (Inam41, 0, 1/2, 2, Iso41) /h  -  hfsEnergyK (Inam41, 0, 1/2, 1, Iso41) /h; 

%%  PLACEHOLDER  -  EVERYONE  ABOVE  THIS  LINE  IS  GOOD  TO  GO  FOR  K! 


%%  Line  Strength  Factors,  SFF 
SFF  =  zeros ( 1 , numTrans ) ; 
for (n  =  l:numTrans) 

SFF (n)  =  S_f un (TransM (n, 3 ) , TransM (n, 4 ) , 0 . 5, TransM (n, 2 ) , TransM (n, 1 ) ) ; 

end 


%%  main  function  output 

abs  -  NumDensityK (T) * Cross Sect ionK (freq, T, PN2, PHe, PCH4, PC2H6, P3He) ; 
end  %end  of  f nAbsCoef f KD2  main  function - 


%%  f nAbsCoef fKD 2  SUBFUNCTIONS 

%%  HYPERFINE  ENERGY  LEVELS  — 
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function  [dEhfs]  =  hf sEnergyK ( Inam, L, J, F, IsoNum) 

%  Calculates  hyperfine  structure  energy  for  Cs  D1  and  D2 . 

o. 

o 

%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http :/ /staff . science . uva . nl/->tgtiecke/PotassiumProperties . pdf 

global  h  Iso39  Iso41; 

%Hyperfine  Structure  Constants  for  K39: 


A39S12  =  h*230 . 8598601e6;  %(Hz) 

A39P12  =  h*27 . 775e6;  %(Hz) 

A39P32  =  h*6.093e6;  %(Hz) 

B39P32  =  h*2 . 786e6;  %(Hz) 

%Hyperfine  Structure  Constants  for  K41: 
A41S12  =  h*12  7 . 0069352e6;  %(Hz) 

A41P12  =  h*15 . 245e6;  %(Hz) 

A41P32  =  h*3 . 363e6;  %(Hz) 

B41P32  =  h*3 . 351e6;  %(Hz) 


K  =  F*(F  +  1)  -  Inam* (Inam  +1)  -  J*(J  +  1); 


if (IsoNum  ==  Iso39)  %  K39  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5*A39S12*K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*A39P12*K; 
elseif (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A39P32*K)  +  B39P32* ( ( 1 . 5*K* (K+l )  -  2*Inam* ( Inam+1 ) * J* ( J+l ) ) 

/ (4*Inam* (2*Inam-l) *J* (2*J-1) ) )  ; 

else 

tempS tr  =  fprintf ( ' lnam=%0 .5g,L=%0.5g,J=%0.5g,F=%0.5g, IsoNum=%0 . 5g ' , . . . 

Inam, L, J, F, IsoNum)  ; 
disp (tempStr ) ; 

tempStr  =  fprintf (' IMPROPER  INPUTS  -  Hyperfine  energy  level  '  ... 

' determination '  )  ; 
disp (tempStr) ; 
dEhfs  =  0; 

end 

elseif ( IsoNum  ==  Iso41)  %  K41  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5  *A4 1S12  *K; 
elseif (L  ==  1  &  &  J  ==  .5) 
dEhfs  =  . 5  *A4 1P12  *K; 
elseif(L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A41P32*K)  +  B4 1P32* ( ( 1 . 5*K* (K+l )  -  2*Inam* ( Inam+1 ) *J*  ... 
(J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) )  ; 

else 

tempStr  =  fprintf (' lnam=%0 . 5g, L=%0 . 5g, J=%0 . 5g, F=%0 . 5g, IsoNum=%0 . 5g ' ,  ... 

Inam, L, J, F, IsoNum)  ; 
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disp (tempStr) ; 

tempStr  =  fprintf (’ IMPROPER  INPUTS  -  Hyperfine  energy  level 
' determination ' )  ; 
disp (tempStr) ; 
dEhfs  =  0; 

end 

else 

tempStr  =  fprintf ( ' lnam=%0 . 5g, L=%0 . 5g, J=%0.5g,F=%0.5g,  IsoNum=%0 . 5g '  ,  .  .  . 

Inam, L, J, F , IsoNum) ; 
disp (tempStr) ; 

tempStr  =  fprintf (' IMPROPER  INPUTS  -  Hyperfine  energy  level  determination'); 
disp  (tempStr) ; 
dEhfs  =  0; 

end 

end  %end  of  function 

%Hyperfine  subfunction 
function  dvisoOut  =  dvlso ( IsoNum) 

%  Provides  Isotope  shift  for  K39  -  all  frequencies  relative  to  K41. 
global  vo39  vo41  Iso39  Iso41; 

dvo  =  vo41  -  vo39;  %  Isotope's  hyperfine  spectrum  rel .  to  same 

%  central  frquency  vo  so  that 
%  voK39  +  dviso  =  voK41 

if (IsoNum  ==  Iso39) 
dvisoOut  -  dvo; 
elseif ( IsoNum  ==  Iso41) 
dvisoOut  =  0; 

else 

fprintf (' \nERROR  —  Isotope  Number  is  outside  allowed  range . \n\n ') ; 

end 

end 


%%  LINE  STRENGTH  FACTORS  - 

function  [S]  =  S_f un (F, Fprime, J, Jprime, I ) 
SI  =  Wigner6 j (J, Jprime, 1, Fprime, F,  I) ; 

S  =  (2*Fprime  +  l) * (2*J+1) * (Sl~2)  ; 

end  %end  of  function 


%%  WIGNER  6 j  FACTORS - 

function  Wigner  =  Wigner6 j ( jl, j2, j3, Jl, J2, J3) 

%  Wigner  6j-symbol  calculator.  Written  by  Amita  B  Deb,  Clarendon  Lab.  2007. 
%  Improved  by  Richard  A.  Holt,  Univ.  of  Western  Ontario,  2009. 

%  Further  improved  to  deal  with  large  input  arguments  by  Lee  Harper, 

%  CRL  Oxford,  2009. 
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%  Builds  on  the  m-file  Wigner 6 jcoef f  posted  by  Richard  A.  Holt. 


%  Calculates  {  jl,  j2  , j 3 }  using  Racah  formula.  See:  Sobelman:  Atomic 
%Spectra  and  Radiative  Transitions. 

%  Jl  J2  J3 

%  Finding  Triangular  coefficients 

tril  =  triangle.coef f ( jl, j2, j3) ; 
tri2  =  triangle.coef f ( jl, J2, J3) ; 
tri3  =  triangle.coef f ( Jl, j2, J3) ; 
tri4  =  triangle.coef f ( Jl, J2, j3) ; 

if  (tril  =  =  0  |  |tri2=  =  0|  |tri3  =  =  0  |  |tri4==0) 

Wigner=0 ; 
return 

end 


%  Finding  the  range  of  summation  in  the  Racah  formula. 


a ( 1 )  =  jl  +  j  2  +  j  3 ; 

a  (2 )  =  jl  +  J2  +  J3 ; 

a  (3)  =  Jl  +  j2  +  J3 ; 

a  (4)  =  Jl  +  J2  +  j3; 

rangei  =  max (a); 


k  ( 1 )  =  jl  +  j2  +  Jl  +  J2; 
k(2)  =  j2  +  j3  +  J2  +  J3; 
k  (3)  =  j3  +  jl  +  J3  +  Jl; 


rangef  =  min(k); 


Wigner  =  0; 


for  t=rangei : rangef 


Wigner  =  Wigner  +  ( (-1 ) "t ) *exp (gammaln (t+2 )  -  fung (t, jl, j2, j3, Jl, J2, J3) ) ; 


end 

Wigner  =  (tril*tri2*tri3*tri4 ) " (0 . 5) *wigner; 


end 

% - Wigner 6 j  subfunctions - 

function  r  -  fung(t,  jl, j2, j3, Jl, J2, J3) 

%  Calculating  the  logarithm  of  the  denominator  in  Racah  Formula,  using 
%  the  gamma  function  in  place  of  the  factorial. 

r  =  sum(gammaln ( [ (t  —  jl  —  j2  —  j3) ;  (t-jl-J2-J3) ;  (t-Jl-j2-J3) ;  (t-Jl-J2-j3) ;  .  .  . 

( jl  +  j2+Jl+J2-t) ;  ( j2  + j3  +  J2+J3-t) ;  ( j3  + jl+J3+Jl-t) ]  +  1) )  ; 
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end 


o_ _ 

o 

function  tri  =  triangle_coef f (a, b, c) 

%  Calculates  triangle  coefficients  for  angular  momenta. 

%  This  version  returns  0  if  the  triangle  inequalities  are  violated.  (RAH) 

if  (a<0  | |  b<0  | |  c<0) 

tri=0 ; 
return 

end 

for  xa  =  abs (a-b) : 1 : a+b 
if  c==xa 

tri  =  factorial (a+b-c) *factorial (a-b+c) *factorial (-a+b+c) /  ... 

(factorial (a+b+c+1) ) ; 
return 

end 

end 

tri=0 ; 

end  %end  of  function 


%%  CROSS  SECTION  AT  LINE  CENTER  - 

function  [sO]  =  sigma_0K (SFF, vD, x) 

%  Calculates  cross  section  at  line  center  for  each  hyperfine  transition 
global  c  vo39  vo41  tau  TransM  Iso39  Iso41; 

if (TransM (x, 5) ==Iso39) 
vo  =  vo39; 

el seif (TransM (x, 5) ==Iso41) 
vo  =  vo41; 

else 

fprintf ( ' \nERR0R  in  sigma_0K  -  -  Inam  is  outside  allowed  range . \n\n ')  ; 
return ; 

end 

dvnat  =  1/ (2*pi*tau) ; 

sO  =  SFF (x) * (1/tau) * ( ( (c*100)  ~2) / (8*pi) ) * (1/ ( (vo+vD (x) ) "2) ) *.  .  . 

( (2*TransM (x, 2) +1) / (2*. 5  +  1) ) * (2/ (pi* dvnat ) )  ; 


end 


%%  NUMBER  DENSITY  - 

function  [ndensity]  =  NumDensityK (T) 

%  calculates  number  density  (cm~-3)  given  Temperature  in  Kelvin 
global  kB; 
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ndensity  =  (KVaporPressure (T) )/ (kB*T) *101325* (le-6) ; 
end 


%%  VAPOR  PRESSURE  - 

function  [Pv]  =  KVaporPressure (T) 

%  Provides  K  vapor  pressure  in  atmospheres  given  T  in  Kelvin 

o, 

o 

%  REFERENCE: 

%  C.B.Alcock,  V.P.Itkin,  and  M.K.Horrigan 

%  Vapour  pressure  equations  for  the  metallic  elements:  298-2500K 
%  Canadian  Metallurgical  Quarterly,  1984,  23,  309-313 

o, 

o 

%  "Precise"  equation  for  K  vapour  pressure. 

%  Uses  melting  point  for  K  at  336. 8K  (63.65  degC) 

%  Valid  for  temperatures  from  melting  point  to  600K. 

o, 

o 

%  For  best  precision,  use  temperatures  between  melting  point  and  600K. 

if  (T  >  336.8  &&  T  <  600) 
a  =  8.233; 
b  =  -4693; 
c  =  -1.2403; 

Pv  =  10.  "(a  +  b*(T,'-l)  +  c*loglO  (T) ) ; 
elseif  (T  >  298  &&  T  <  336.8) 
a  =  4.961; 
b  =  -4646; 

Pv  =  10 . ~  (a  +  b* ( T  ~  — 1 ) ) ; 
elseif  (T  <  298  | |  T  >  600) 

fprintf ( ' \nTemperature  is  outside  of  allowed  range. \n\n') 
return 

end 

end 


%%  CROSS  SECTION  - 

function  [cs]  =  CrossSectionK (x, T, PN2 , PHe, PCH4 , PC2H6 , P3He) 

%  Calculates  cross  section  as  function  of  frequency,  x 
global  tau  sO  vo39  vo41  TransM  vHF  numTrans  Iso39  Iso41  c  SFF; 


dvnat  =  1/ (2*pi*tau) ; 
csHF  =  zeros (1, numTrans) ; 


%  natural  linewidth,  a  nu  natural 
%  Cross  section  for  an  individual  HF  transition 


temp  =  0; 

for  n  =  1: numTrans 

Flwr  =  TransM(n,3); 

IsoNum  =  TransM(n,5); 

[gamma, shift , dvD, dvL]  =  Buf ferGasK (T, IsoNum, PN2 , PHe, PCH4 , PC2H6, P3He) ; 
if (IsoNum  ==  Iso39)  %  Set  central  frequency  based  on  isotope  (IsoNum) 

vo  =  vo39; 

elseif (IsoNum  ==  Iso41) 
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vo  =  vo41; 


else 

fprint f ( ' \nERROR  in  CrossSect ionK-Isotope  Number  is  outside 
'allowed  range . \n\n ') ; 
return; 

end 

sO  =  SFF  (n) * (1/tau) *(((c*100)~2)/ (4*pi) ) * (1/ ( (vo+vHF (n) +shift) A2) ) *.  . . 

( (2*TransM(n,2)+l)/(2*.5  +  l)  )  ; 

csHF  =  KNat Abundance (I soNum)  *  BoltzRatioK (IsoNum, Flwr, T)  *  sO  *  ... 

( (x  +  vHF (n)  +  shift  +  vo)  /  (vHF (n)  +  shift  +  vo) )  *  . . . 
gVoigt ( (x  -  vHF (n) -shift ), dvL, dvD) ; 
temp  =  temp  +  csHF; 

end 

cs  =  temp; 
end 

o,  o. 
o  o 

%  Need  Total  Nuclear  Ang.  Momentum 
function  eta  =  KNatAbundance ( IsoNum) 

%  Return  Potassium  natural  abundance  based  on  Total  nuclear  angular  momentum: 
global  Iso39  Iso41; 

%  REFERENCE: 

%  T.G.  Tiecke 
%  Properties  of  Potassium 

%  http :/ /staff . science . uva . nl/-itgtiecke/PotassiumProperties . pdf 

%Natural  Abundance  Ratios: 

Eta39  =  0.932581;  %Nuclear  Spin  =  3/2 
Eta41  =  0.067302;  %Nuclear  Spin  =  3/2 

%Note:  Eta40  =  0.000117;  %Nuclear  Spin  =  4  and  is  ignored  here  b/c  of  low  natural 
%abundance . 

if  (IsoNum  ==  Iso39) 
eta  =  Eta39; 
return ; 

elseif  (IsoNum  ==  Iso41) 
eta  =  Eta41; 
return ; 

else 

fprintf ( ' \nIsotope  Number  is  outside  of  allowed  range  for  K39  or  K41.\n\n'); 
return ; 

end 

end 

%%  BUFFER  GAS  EFFECTS  - 
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function  [gamma, shift ,  dvD, dvL]  =  Buf ferGasK (T, IsoNum, PN2 , PHe, PCH4 , PC2H6, P3He) 
%  Rb  buffer  gas  effects  -  pressure  broadening  and  frequency  shift 

o, 

o 

%  REFERENCES: 

%  Nyunt  Lwin  and  D  G  McCartan 

%  Collision  broadening  of  the  potassium  resonance  lines  by  noble  gases 
%  J.  Phys .  B:  Atom.  Molec.  Phys . ,  Vol .  11,  No.  22.  1978. 

o, 

o 

%  Pitz,  G.  A. 

%  Pressure  broadening  and  shift  of  the  potassium  D1  transition  by  the 
%  noble  gasses  and  N2,  H2,  HD,  D2,  CH4 ,  C2H6,  C38,  and  n-C4H10 
%  JQSRT,  113,  387-395  2012. 

o. 

o 

%  Pitz,  G.  A. 

%  Private  communication  for  D2  line  broadening  and  shift. 

%  2012-MAY-01. 


global  tau  c  kB  Iso39  Iso41  MK39  MK41  vo39  vo41; 


%Calculate  Broadening: 

gamma  =  (18 . 84e6) *  (  (T/400) ~ (1/2) ) *PN2  ...  % 

+  (19 . 84e6) * ( (T/318) " (1/2) ) *PHe  ...  % 

+  (27.78e6) * ( (T/318) ~ (1/2) ) *PCH4  ...  % 

+  (27 . 60e6) *  ( (T/318) " (1/2) )*PC2H6. . .  % 

+  (21 . 14e6) * ( (T/328) " (1/2) )*P3He;  % 

%D 1  Pitz  # s ) 

%  gamma  =  0;  %for  test  purposes  to  check  the 


MHz/Torr 

MHz/Torr 

MHz/Torr 

MHz/Torr 

MHz/Torr 


(Lwin) 

(Pitz  Commun . ) 
(Pitz  Commun . ) 
(Pitz  Commun . ) 
(Extrap.  Guess 


buffergas  shift 


using 


%Calculate  Shift: 

shift  =  (-10 . 14e6) *( (T/400) " (1/2) ) *PN2  ... 
+  (  0 . 52e6) *( (T/318) ~ (1/2) ) *PHe  ... 

+  (  -8.38e6) * ( (T/318) " (1/2) ) *PCH4  ... 
+  (  -8.04e6) *( (T/318) ~ (1/2) )*PC2H6.  .  . 
+  (  7 . 72e6) * ( (T/328) " (1/2) ) *P3He; 

%D 1  Pitz  # s ) 


o. 

o 

MHz/Torr 

(Lwin) 

o, 

o 

MHz/Torr 

(Pitz 

Commun . ) 

o. 

o 

MHz/Torr 

(Pitz 

Commun . ) 

o, 

o 

MHz/Torr 

(Pitz 

Commun . ) 

o_ 

o 

MHz/Torr 

(Extrap.  Guess 

%  Set  central  frequency  and  mass  based  on  isotope 
if  (IsoNum  ==  Iso39) 
vo  =  vo39; 

MK  =  MK39 ; 

elseif ( IsoNum  ==  Iso41) 
vo  =  vo41; 

MK  =  MK4 1 ; 

else 

fprintf ( ' \nERROR  in  Buf ferGasK-Iso  is  outside  range . \n\n ') ; 
return ; 

end 


%Calculate  Doppler  and  Lorentz: 

dvD  =  (vo  +  shift)*((8*kB*T*log(2))/ (MK* (c~2) ) ) " (1/2) ; 
dvL  =  (1/ (2*pi*tau) )  +  gamma; 
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end 


%%  BOLTZMAN  RATIO  - 

function  [BR]  =  Bolt zRat ioK ( IsoNum, F, T) 
global  kB  h  dv39S12  dv41S12  Iso39  Iso41; 

if (IsoNum  ==  Iso39) 
if (F  ==  1) 

BR  -  5 . / (5  +  7*exp (- (h*dv39S12) . / (kB. *T) ) ) ; 
elseif (F  ==  2) 

BR  =  (7*exp (- (h*dv39S12 ) . / (kB . *T) ) ) . / (5  +  7*exp (- (h*dv39S12 ) . / (kB . *T) ) ) ; 
elseif  (F  ^  1  &&  F  ^  2) 

fprintf ( ' \nIMPROPER  INPUTS\n\n’) 

BR  =  0 ; 

end 

elseif ( IsoNum  ==  Iso41) 
if (F  ==  1) 

BR  =  3 . /  ( 3  +  5*exp (- (h*dv41S12)  . / (kB. *T)  )  )  ; 
elseif (F  ==  2) 

BR  =  (5*exp (- (h*dv41S12) . / (kB. *T) ) ) . / (3  +  5*exp (- (h*dv4 1S12 ) . / (kB . *T) ) ) ; 
elseif (F  ^  1  &&  F  ^  2) 

fprintf (' \nIMPROPER  INPUTS\n\n') 

BR  =  0 ; 

end 

else 

fprintf (' \nIMPROPER  INPUTS\n\n') 

BR  =  0; 

end 

end 


%%  VOIGT  LINESHAPE  - 

function  [y]  =  gVoigt (v, dvL, dvD) 

%  A  useful  approximation  for  the  Voigt  convolution  is  given  by  Thompson 

o, 

o 

%  REFERENCE: 

%  William  J.  Thompson 

%  Numerous  neat  algorithms  for  the  Voigt  profile  function. 

%  Computers  in  Physics,  Vol  7,  No.  6,  Nov/Dec  1993. 

o, 

o 

%  gV(u,a)  =  Re [Exp ( z ~ 2 ) erf c ( z ) ] 

%  z  =  a  +  iu 

%  a  =  sqrt(ln(2))  vLorentzian/vDoppler 

%  u  =  2  sqrt(ln(2))  (v-vo) /vDoppler 

%  NOTE:  Not  valid  for  Doppler  free  spectra  where  vDoppler  =  0 
%  NOTE:  erfc(z)  =  1  -  erf(z)  (related  to  the  Faddeeva  function) 

a  =  sqrt (log (2) ) * (dvL/dvD) ; 
u  =  2*sqrt (log (2) ) * (v/dvD)  ; 
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z  =  a  +  li*u; 

y  =  2*sqrt (log (2) /pi) * (1/dvD) *real (faddeeva  (li*z) ) ; 
end 


%  -  gVoigt  subfunctions  - 

function  w  =  f addeeva ( z , N) 

%  FADDEEVA  Faddeeva  function 

%  W  =  FADDEEVA (Z)  is  the  Faddeeva  function,  aka  the  plasma  dispersion 

%  function,  for  each  element  of  Z.  The  Faddeeva  function  is  defined  as: 

o. 

o 

%  w(z)  =  exp(-z~2)  *  erfc(-j*z) 

o, 

o 

%  where  erfc(x)  is  the  complex  complementary  error  function. 

o, 

o 

%  W  =  FADDEEVA (Z, N)  can  be  used  to  explicitly  specify  the  number  of  terms 
%  to  truncate  the  expansion  (see  (13)  in  [1]) .  N  =  16  is  used  as  default. 


Example : 

x  =  linspace (-10, 10, 1001) ;  [X,Y]  =  meshgrid (x, x) ; 

W  =  faddeeva (complex (X, Y) ) ; 
figure ; 

subplot (121)  ;  imagesc (x, x, real (W) ) ;  axis  xy  square;  caxis([-l  1]); 
title (' re ( faddeeva ( z ) ) ' ) ;  xlabel ( ' re ( z ) ' ) ;  ylabel ( ' im ( z ) ' ) ; 
subplot ( 122 ) ;  imagesc (x, x, imag (W) ) ;  axis  xy  square;  caxis([-l  1]); 
title ( ' im ( faddeeva ( z ) )  ' )  ;  xlabel ( ' re ( z )  ' ) ;  ylabel ( ' im ( z )  ' ) ; 

Reference : 

[1]  J.A.C.  Weideman,  "Computation  of  the  Complex  Error  Function,"  SIAM 
J.  Numerical  Analysis,  pp.  1497-1518,  No.  5,  Vol.  31,  Oct.,  1994 
Available  Online:  http://www.jstor.org/stable/2158232 


if  nargin<2,  N  =  [];  end 
if  isempty(N),  N  =  16;  end 

w  =  zeros (size  (z) ) ;  %  initialize  output 


o,  o,  o.  o,  o, 
o  o  o  o  o 

%  for  purely  imaginary-valued  inputs,  use  erf  as  is  if  z  is  real 
idx  =  real(z)==0;  % 

w (idx)  =  exp ( — z (idx) . ~ 2)  . *erfc (imag (z (idx) ) ) ; 

if  all (idx),  return;  end 
idx  =  -iidx; 


o,  o_  o,  o.  o_ 
o  o  o  o  o 

%  for  complex-valued  inputs 

%  make  sure  all  points  are  in  the  upper  half-plane  (positive  imag.  values) 
idxl  =  idx  &  imag(z)<0; 
z (idxl)  =  con j (z (idxl) ) ; 
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M  =  2  *N; 

M2  =  2  *M; 

k  =  (-M+1:1:M-1) 1 ;  %  M2  =  no .  of  sampling  points. 

L  =  sqrt (N/sqrt (2 ) ) ;  %  Optimal  choice  of  L. 

theta  =  k*pi/M; 

t  =  L*tan (theta/2 ) ;  %  Variables  theta  and  t. 
f  -  exp (-t . " 2 )  . *  (L~2+t . ~ 2 ) ; 
f  =  [0;  f] ;  %  Function  to  be  transformed. 

a  =  real ( f ft ( f ft shift ( f) )) /M2 ;  %  Coefficients  of  transform, 
a  =  f lipud (a (2 : N+l ) ) ;  %  Reorder  coefficients. 

Z  =  (L+li*z (idx) )  . / (L-li*z (idx)  )  ; 
p  =  polyval (a, Z ) ;  %  Polynomial  evaluation. 

w(idx)  =  2*p . / (L-li*z (idx) )  . ~2  +  (1/sqrt (pi) )  . / (L-li*z (idx) )  ; 

%  convert  the  upper  half-plane  results  to  the  lower  half-plane 
w(idxl)  =  con j (2*exp (-z (idxl ) . "2 )  -  w(idxl)); 

end 


Evaluate  w ( z ) . 
if  necessary 
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Use  of  K  Di  and  D2  Absorption  Lines 


The  following  code  sample  illustrates  use  of  Matlab®  functions  fnAbsCoef  fKDl  .m, 
fnAbsCoef  fKD2  .m  creating  and  output  for  both  the  D\  and  D2  K  lines.  The  code 
below  was  used  to  create  the  spectral  profiles  for  the  pressure-broadened  DPAL  gain 
cells  used  to  calculate  atmospheric  molecular  absorption. 


%  Adapted  -  2012  -  CAR  -  From  PMonte  Anderson's  original  code 

%  m-File  tha  computes  Absorption  Coefficient  Dl  and  D2  using  fnAbsCoef fCsDl 
%  and  fnAbsCoef fCsD2 . 

%% 

%  Cleanup  Environment: 

close  all 
clear  all 
clc 


%  General  Settings: 

%Start  Clock: 

tic; 

%Vary  One  Of  The  Input  Parameters: 

numPoints  =  10000; 

VarPMat  =  7 60 . * [ 1 , 5 , 10 , 20 ] ; 

numRuns  =  length (VarPMat ) ; 

%Pre-create  Each  Outout  Matrix: 

absCoef fDIMat  =  zeros (numPoints , numRuns ) ; 

absCoef fD2Mat  =  zeros (numPoints , numRuns ) ; 

%% 

%Run  Sim  For  Each: 

for (m  =  1: numRuns) 

%%  calcuations 
c  =  2 . 99792458e8;  %m/s 

voDl  =  389286058.716;  %MHz 
voD2  =  391016170.03;  %MHz 
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%Detune  Freq  (GHz) : 

nu  =  linspace (-5000e9, 5000e9, numPoints) ; 

T  =  60  +  273.15;  %T  in  K  -  What  should  be  used  here 
PN2  =  0;  %0  (Torr) 

PHe  =  VarPMat (m) ;  %200  %Typical  DPALs  Helium  pressure  (Torr) 

PCH4  =  0;  %400  %Typical  DPALs  methane  pressure  (Torr) 

PC2H6  =  100;  %100  %Typical  DPALs  ethane  pressure  (Torr) 

P3He  =  0; 
cell_length  =  5; 

absCoeffDl  =  zeros (length (nu) , 1) ; 
absCoeffD2  =  zeros ( length (nu) , 1 ) ; 
lambdaDl  =  zeros (length (nu) f 1) ; 
lambdaD2  =  zeros ( length (nu)  ,  1 )  ; 
nuDl  =  zeros  (length (nu) , 1) ; 
nuD2  =  zeros (length (nu) , 1) ; 

for  i=l : length (nu) 

absCoeffDl (i)  -  f nAbsCoef fKDl (nu(i) , T, PN2, PHe, PCH4, PC2H6, P3He) ; 

absCoef fD2 (i)  =  f nAbsCoef fKD 2 (nu (i) , T, PN2 , PHe, PCH4 , PC2H6 , P3He ) ; 

lambdaDl (i)  =  c ./ (nu (i) * (le-9) +voDl* (le-3) ) ; 

lambdaD2 (i)  =  c ./ (nu (i) * (le-9) +voD2*  (le-3) ) ; 

nuDl (i)  =  (nu (i) tvoDl* (le6) )  . /  (c*100) ; 

nuD2 (i)  =  (nu  (i) +voD2* (le6) )  . / (c*100) ; 

end 

%Store  the  Data: 

absCoef fDIMat (:  ,  m)  =  absCoeffDl; 
absCoef fD2Mat (: , m)  =  absCoeffD2; 

absCoef fDIMatNorm (:, m)  =  absCoef fDl . /max (absCoef fDl ) ; 
absCoef fD2MatNorm (:, m)  =  absCoef fD2 . /max (absCoef fD2 )  ; 

%End  Of  Calc: 
toe 

end 

%%  — l — l — l — l — l — l — l — l — l — V 

%  Plotting: 

close  all; 

opts  =  [ ' f ont_size , 1 , f ont_scale , 1 . 15 ’ ] ; 
fSz  =  1; 
f Sc  =  1.15; 

figurePos  =  [10  60  1000* (1+sqrt (5) ) /2  1000]; 

%Plot  Each  Line  in  Wavenumber: 
figured- 

plot  (nuDl, absCoef fDIMat ) 

xlabel ( ’ Wavenumber,  $cmA {- 1 } $  1 ,  ' Interpreter '  ,  'latex ')  ; 
ylabel('Abs  Coeff  $ (cm" {-1} )$',' Interpreter 'latex'); 
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title(['$D_l$'],  ' Interpreter ' , 'latex ' ) 

legend (num2str (VarPMat ' ) ) 

set (geo,  ' Interpreter ' ,  ' latex ' ) ; 

nPlot (gef , ' f ont_size ' , fSz, 1 font_scale ' , fSc) ; drawnow; 
set (gef,  'Position' , figurePos) ; 


%%  Plot  Each  Line  in  Relative  Freq: 
figure; 

temp  =  linspace ( 0 , 0 . 9 , numRuns ) ; 
cMat  =  [temp; temp; temp] ; 
for  (m  =  1: numRuns) 

plot  (nu*le-9, absCoef fDIMat ( : , m) ,  'Color ' , cMat ( : , m)  '  ) 
hold  on; 

end 

xlabel (' Detuning  Freq.  (GHz)'); 
ylabel('Abs.  Coeff.  (crrT{-l})  ' )  ; 

return; 
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Frequency  Calibration  of  Tunable  Laser  Data  Using  an  Etelon 

The  following  code  stitches  together  spectrally  adjacent  laser  scans  of  the  New- 
Focus  tunable  diode  lasers  and  calibrates  the  x-axis  of  the  final  data  using  an  etelon 
in  the  absence  of  a  high  resolution  wavemeter.  The  algorithm  first  linearizes  each 
individual  scan  using  the  etelon  response  by  finding  the  peaks  given  the  approxi¬ 
mate  etelon  free  spectral  range.  Then,  detector  response  for  the  laser  power  reference 
signal  is  used  for  adjacent  scans  to  find  where  the  two  laser  scans  should  be  joined 
together.  Additional  signals  can  be  used  to  improve  the  estimate,  however  the  etleon 
and  laser  power  reference  signal  typically  work  best.  Natural  eteloning  in  the  system 
also  improves  the  process  of  joining  adjacent  free  spectral  ranges. 


function  [completeDataOut ]  =  mergeFSREtelon (fullData) ; 

%%  — ^ — l — l — l — l — l — l — l — l — ^ — l — ^ — l — l — l — l — l — l — l — l — I — l — ^ — ( — l — l — l — l — l — l — l — ^ ^ — l — l — l- 
%  Created  -  2012-03-22  -  CAR  -  Chris  Rice  -  crice@afit.edu 
%  Updated  -  2012-03-22  -  CAR  -  Adapted  from  merge  FSR. 

%  m-File  that  takes  the  full  data  structure  and  creates  a  full  data  output. 

%%  Init: 

%Verbose  Mode: 
verbose  =  1; 

%Input  Structure  Should  Look  Like: 
szln  =  size (fullData) ; 

%Number  Of  Points  In  Each  Collection: 
numPts  =  szln(l); 

%Input  Vector  is  5  columns: 

%  Setup: 

%  SRS510  I - >DAQ  AnlnO 

%  SRS510  10 - >DAQ  Anlnl 

%  SRS510  Ref - >DAQ  Anln2 

%  SRS530  Etelon->DAQ  Anln3 

%  NF  Lambda - >DAQ  Anln4 

numChan  =  szln  (2); 

%Number  of  Scans  In  Collection: 
numScans  =  szln (3); 
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%%  Linearize  The  Data: 
tempSz  =  size  ( fullData) ; 

%Pre-create  the  output: 

fullDataLinear  =  zeros (tempSz (1) , tempSz (2) +1, tempSz (3) ) ; 

%For  each  scan,  use  FSR  to  linearize  data: 
for (m  =  l:numScans) 

%Find  Freq.  Axis: 

[ f reqOut , f reqln]  =  linearizeEtelonScan ( fullData (:, 4 , m) ) ; 

%Use  Old  and  New  Freq.  Matricies  to  re-scale  data: 
for  (n  =  l:numChan) 

fullDataLinear ( : , n, m)  =  interpl ( f reqln, fullData ( : , n, m) , f reqOut ) ; 

end 

fullDataLinear (:, 7 , m)  =  freqOut; 

end 

%%  Heavy  Lifting: 
for (m  =  2 :numScans) ; 

if (verbose) ; disp ([' Joining  Scan  #  =  ' , num2str (m) ] ) ; end 
%Pull  Out  Data: 

leftFSR  =  fullDataLinear (:,:, m-1 ) ; 
rightFSR  =  fullDataLinear (:,:, m) ; 

%Pass  the  left  FSR  and  the  right  FSR  Scans: 

[LFSRInd (m) , RFSRInd (m) ]  =  mergeTwoFSREtelon (leftFSR, rightFSR) ; 

end 

%%  Merge  The  Data  Into  One  Marix: 
for (m  =  1 :numScans) 
if (m^l) 

completeDataOut  =  [ completeDataOut ; fullDataLinear (RFSRInd (m) :end, :,m)]; 

else 

completeDataOut  =  fullDataLinear ( 1 : LFSRInd (m+1 ),:, m) ; 

end 

end 

%%  Check  The  Relative  Freq.  Channel,  Make  Continuous: 

CDODiff  =  diff (completeDataOut (:,  7 ))  ; 

%Find  All  Positive  Movement: 
poslnd  =  f ind (CDODif f>0) ; 
neglnd  =  find (CDODif f<0 )  ; 

%Find  Where  Points  Lie: 

CDOMean  =  mean (CDODif f (poslnd)  )  ; 

CDOSTD  =  std (CDODiff (poslnd)  )  ; 

%Cycle  Through  Scans : 

currMaxFreq  =  completeDataOut (neglnd ( 1 ),  7 )  ; 
freqNew  =  zeros  ( size (completeDataOut (:,  7 )))  ; 
f reqNew ( 1 : neglnd ( 1 ) )  =  completeDataOut ( 1 : neglnd ( 1 ),  7 )  ; 
for (m  =  2 : length (neglnd) ) 

freqNew (neglnd (m-1 ) +1 : neglnd (m) )  =  completeDataOut (neglnd (m-1 ) +1 : neglnd (m) , . . . 

7 ) +currMaxFreq-min (completeDataOut (neglnd (m-1 ) +1 : neglnd (m) , 7 ) ) ; 

%Current  Max  Freq: 
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currMaxFreq  =  f reqNew (neglnd (m) ) ; 

end 

%Finish  Off  Last  Scan: 

f reqNew (neglnd (m) +1 : end)  =  completeDataOut (neglnd (m) +1 : end, 7 ) +currMaxFreq- . . . 
min (completeDataOut (neglnd (m) +1 : end,  7 ) ) ; 

%Re-write  Old  Matrix: 
completeDataOut (:, 7 )  =  f reqNew; 


%%  Cleanup: 
%None . 
return; 
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Bristol  Wavemeter  control  by  Matlab® 


The  following  code  demonstrates  the  use  of  a  Bristol  621 A  wavemeter  with  Mat- 
lab®.  In  order  to  use  the  following  code,  the  Bristol  install  package  must  have  been 
completed,  and  the  hies  CLDevIfaee.h,  CLDevIface.dll,  and  CLDevIface.lib  must  be 
in  the  MatLab  active  path.  The  Bristol  wavemeter  has  a  relatively  slow  and  inaccu¬ 
rate  response  compared  the  the  High  Finesse  WSU2,  but  its  use  during  the  course  of 
an  experiment  is  still  preferred  to  no  spectral  reference  and  is  an  excellent  diagnostic 
tool. 


function  [libName, BristolDevHandle]  =  initBristol ( comNumber ) ; 

%% 

%  Created  -  2009-02-05  -  CAR  -  Chris  Rice  -  crice@afit.edu 
%  Updated  -  2009-02-12  -  CAR  -  Adapted  from  test  code  of  runWavemeter 

%  Loads  Bristol  SI  Labs  CP2101  USB-2-UART  Library  for  access  to  Bristol 
%  621A  from  MatLab 

%  NOTE  -  Currently  closure  of  the  Bristol  dynamic-linked  library  causes 
%  program  termination  -  do  not  close  library  until  this  is  fixed  by  Bristol 

%% 

%  Attempt  to  Load  Library: 

libName  =  ' BristolLibrary ' ; 

%Check  For  Load: 

libLoaded  =  libisloaded (libName) 
if (libLoaded) 

%If  Loaded,  Unload  and  re-load: 

% (Some  versions  of  Britol  Lib.  Fail  During  MatLab  Unload  Calls) : 
unloadlibrary (libName) 

end 

%Loads  Library: 

loadlibrary ( ' CLDevI f ace .dll',  ' cldevdll .h',  'alias', libName )  ; 

%Display  Functions  In  Library: 
libfunctions (libName) ; 

%View  Functions  of  Library: 
libfunctionsview (libName) 
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%View  Version  of  Library: 

call lib ( libName, ' CLGetDllVersion ' ) ; 

%% 

%  Open  The  Device: 

BristolDevHandle  =  calllib (libName, 1 CLOpenUSBSerialDevice f , comNumber ) ; 

%% 

%  Test  A  Few  Calls  To  Library: 

%— H — I — l — l — l — l — I — l — l — l — l — l — l — l — h 

Blambda  =  calllib (libName, 1 CLGetLambdaReading ' , BristolDevHandle) 

Bpower  =  calllib (libName, ' CLGetPowerReading ', BristolDevHandle) 

%% 

%  Attempt  to  Cleanup: 

%— H — l — l — l — l — l — l — l — l — I — l — l — h 

%Close  the  Device: 

% (Some  versions  of  Britol  Lib.  Fail  During  MatLab  Unload  Calls) : 
test  =  calllib ( libName, ' CLCloseDevice ', BristolDevHandle ) 

%Unload  the  Library: 

% (Some  versions  of  Britol  Lib.  Fail  During  MatLab  Unload  Calls)  : 
unloadlibrary (libName) 

return 
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High  Finesse  Wavemeter  control  by  Matlab® 


The  following  code  demonstrates  the  use  of  a  High  Finesse  WSU2  wavemeter  with 
Matlab®.  In  order  to  use  the  following  code,  the  High  Finesse  install  package  must 
have  been  completed,  and  the  hie  wlmData.dll  must  be  in  the  MatLab  active  path  or 
Windows  system  32  library  folder  along  with  its  associated  header  hie.  The  routine 
also  demonstrates  how  public  functions  of  the  shared  library  can  be  investigated, 
leading  future  programmers  towards  other  development  options. 


function  [libName]  =  initHFWM; 

%% 

%  Created  -  2010-05-01  -  CAR  -  Chris  Rice  -  crice@afit.edu 

%  Updated  -  2010-05-04  -  CAR  -  Changed  Source  Origin  to  \{OS}\system32 

%  Updated  -  2010-05-06  -  CAR  -  Altered  function.  Integrated  into  system 

%  Updated  -  2011-05-23  -  CAR  -  Added  more  comments  for  testing  with  Stefan 

%  Burgardt,  inserted  more  'save'  commands  for  testing. 

%  Loads  HighFinesse  WS-U2  Access  dll  for  control  of  wavemeter  from  MatLab 
%  Proprietary  information!  Copyright  Chris  Rice.  Proprietary  information! 

%% 

%  Attempt  to  Load  Library: 

libName  =  ' WSULibrary ' ; 

%Test  to  see  if  'save'  command  still  works: 
save ( ' Testl .mat ' )  ; 

%Check  For  Load: 

libLoaded  =  libisloaded (libName) ; 
if (libLoaded) 

%If  Loaded,  Unload  and  re-load:  (Should  We  Reload  or  let  it  be?) 
disp (' HighFinesse  Library  Already  Loaded  -  Unloading...'); 
unloadlibrary (libName)  ; 

disp (' HighFinesse  Library  Unload  Successfull .  '  )  ; 

end 

%Test  to  see  if  'save'  command  still  works: 
save ( ' Test 2 . mat ' ) 

%Loads  Library 
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% IMPORTANT  -  ONCE  THIS  LIBRARY  IS  OPEN,  THE  SHARED  MEMORY  IS  USED  BY  THIS  % 
%LIBRARY  ONLY,  OTHER  COPIES  WILL  NOT  WORK  WHILE  THIS  IS  OPEN: 

loadlibrary ( ' C : \windows\ system32\wlmData . dll ' , ' C : \windows\ system32\wlmData . h 1 , . . . 
' alias ' , libName) ; 

disp ( ' HighFinesse  Library  Load  Successfull . ' ) ; 

%Test  to  see  if  'save'  command  still  works: 
save ( ' Test 3 .mat ' ) 

%Display  Functions  In  Library  (Not  Needed  as  a  function) : 

%  libfunctions (libName, 1 -full ' ) ; 

%GUI  View  Functions  and  Version  of  Library  (Not  Needed  as  a  function) : 

%  libfunct ionsview (libName) 

%  Open  The  Device: 

%  Instantiate  -  Check  for  Server,  If  not  start,  start  it:  (Is  This  Required?) 
WLMServerOn  =  0; 

% Inst antiate (ReasonForCall, Mode, PI ,  P2 ) 

%RFC  [-l=Check  For  Server  |  0=Reset  |  l=Insall  or  Remove  Callback  for  Export 
%  |  2=Copy  Pattern  |  3=InstControlWLM  (Obsolete)  |  4=Delay  (Obsolete)  | 

%  5=Priority  (Obsolete) ] 

%Mode  [ 0=Not if ylnstallCallback  |  l=Not ifyRemoveCallback  | 

%  2=Not ifylnstallWaitEvent  |  3=Not ifyRemoveWaitEvent  | 

%  4=Not ifylnstallCallbackEx  |  5=Notif ylnstallWaitEventEx] 

%P1  [Callback  Address  for  install  of  callback  or  timeout  for  wait-event] 

%P2  [Priority  of  Callback  (0=Standard) ] 

WLMServerOn  =  calllib (libName,  1  Instantiate ', -1, 0, 0, 0) ; 

%Test  to  see  if  'save'  command  still  works: 
save ( 1  Test 4 . mat 1 ) 

if (WLMServerOn==0 ) 

%Test  to  see  if  'save1  command  still  works: 
save ( ' Test5 .mat ' ) 

%Start  The  WLM  Server: 

%ControlWLM (Action, App, Ver ) 

%Act  [ l=WLMShow  |  2=WLMHide  |  3=WLMExit  |  10=WLMWait  | 

%20=WLMStartSilent  |  40=WLMSilent ] 

%App  [Pointer  to  String  of  Application  to  Start  ( 0=NoAppToStart ) ] 

%Ver  [0=Last  Used  Server] 

WLMServerOn  =  calllib (libName,  ' ControlWLM'  ,  1,  0, 0) ; 

end 

%Test  to  see  if  'save'  command  still  works: 
save ( ' Test 6 .mat ' ) 

%Display  Wavemeter  Version: 
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WLMVersion  =  calllib ( libName, ' GetWLMVersion 1 , 1 ) ; 
disp ( ['Using  WSU-2  Version  ' , num2str (WLMVersion) ]) ; 

%Set  the  Result  Mode  UNIT-> [ 0=WLVac,  l=WLAir,  2=Freq,  3=WN,  4=PhotEnergy ] : 
theError  =  calllib ( libName, ' SetResultMode 1 ,  3 ) ; 

%Test  to  see  if  'save'  command  still  works: 
save ( 1 Test 7 .mat ' ) 

%Set  Exposure  Mode  [O=0ff,  l=Auto] : 

theError  =  calllib ( libName, ' SetExposureMode f , 1 ) ; 

%Test  to  see  if  'save'  command  still  works: 
save ( 1  Test 8 . mat 1 ) 

%Start  Measurement : 
if (WLMServerOn==l ) 

%Start  a  measurement  operation  [ 0x00=Stop, 0x01=Ad justment , 

%  0x02=Measurement , 0x04=Record, 0x08=Replay]  ; 

WLMOperation  =  calllib (libName,  f Operation 2 ) ; 

end 

%Take  A  Measurement : 

tempWave  =  calllib ( libName, ’ GetWavelengthNum ' , 1 , 0 ) ; 

% — l — l — l — l — l — l — I — l — l — l — l — l — I — l — l — I — l — I — \ — , — l — l — I — I — l — l — l — l — l — l — , — l — l — I — l — l — (- 
%A11  code  in  this  document  is  Proprietary  information!  Copyright  Chris  Rice 
% — I — l — I — l — l — l — l — l — I — I — I — I — l — I — l — l — l — l — I — I — I — I — l — l — l — l — l — l — l — I — I — I — I — \ — l — I — h 
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High  Finesse  Wavemeter  Calibration  and  Use 


The  following  code  calibrates  a  High  Finesse  WSU2  wavemeter  using  a  frequency 
stabilized  HeNe  laser  using  Matlab®.  In  order  to  use  the  following  code,  the  High 
Finesse  install  package  must  have  been  completed  and  the  file  wlmData.dll  must 
be  in  the  Matlab®  active  path  or  Windows  system  32  library  folder  along  with 
its  associated  header  hie.  This  only  works  with  High  Finesse  wavemeters  with  the 
optional  calibration  laser  port  located  on  the  device  -  the  code  would  need  to  be  altered 
to  use  the  external  fiber  switch,  a  popular  option  on  High  Finesse  wavemeters. 


function  [theError]  =  calHFWM (HFWMLib) ; 

%Chris  Rice 

%Created  -  2010-11-10  -  CAR  -  crice@afit.edu 

%Updated  -  2010-11-10  -  CAR  -  Used  long  string  of  inputs  to  make  this  cal 
%and  cleaned  up  code  outside  of  this  function. 

%Switch  To  Rear  Port  (Mode-> [l=Serial, 2=Separated  Low, 3=Seperated  Specific] 

%, Port-> [Front  =  l , Back=2 ] , Channel, Not  Used)  : 

theErrorl  =  calllib (HFWMLib,  ' SetActiveChannel ’ , 3, 2, 5, 0) ; 

%Stop  a  measurement  operation  [0x00=Stop, 0x01=Ad justment, 0x02=Measurement , 
%0x04=Record, 0x08=Replay]  ; 

WLMOperation  =  calllib (HFWMLib,  ' Operation ',  0 )  ; 

%Calibrate  the  Wavemeter  (Type-> [0=HeNe] , UNIT-> [0=WLVac,  l=WLAir,  2=Freq, 
%3=WN,  4=PhotEnergy] , Value, Channel) : 

theError2  =  calllib (HFWMLib, ' Calibration ', 0, 0, 632 . 9910760, 5) ; 

%Re-Start  a  measurement  operation  [0x00=Stop, 0x01=Adjustment, 
%0x02=Measurement , 0x04=Record, 0x08=Replay]  ; 

WLMOperation  =  calllib (HFWMLib, ' Operation ', 2 ) ; 

%Switch  To  Front  Port  (Mode-> [ l=Serial , 2=Seperated  Low, 3=Seperated 
%Specific] , Port-> [Front  =  l , Back=2 ],  Channel ,  Not  Used)  : 
theError3  =  calllib (HFWMLib,  ' SetActiveChannel ', 3, 1, 1, 0) ; 

%Combine  Errors : 

theError  =  theErrorl | theError2 | theError3 ; 

This  Matlab® function  can  be  called  at  any  time  the  High  Finesse  WSU2  is  not 
being  calibrated  to  retrieve  the  current  frequency  of  the  laser  being  investigated. 
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function  [output]  =  getHFWM (HFWMLib, cmd) ; 

% — I — l — I — I — I — I — l — l — l — I — I — I — l — I — I — I — I — l — I — I — I — I — l — l — I — I — I — I — l — l — I — I — I — l — I — I — f 
%Chris  Rice 

%Created  -  2010-11-11  -  CAR  -  crice@afit.edu 

%Updated  -  2010-11-11  -  CAR  -  Used  to  find  various  params  from  WM 

%  — H - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - h 

%-,WL'  -  returns  an  average  for  vaccum  wavelength  in  nm 

%  — H - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - h 

if (strcmp (cmd, 1 WL ' ) ) 

%Get  An  Average  Value,  check  each  one  for  bad  values: 
for (m  =  1:10) 

theFreqs (m)  =  calllib (HFWMLib, ' GetFrequencyNum ', 1 , 0 ) ; 
if  (theFreqs (m) <0) 
sampTry  =  0; 
while (theFreqs (m) <0) 
pause (0.001) ; 

theFreqs (m)  =  calllib (HFWMLib, 1 GetFrequencyNum ', 1 , 0 ) ; 
sampTry=sampTry+l ; 
if (sampTry>5000) 

error ( 1  Cannot  Get  Good  Sample  from  Wavemeter ' ) ; 
break 

end 

end 

end 

end 

output  =  299792 . 458/mean (theFreqs)  ; 

end 

%Return  Output : 
return 
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HELEEOS  Batch  Runs 


The  following  code  demonstrates  the  routine  used  to  run  HELEEOS  while  varying 
the  scenario  (pre-created  with  the  HELEEOS  GUI),  laser  type,  and  if  a  DPAL,  gain 
cell  pressure.  The  routine  also  saves  the  data  for  future  analysis,  and  the  plotting 
routines  plot  and  sort  the  associated  bar  graphs  so  that  the  bar  graph  order  is  plotted 
with  the  shorted  bars  in  the  front  and  the  tallest  bars  in  the  back. 


%%  +-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 

%  Adapted  -  2012-05-22  -  CAR  -  Chris  Rice  -  crice@afit.edu  -  From 
%  sample_ExPERT_script . m 

%  Updated  -  2012-05-25  -  CAR  -  From  PostProcDataNLinFit . m 
%  Removed  all  MatLab  nLinFit  calls,  will  use  PeakFit  for  fitting 

%  m-File  that  runs  HELEEOS  through  a  variety  of  simulations: 

%%  -+-+-+-+-+-+-+-+-+-+-+-+-+-+ 

%  Setup  Initial  Variables: 

tic; 

global  inputs  outputs 

%Add  path  for  HELEEOS  Scripts: 

%change  to  your  path  settings: 

addpath (genpath ( ' D : \School\Doctorate\HELEEOS\HELEEOS ' ) ) 

%Load  HELEEOS  Pre-Created  Profiles: 

load ( ' D : \ School \Doct or at e\ HELEEOS \HELEEOS\ GUI s\ HELEEOS \mainHeleeos 3 Library . mat ' ) 

%Find  Out  How  Many  Scinarios  There  Are  -  Pull  Out  Settings: 

[numGeometries,^]  =  size ( savedStructures )  ; 
for (m  =  1 : numGeometries ) 

runlnputsDescfm}  =  savedStructuresfm, 1} ; 
runlnputsfm}  =  savedStructuresfm, 2} ; 

end 

%Pick  Out  Laser  Types: 

numLasers  =  4 ;  %K  DPAL,  Rb  DPAL,  Cs  DPAL,  COIL 
numLasersLabels  =  { ' K ' , ' Rb ' , ' Cs ' , ' COIL ' } ; 
useLBLRTMs  =  [1  1  1  1]; 

laserWavelengths  =  [0.7701  0.7949  0.8945  1 . 31525] . *le-006; 
laserTypes  =  [1234]; 
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%Number  Of  Pressures: 
pressures  =  [760  760*10  760*20]; 
numPressures  =  length (pressures ) ; 

%Clear  Out  Saved  Structures  -  No  Longer  Needed: 
clear  savedStructures 

%%  — l — l — l — l — l — l — l — l — l — l — l — l — l — l- 
%  Heavy  Lifting: 

%Initalize  Counters: 
counts  =  0; 

%Step  through  each  scenario  in  HELEEOS  and  calculate  irradiance: 
for (m  =  1 : numGeometries ) 
for (n  =  l:numLasers) 

for (mm  =  1 : numPressures ) 

%Pick  An  Input: 
inputsStruct  =  runlnput s{m} ; 

%Check  Time  of  Calculation: 
tinit  =  toe; 

%Re-structure  Inputs  For  HELEEOS: 

t ransf ormedlnput s  =  GuiStruct2DataStruct ( input sStruct ) ; 

%Need  Assignments  For  The  Following  Variables: 
t ransf ormedlnput s . DPAL . useLBLRTM  =  1; 

%Switch  For  Laser  Wavelength: 

transf ormedlnput s . DPAL . laserWavelength  =  laserWavelengths (n) ; 
transf ormedlnput s . Laser . wavelength  =  laserWavelengths (n) ; 

%Switch  For  Laser  Type: 

transf ormedlnput s . DPAL . laserType  =  laserTypes (n) ; 

%Laser  Pressure  (Torr) : 

transf ormedlnput s . DPAL . laserPressure  =  pressures (mm) ; 
transf ormedlnput s . DPAL . runDesc  =  runlnput sDescjm} ; 

%Execute  calculation: 

outputs  =  heleeosCalcTotalScenario (transf ormedlnputs , false) ; 

%Store  All  Inputs  And  Outputs: 

Alllnputsjm, n, mm}  =  transf ormedlnput s ; 

A110utputs{m, n, mm}=  outputs; 

%Store  Outputs  (Peak  Irradiance,  All  Effects) : 
out{m,n,mm}  =  output s  .  avg .  irrAHEf  feet sPeak; 

%Update  Cycle: 
counts  =  counts  +  1; 

%Update  Display: 
avgTime  =  toe/  (counts) ; 

t imeRemaining  =  (numGeometries*numLasers*numPressures-counts ) * . . . 
avgTime ; 

disp(['Done  with  Site  ' , num2str (counts ),'/', num2str (numGeometries* .. . 
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numLasers*numPressures) ] ) ; 

disp(['Time  Remaining  ' ,  num2str (t imeRemaining) , '  Avg  Time: 

num2str (avgTime) ] ) ; 
pause ( . 01 ) ; 

end 

end 

end 

%% 

%  Save  Those  Outputs: 

tempStr  =  [ ' HELEEOSRuns- ' , datestr (now) , ' . mat ' ] ; 

tempStr  =  strrep (tempStr,  '  :  ' ,  ' -  1 ) ; 

save (tempStr, 'AllOutputs ' , 'Allinputs ' ) ; 

tempStr  =  [' HELEEOSRunsAll- ', datestr (now) ,  '. mat '] ; 

tempStr  =  strrep (tempStr, 1 ; 

save (tempStr) ; 

%% 

%  Plotting  Molec  Abs : 

%Close  Exiting  Figures: 
close  all; 

%Make  a  Figure  For  Each  Scenario: 

figHands  =  []; 

for (m  =  1 : numGeometries ) 

figHands (end+1 )  =  f igure ( ' Name ' , [ runlnput sDescjm} , ' -IntAbs 1 ] ) ; 
for (n  =  l:numLasers) 

for (mm  =  1 : numPressures ) 

%Collect  Settings: 

templn  =  Allinput s{m, n, (mm)}; 

tempOut  =  AllOutputsjm, n, (mm)}; 

%Plot  Results: 

varToPlot  =  ' tempOut . Atmosphere . molecAbs ' ; 
if  (n==4 ) 

faceColorVal  =  [1  1  1]; 

else 

faceColorVal  =  [mm/numPressures  0  1 -mm/ numPressures ] ; 

end 

tempStr  =  [ 'bar (n, sum ( ' , varToPlot , ' ) , ' ' FaceColor ' ' , faceColorVal) ;']; 
tempH (mm)  =  eval(tempStr); 

tempVal (mm)  =  eval ([' sum (', varToPlot ,');']) ; 
tempF  =  get (tempH (mm) ,  1  Child ') ; 
set (tempF, ' FaceAlpha ' , 1 ) ; 
hold  on; 

end 

%Stack  The  Three  Objects,  Highest  In  Back: 

[-i,  tempIndMax]  =  max  (tempVal )  ; 

[-i,  tempIndMin]  =  min  (tempVal )  ; 
uistack (tempH (tempIndMax) , 'bottom' ) ; 
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uistack (tempH (tempIndMin) ,  '  top ’ ) ; 

end 

set (gca, ' XTick ' , [ 1 : numLasers ] ) ; 
set (gca, 1 XTickLabel ' , numLasersLabels ) ; 
ylabel (varToPlot) ; 
title  (runInputsDesc{m} )  ; 

%Make  Figures  Nice: 


end 


%  Plotting  Strehl: 

%Close  Exiting  Figures: 
close  all; 

%Make  a  Figure  For  Each  Scenario: 
varToPlot  =  'beamStrehl'; 
figHands  =  []; 
for (m  =  1 : numGeometries ) 

figHands (end+1 )  =  figure ( ' Name  f ,  [ run Input sDescjm} ,  ' - ' ,  ' BeamStrehl ' ] ) ; 
for  (n  =  l:numLasers) 

for (mm  =  1 : numPressures ) 

%Collect  Settings: 

templn  =  Allinput s{m, n, (mm)}; 

tempOut  =  AllOutputsjm, n, (mm)}; 

%Plot  Results: 

irrDiff  =  tempOut . avg . irrDif fPeak; 
irrDif fBQTurb  =  tempOut . avg . irrDif fBQTurbPeak; 
irrDif fBQTurbTB  =  tempOut . avg . irrDif fBQTurbTBPeak; 
irrAHEff  =  tempOut .  avg .  irrAHEf  feet sPeak; 
beamStrehl  =  irrAHEff /irrDif  f; 
if  (n==4 ) 

faceColorVal  =  [1  1  1]; 

else 

faceColorVal  =  [mm/numPressures  0  1 -mm/ numPressures ] ; 

end 

tempStr  =  [ 1  bar (n, sum ( 1 , varToPlot,  ' ) ,  '  ' FaceColor 1  ' , faceColorVal) ;  ' ] ; 
tempH (mm)  =  eval(tempStr); 

tempVal (mm)  =  eval ([' sum (’, varToPlot ,');']) ; 
tempF  =  get (tempH (mm) ,' Child ')  ; 
set (tempF, ' FaceAlpha 1 , 1 ) ; 
hold  on; 

end 

%Stack  The  Three  Objects: 

[-i,  tempIndMax]  =  max  (tempVal )  ; 

[-i,  tempIndMin]  =  min  (tempVal )  ; 
uistack (tempH (tempIndMax) , 'bottom' ) ; 
uistack (tempH (tempIndMin) , ' top ' ) ; 
disp  (num2str (tempVal) )  ; 

tempValues (m, n)  =  ( 1-min (tempVal ) /max (tempVal )) ; 
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disp(['LT  =  ' , num2str (n) ,  ' %  Diff  =  ' , num2str ( 1-min (tempVal ) / . . . 
max (tempVal) ) ] ) ; 

end 

set (gca, ' XTick ' ,  [ 1 : numLasers ] ) ; 
set (gca, ' XTickLabel ' , numLasersLabels ) ; 
ylabel (varToPlot) ; 
title  (runInputsDesc{m} )  ; 

end 

%% 

%  Plotting  Dwell  Time: 

%Close  Exiting  Figures: 
close  all; 

%Make  a  Figure  For  Each  Scenario: 
varToPlot  =  ' tempOut . dwellTime ' ; 
figHands  =  []; 
for (m  =  1 : numGeometries ) 

figHands (end+1 )  =  figure ( ' Name  1 ,  [ run Input sDescjm} ,  1  - ' ,  1 DwellTime ' ] ) ; 
for (n  =  l:numLasers) 

for (mm  =  1 : numPressures ) 

%Collect  Settings: 

templn  =  Allinput s{m, n, (mm)}; 

tempOut  =  AllOutputsjm, n, (mm)}; 

%Plot  Results: 
if  (n==4 ) 

faceColorVal  =  [1  1  1]; 

else 

faceColorVal  =  [mm/numPressures  0  1 -mm/ numPressures ] ; 

end 

tempStr  =  [ 'bar (n, sum ( ' , varToPlot, ' ) , ' ' FaceColor 1 ' , faceColorVal) ;']; 
tempH (mm)  =  eval(tempStr); 

tempVal (mm)  =  eval ([' sum (’, varToPlot ,');']) ; 
tempF  =  get (tempH (mm) ,' Child ')  ; 
set (tempF, 1 FaceAlpha 1 , 1 ) ; 
hold  on; 

end 

%Stack  The  Three  Objects,  Highest  In  Back: 

[-i,  tempIndMax]  =  max  (tempVal )  ; 

[-i,  tempIndMin]  =  min  (tempVal )  ; 
uistack (tempH (tempIndMax) , 'bottom' ) ; 
uistack (tempH (tempIndMin) , ' top ' ) ; 

end 

set (gca, 'XTick ' , [ 1 : numLasers ] ) ; 
set (gca, ' XTickLabel ' , numLasersLabels ) ; 
ylabel (varToPlot) ; 
title  (runlnputsDescjm} )  ; 

%Make  Figures  Nice: 

end 
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